Bisulfite analysis was performed using Bismark-bowtie 2 (https://www.bioinformatics.babraham.ac.uk/projects/bismark/). We used the default Bowtie2 implementation of Bismark with the specifications that only uniquely mapped reads should be aligned. All alignments were performed with high stringency allowing for only one base mismatch (n = 1). We also clearly specified that no discordant pairs of the pair-end reads should be aligned. DNA methylation in the CG, CHG, and CHH contexts was calculated by dividing the total number of aligned methylated reads by the total number of methylated plus un-methylated reads.
Free full text: Click here