Using Raman microscopy, we investigated epithelial cervical cells collected from 96 women with squamous cell carcinoma (SCC) or belonging to groups I, IIa, IIID-1 and IIID-2 according to Munich III classification (IIID-1 and IIID-2 corresponding to Bethesda LSIL and HSIL groups, respectively). All women were tested for human papillomavirus (HPV) infection using PCR. Subcellular resolution of Raman microscopy enabled to understand phenotypic differences in a heterogeneous population of cervical cells in the following groups: I/HPV−, IIa/HPV−, IIa/HPV−, LSIL/HPV−, LSIL/HPV+, HSIL/HPV−, HSIL/HPV+ and cancer cells (SCC/HPV+). We showed for the first time that the glycogen content in the cytoplasm decreased with the nucleus size of cervical cells in all studied groups apart from the cancer group. For the subpopulation of large-nucleus cells HPV infection resulted in considerable glycogen depletion compared to HPV negative cells in IIa, LSIL (for both statistical significance, ca. 45%) and HSIL (trend, 37%) groups. We hypothesize that accelerated glycogenolysis in large-nucleus cells may be associated with the increased protein metabolism for HPV positive cells. Our work underlines unique capabilities of Raman microscopy in single cell studies and demonstrate potential of Raman-based methods in HPV diagnostics.