^aThe complete mtDNA genomes of Le131, Le834 and Le1143 contained no private non-synonymous and mt-tRNA variants and was not included in the table. ^bThe uniqueness of each variant was searched according to the described strategy [26] on 10 August, 2012 (e.g., both “G3635A mtDNA” and “3635G>A mtDNA” were queried). ^cThe haplogroup-specific variant was determined according to the available global mtDNA phylogenetic tree at the Phylotree (mtDNA tree Build 14, 5 Apr 2012). The haplogroup status as it defined in that tree was indicated in round brackets. ^d The variant frequency was calculated as the number of occurrences of each variant in 15,859 complete or near complete mtDNA sequences that were summarized by MitoTool project [24]. ^eThe conservation index which was defined by Ruiz-Pesini et al. [30], was calculated by using MitoTool [24] concerning 43 primate species. A CI value of 0.419 meant 41.9% of 43 primate species share the same allele with the revised Cambridge Reference Sequence (rCRS) [18] (GenBank accession number NC_012920). ^fThe pathogenicity score was referred to the recent study of Pereira et al. [27]. The score ranges from 0 to 1. Higher pathogenicity score is associated with greater possibility that the variant is pathogenic. ^gAll samples contain m.3635G>A.