“…Several nuclear and mitochondrial DNA markers (cytochrome b, COII, gp63, 18S-rRNA, mini-exon, HSP-70, and ribosomal internal transcribed spacer regions (ITS-rDNA)) as well as different molecular techniques such as polymerase chain reaction (PCR) [9], multilocus sequence typing (MLST) [10], multilocus enzyme electrophoresis (MLEE) [5,11], random amplified polymorphism (RAPD) [12], amplified fragment length polymorphisms (AFLP) [13] and microsatellites [14], have been used to evaluate the inter and intraspecific genetic variability of different Leishmania species (L. tropica, L. major, L. donovani, L. braziliensis, L. peruviana, L. guyanensis, and L. panamensis) [13][14][15][16][17][18][19]. However, the identification of such genetic variants across the whole genome has become possible only with the arrival of DNA-seq.…”