Fifty-seven genotypes from eight population of Satureja bachtiarica was evaluated using fifteen ISSR and eleven RAPD markers. DNA profiling using RAPD primers amplified 84 loci, among which 81 were polymorphic with an average of 7.36 polymorphic fragments per locus. Also, using RAPD markers maximum and minimum polymorphic bands observed for Semyrom (77.38 %) and Farsan (40.48 %) populations, respectively. Semyrom population recorded the highest unbiased expected heterozygosity (0.259) and Shannon's Indices (0.38). While, the lowest values of unbiased expected heterozygosity (0.172) and Shannon's Index (0.245) were recorded for Eghlid and Farsan populations, respectively. On the other hand, ISSR primers produced 136 bands, from which 134 were polymorphic with an average of 9.06 polymorphic fragments per primer (98.52 %). The ISSR markers evaluation revealed that maximum and minimum polymorphic bands observed for Semyrom (66.18 %) and Farsan (31.62 %), respectively. Shahrekorud population recorded the highest unbiased expected heterozygosity (0.211) and Shannon's Indices (0.301). While, the lowest value of unbiased expected heterozygosity (0.175) observed for Farsan and Yazd populations and the lowest Shannon's Index (0.191) recorded by Farsan population. The overall results of the study revealed that both ISSR and RAPD markers were effective for evaluation of genetic variation of S. bachtiarica.