Ben (previous postdoc in the lab) recently discovered that MethylKit automatically filters data based on coverage when the files are first read into R. We usually destrand our data (i.e. merge the reads from the positive and negative DNA strand to increase coverage at each CpG site) after we have read the data in and … Continue reading Ensuring correct coverage in MethylKit →