Domestic donkey plays a key role as a draft animal in rural economy of Pakistan where its population is increasing every year. The complete mtDNA control region of forty randomly sampled donkeys was PCR- amplified and sequenced bi-directionally using specific primers. Distinct mtDNA haplotypes obtained in the current study (KY446001-KY446011) were subjected to haplotype (h) and nucleotide diversity (π) measures using DnaS as well as to phylogenetic, Network, and AMOVA analyses. There were a total 27 polymorphic sites present within 11 unique mtDNA haplotypes from the studied 40 animals from different regions. Neighbor-joining network and median-joining network both illustrated the splitting of all these haplotypes into two well-defined Nubian and Somali lineages, confirming African maternal origin of Pakistani domestic donkey. Diversity parameters h (0.967± 0.037) and π (0.02917± 0.00307) were found to reveal high levels of genetic diversity in Pakistani donkeys. AMOVA demonstrated only 1% of genetic differences between two mtDNA maternal lineages, pointing to lack of population substructure in Pakistani donkeys as is the case with worldwide domestic donkey population. Pakistani donkeys have African maternal origin and high levels of mtDNA diversity. High genetic diversity may be due to non-selective breeding and heteroplasmy. We herein provide the first report on mtDNA diversity of control region in Pakistani domestic donkey.