Fig. 3

Phylogenetic analysis of HPV16 E6 (a) and HPV18 E6 (b) gene regions was conducted using the maximum-likelihood method based on the Kimura 2-parameter model with bootstrap resampling (1000 replicates) by the MEGA 6 package [37]. Numbers above the branches indicate the bootstrap values. Ninety-six different nucleotide patterns of studied sequences were indicated by black triangles (GenBank accession numbers OP572427 through OP572474 for HPV16 and OP572475 through OP572522 for HPV18). The accession number of reference sequences of each sublineage used for phylogenetic analysis in this study was indicated by white triangles