Here , we present a highly accurate probabilistic algorithm , which is an extension of the Genomic Next-generation Universal MAPper to accommodate bisulfite sequencing data ( GNUMAP-bs ) , that addresses the computational problems associated with aligning bisulfite sequencing data to a reference genome .