This study aimed to explore geographic distribution and molecular characterization of cutaneous leishmaniasis (CL) species by amplifying two popular markers in kinetoplast DNA and internal transcribed spacer 1 loci by nested-PCR, and characterized by sequencing and phylogenetic analyses. Findings demonstrated that two species co-existed in the province: L. tropica (88.5%) and L. major (11.5%). All gender and age groups were equally infected, although males, 21-30 years old, exhibited a significantly higher infection. Sequencing and phylogenetic analyses of 34 randomly selected samples showed that L. tropica isolates exhibited some degree of heterogeneity. Both anthroponotic CL and zoonotic CL are present in south-eastern Iran with predominance of L. tropica species. Some level of heterogeneity was observed in L. tropica isolates which possibly reflects different colonies in the area. Implementation of diagnostic tools directly from clinical samples could be an important strategic approach for exploration of spatial distribution, molecular characterization and phylogenetic analyses.