Fig. 3

The evolutionary history was inferred by using the Maximum Likelihood method and Kimura 2-parameter model [1]. The tree with the highest log likelihood (-1872.22) is shown. 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 Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach and then selecting the topology with a superior log-likelihood value. A discrete Gamma distribution was used to model evolutionary rate differences among sites (5 categories (+G, parameter = 0.7700). The rate variation model allowed for some sites to be evolutionarily invariable ([+I], 43.64% sites). The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. This analysis involved 15 nucleotide sequences. There were a total of 873 positions in the final dataset. Evolutionary analyses were conducted in MEGA X [2]. Red and green dots indicate haplotypes identified in this study; squares with different colors represent reference haplotypes extracted from GenBank