A number of sequences (Supplementary Table 1) of serotypes O and A FMDV circulating in Bangladesh in 2012–13 (KC795948–KJ175182) and reference sequences (NC011450 and NC004004) from NCBI database were used to predict the amino acid sequence variability of VP1 region. EBI, EMBOSS Bioinformatics Tools (Transeq) were used to translate the nucleic acid sequences into respective amino acid sequences. VP1 amino acid sequences and reference sequences were aligned by ClustalX (1.81) and analyzed using SeqMan II (Lasergene 8.0; DNAStar Inc., WI, USA) for discerning amino acids sequence variability. Protein variability server (PVS) was used to calculate protein variability index using Wu–Kabat variability coefficient.19 (link) The variability coefficient is computed using the following formula: variability = N*k/n, where, N is the number of sequences in the alignment, k is the number of different amino acids at a given position, and n is the times that the most common amino acid at that position is present.20 (link)