This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Plotting MBD-seq results

Does anyone know of any programs or R packages to plot MBD-seq methylation status results like this? Thanks!

I'm using the R package MEDIPS for my analysis.

enter image description here

methylation mbd-seq

1 answer

The deepTools package can do things like this

https://deeptools.readthedocs.io/en/develop/content/tools/plotProfile.html

You would prepare coverage files in the form of bigwigs with computeMatrix

https://deeptools.readthedocs.io/en/develop/content/tools/computeMatrix.html\

Converting aligned reads to bigwig is generally done from things like macs, which automatically makes wig files, or bedtools genomecov which creates a bedgraph file from a inputted bam, and then run bedGraphToBigWig

Very helpful, thanks! I wonder if you'd also know what package to use for performing GO annotations on MEDIPS results? I'd like to do GO annotations so I can visualize the results with clusterProfiler.

Log in to answer this question.