We used Schmutzi (35) to determine the endogenous consensus mtDNA sequences. Reads were mapped to the Cambridge reference sequence and filtered for MAPQ ≥ 30. Haploid variants were called using the endoCaller program implemented in Schmutzi (35), and only variants with a posterior probability exceeding 50 on the PHRED scale (probability of error: 1/100,000) were retained. We then used Haplogrep 2.2 (36) to determine the mtDNA haplogroups, specifying PhyloTree (build 17) as the reference phylogeny (37).
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this article to respond.