The purpose of the present study was to examine the genetic variation in Salamandra infraimmaculata from different breeding site habitats using the Amplified Fragment Length Polymorphism (AFLP) method. The results of the dendogram from a hierarchical cluster analysis show that the grouping of S. infraimmaculata as cluster 5 differs from all the other clusters, including the St1 (Tel-Dan stream) population, which was the most predictable. Five Haplogroups (Hg) were characterized. The mean number of alleles per locus in each population (Ne) ranged from 10.566 (Sp1) to 2.720 (Po6). An average estimated heterozygosity (He) by population ranged from 0.100 (Po6) to 0.186 (St1). Population St1, a permanent breeding site where water was available all year round, exhibited the highest level of polymorphism, while population Po6, from the ephemeral breeding site, exhibited the lowest level of polymorphism. Gene flow between clusters showed that clusters 3 and 4 are sources of migrants and also receive gene flow, while clusters 1 and 2 may be a source of migrants but may not receive much gene flow. A phylogenetic analysis, based on clustering using Nei’s genetic distance, demonstrated that the Tel-Dan population is located on a separate branch within its sub-population. The conclusion of the present study shows that the genetic divergence among isolated populations is not correlated to distance but is affected by the variation of habitats.