Phylogenetic analysis
From the 234 sequences dataset, 26 HA partial or identical sequences were removed. A total of 208 complete HA gene sequences were aligned using the MUSCLE algorithm (Edgar, 2004). A circular phylogenetic tree was constructed by the Maximum-likelihood method with 1000 bootstrap replicates using the best-fit nucleotide substitution model in MEGA X (Kumar et al., 2018) and was edited with FigTree v1.4.4 software (Rambaut, 2018). The nucleotide substitution model was selected using the Akaike and Bayesian information criteria (AIC and BIC) (Kumar et al., 2018).