Identifying differentially methylated cytosines in bisulfite sequencing data via a Bayesian regression model