Viral sequences were aligned using the E-INS-i algorithm implemented in MAFFT version 7 (Katoh and Standley 2013). Ambiguously aligned regions were then removed using the TrimAl program (Capella-Gutierrez et al. 2009). The best-fit model (LG+Γ) of aa sequence evolution was estimated using MEGA version 7.0 (Kumar et al. 2016). Phylogenetic trees were then estimated using the maximum likelihood (ML) method implemented in PhyML version 3.0 (Guindon et al. 2010) with bootstrap support values calculated from 1,000 replicate trees. Bootstrap values >70 per cent were considered significant. Identities among nt and aa sequences were calculated using the MegAlign program implemented in the Lasergene software package v5.0 (DNAstar).