Tagged with #beagleoutputtovcf
1 documentation article | 0 announcements | 1 forum discussion


No posts found with the requested search criteria.
Comments (15)

I'm having a problem running BeagleOutputToVCF, showing an error message:

##### ERROR stack trace
java.lang.IllegalStateException: Allele in genotype AT not in the variant context [A*, G]
        at org.broadinstitute.sting.utils.variantcontext.VariantContext.validateGenotypes(VariantContext.java:1197)
        at org.broadinstitute.sting.utils.variantcontext.VariantContext.validate(VariantContext.java:1137)
...

Is this because my Beagle output file include INDELs?

I'm using GenomeAnalysisTK-2.2-10-gbbafb72, and the command line is like:

> java -Xmx4g -jar /raid/software/src/GenomeAnalysisTK-2.3-5-g49ed93c/GenomeAnalysisTK.jar \
-R hg19_lite.fa \
-T BeagleOutputToVCF \
-V chr1.01.vcf \
-beagleR2:BEAGLE beagle/chr1.01.r2 \
-beaglePhased:BEAGLE chr1.01.phased \
-beagleProbs:BEAGLE chr1.01.gprobs \
-o beagle_vcf/chr1.01.vcf \
-L 1:1-10000000

Thanks, Hiroko Matsui