Sequencing and phylogenetic analysis
Amplified DNA fragments of 1,210bp were used forthe dideoxynucleotide sequencing reactions. PCR-amplified products were analyzed on an ABI Prism 3700 automated DNA sequencer (Applied Biosystems, Foster City, California). Low quality sequences were trimmed with Chromas version 2.6.2 for windows licensed by Technelysium Pty Ltd 1998-2016 (www.technelysium.com.au). The sequences were combined with those currently available from GenBank (as shown in Table 2). Pairwise and multiple genomic alignments were done with Clustal W13 alignment programs. The evolutionary history was inferred by using the Maximum Likelihood (ML) method and Kimura 2-parameter model14 . The percentage of trees in which the associated taxa clustered together is shown next to the branches. Initial tree(s) for the heuristic search were obtained automatically by applying ML algorithm with a GTR+G model, and 1000 bootstrap resampling. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. This analysis involved 55 nucleotide sequences with 556 positions in the final dataset. Evolutionary analyses were conducted in iqtree15 . Sequences were submitted to the GenBank and accession numbers MW916908-MW916918 were assigned for 11 isolates in this present study. Information regarding the study sequences can be accessed in supplementary table 1.