This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!
I've run the GATK best practices pipeline up through VQSR and have the recalibrated variants in a VCF file. Because I'm analyzing pedigree samples, I'm now attempting to run CalculateGenotypePosteriors. When I do this, I get the following error:
ERROR MESSAGE: Variant does not contain the same number of MLE allele counts as alternate alleles for record at 1:768589
However, when I looked at that variant in the recalibrated VCF file, I see that MLEAC=1, AC=1, and I can confirm that there is only one sample that is heterozygous, so both AC and MLEAC are correct.
Might this be a bug?