Table 1

List of overlapping primer pairs used for mtDNA sequencing

Gene/region	Forward primer	Reverse primer	mtDNA region	PCR product size
ND1	ND1F1 TACTTCACAAAGCGCCTTCC	ND1R1 ATGAAGAATAGGGCGAAGGG	3150–3980	830bP
	ND1F2 TGGCTCCTTTAACCTCTCCA	ND1R2 AAGGATTATGGATGCGGTTG	3777–4660	883bP
ND4L	ND4LF TCTCCATCTATTGATGAGGGTCT	ND4LR AATTAGGCTGTGGGTGGTTG	9967–10 859	892bP
ND4	ND4F1 GCCATACTAGTCTTTGCCGC	ND4R1 TTGAGAATGAGTGTGAGGCG	10 653–11 511	858bP
	ND4F2 TCACTCTCACTGCCCAAGAA	ND4R2 GGAGAATGGGGGATAGGTGT	11 314–12 096	782bP
	ND4F3 TATCACTCTCCTACTTACAG	ND4R2 AGAAGGTTATAATTCCTACG	11 929–12 793	864bP
ND6	ND6F GCATAATTAAACTTTACTTC	ND6R AGAATATTGAGGCGCCATTG	14 081–15 017	936bP
D-loop	HSV I F1 TCATTGGACAAGTAGCATCC	HSV I R1 GAGTGGTTAATAGGGTGATAG	15 792–16 400& 1–31	809bP
	HSV II F2 GAGTGGTTAATAGGGTGATAG	HSV II R2 AGGCTAAGCGTTTTGAGCTG	10–794	784bP
M13 5′ Tag	M13F GTAAAACGACGGCCAGT	M13R CAGGAAACAGCTATGACC	–	–

Six primer pairs were used to sequence ND1, ND4, ND4L and ND6 genes; two primer pairs were also used to sequence the hypervariable region of the mtDNA. M13 generic tag sequences were engineered into the 5′ end of each primer for sequencing reactions robustness.