Gene sequencing was performed using CLC Main Workbench 7.9.1. Alignment was conducted using human and animal coronavirus sequences registered in Global Initiative on Sharing All Influenza Data (GISAID) and NCBI GenBank. The phylogenetic tree was analyzed using MEGA6 with the neighbor-joining method, maximum composite likelihood-parameter distance matrix, and bootstrap values of 1,000 replicates.