A probabilistic generative model for quantification of DNA mudifications enables analysis of demethylation pathways