Phylogenetic analysis

Before constructing the ML phylogenetic trees in this study, sequence saturation was evaluated using DAMBE software (51, 52). The analysis revealed that all sequences of ML trees exhibited little saturation. Phylogenetic analysis of the complete genome and VP1 genes was performed using the ML method available in MEGA version 11 (53) with 1,000 bootstrap replicates. Complete genome tree employing an optimal nucleotide substitution model of Tamura-Nei distances and a gamma distribution with invariant sites (G + I) and VP1 tree employing the GTR nucleotide substitution model and a gamma distribution of rate variation among sites (G + I) were recommended by the models function in the MEGA version 11. All resulting phylogenetic trees were visualized using iTOL (54).