All the sequence alignments of the various HAs and NAs were performed using Clustal 2.0. Repetitive, incomplete and mixed sequences were removed. For the purpose of convenience, one representative sequence with the generally longest length from each subtype was chosen for further description (Table S1 and S2).