Bug Bulletin: The recent 3.2 release fixes many issues. If you run into a problem, please try the latest version before posting a bug report, as your problem may already have been solved.

Error in BeagleOutputToVCF

hirokomatsuihirokomatsui Posts: 8Member
edited January 2013 in Ask the GATK team

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

Post edited by Geraldine_VdAuwera on

Answers

Sign In or Register to comment.