substitution model, HKY+Gamma (Gamma =0.1869), was selected
by MEGA 5 [22] and applied to neighbor- joining tree
construction in MEGA 5. In addition, RAxML-HPC [23] was
used to perform phylogenetic analysis using maximum-likelihood
methods. To evaluate the robustness of the internal branches of
the ML tree, 500 bootstrap replications (BS) were run under the
GTR+ Gamma model. Sequence divergences between species
were calculated by MEGA 5 using