I have analyzed the methylation status of biological samples (promoter region of target genes) using bisulfite sequencing. I have a small N and am interested in the effect nutrition has on the methylation status of these genes, therefore the differences are probably subtle. I have browsed through previous Q&A on this topic, and everybody recommends a different method. I was hoping someone could shed some light on this, so I can make an appropriate decision. My data is currently in binary (1=meth, 0=non methy) for every CpG found, per all samples. Thanks!