Background: Molecular diversity of Leishmania major and its morphological changes have become a controversial issue among researchers. Some aspects of polymorphic shapes of amastigotes in clinical manifestations along with molecular variation were evaluated among suspected patients of some exceptional zoonotic cutaneous leishmaniasis locations in Northern Khuzestan, Southwestern Iran. Methods: Suspected patients (n = 165) were sampled in zoonotic cutaneous leishmaniasis foci over two consecutive years during 2012-2014. Prepared smears were stained, scaled and measured by ocular micrometer. DNA was extracted from smears ITS-rDNA and Cytochrome b (Cyt b) markers were amplified, and PCR products were digested by BsuR1 restriction enzyme. Then the RFLP and sequencing were employed. Results: Only L. major was identified in patients containing regular amastigotes' shapes (oval or round) with a size of 2-4 µm in each of classical wet, dry, mixed lesions. Meanwhile, irregular shapes (spindle, pear, or cigarette) were observed separately in non-classical wet lesions with more than 4 µm. Interestingly, a few amastigotes with an external flagellum were observed in some lesions. All sequenced ITS-rDNA and Cyt b genes of L. major did not show any molecular variation (&chi 2 P > 0.05), including only one common haplotype (GenBank access no. EF413075). Conclusion: Findings proved that unlike other endemic foci, there is not a meaningful correlation between phenotypic and genotypic features of L. major isolates. This study is considered as the first comprehensive report to incriminate morphometric shapes of L. major amastigotes, which enhances our knowledge concerning their relevance with various clinical appearances and genotypic traits.