The emerging popularity of single-cell sequencing technologies has revealed unexpected heterogeneity of chromatin states and gene expression in seemingly homogeneous cell populations. Similarly, single-cell proteomics has proven capable of quantifying heterogeneity of proteomes amongst single cells, yet, due to their low stoichiometry, it remains limited in investigating the heterogeneity of protein post-translational modifications (PTMs). Here, we present a robust mass spectrometry-based method for the unbiased analysis of histone PTMs (hPTMs) at single-cell level, in order to investigate the complex relationship between the regulation of epigenetic marks and cell identity. Our method identified both single and combinatorial hPTMs (68 peptidoforms in total), including nearly all of the most studied hPTMs. We demonstrate technical reproducibility comparable to traditional ‘bulk’ experiments. Furthermore, our platform distinguishes/removes technical noise from true biological cell-to-cell variation of hPTM abundances. As a proof of concept, we treated cells with sodium butyrate, a histone deacetylase inhibitor, and demonstrated that our method can i) distinguish between treated and non-treated cells, ii) identify sub-populations of cells that respond differently to treatment, and iii) reveal differential co-regulation of hPTMs in the context of drug treatment. Overall, these experiments establish the applicability of our method to investigate chromatin heterogeneity at the single-cell resolution, which has important implications for understanding complex conditions like cancer and aging. In addition, this method provides new opportunities for analyzing the ‘histone code’ by utilizing cellular heterogeneity in covariation analysis.