The aim of this study was to identify and associate the sequence variations of human Papillomavirus 16 (HPV16) genes from women who live in two different areas of Romania and associate them with malignant progression. One hundred twenty-four HPV16-positive cervical isolates were collected, and the E2, E4, E5, E6 and E7 viral genes were sequenced. Two new missense mutations in the E6 gene (C279G and A305C) were found (together or alone, in association with other mutations) in 44 of 124 cases. The most frequently simultaneously mutated genes were E4/E2 hinge, E5 and E6 (p = 0.0004) in squamous cell carcinoma (SCC) samples. Also, for SCC patients, the best-correlated mutation patterns were obtained for E4/E2 hinge-E5 (r = 0.7984; p < 0.0001). No sample was found to have all of the investigated viral genes concurrently mutated. Phylogenetic analysis was performed to characterize the viral variants. Similar results were found for SCC and cervical intraepithelial neoplasia III (CINIII) cases. After all of the target gene sequences were assembled, all patients were found to be infected with viruses of the HPV16- European-German (EG) lineage, and two clusters were identified, the first (55/96 variants) from Moldavia and the second (41/96 variants) from Bucharest. The distinct cluster derived from EG in Moldavia could partially explain the increased frequency of SCC in this area. This study has generated a comprehensive set of sequence variation data on HPV16 circulating in Romania to join the existing data and highlight the important role of HPV16 variants during cervical carcinogenesis.