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