This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Allele count of 2 for homoplasmic MT variants in VCF

I have noticed that homoplasmic variants in VCF files generated by GATK HC have allele count of 2. I am wondering if this is a common practice in generating VCF files for mitochondrial variants?

variant-calling gatk vcf

Please stop using bioinformatics as a tag unless your post is about the field of bioinformatics itself.

1 answer

The GATK (and afaik every caller that I know) treats every contig of the reference as having the same ploidy (so all haploid, all diploid, etc). You will also see diploid calls on X and Y chromosomes unless you call them separately. The mitochondria is simply another contig and is treated the same way.

The best practices (according to GATK) for processing the MT is here: https://gatk.broadinstitute.org/hc/en-us/articles/4403870837275-Mitochondrial-short-variant-discovery-SNVs-Indels

Thank you for the helpful answer!

Please accept the answer (green check mark) to provide closure to this thread.

Log in to answer this question.