This network was redrawn manually from the output of median joining/reduced network obtained using NETWORK program (version 4.1) (www.fluxus-engineering.com). Each sample represented on the diagram has been sequenced for the HVS-I region and genotyped for the coding region mutations that are indicated. Circle sizes are proportional to the number of mtDNAs with that haplotype. Recurrent mutations are underlined. Coalescent times were calculated by a calibration method described elsewhere [33]. 16182C, 16183C and 16519 polymorphisms were omitted. Suffixes A, C, G, and T indicate transversions.