We use data from whole genome or reduced representation bisulfite sequencing to test for differential methylation between experimental groups. Our workflow relies on the nf-core methylseq pipeline to perform the alignment and methylation calling. This is followed by an in-house pipeline, which tests for differential methylation at the level of individual CpGs and larger genomic regions.