I am working with WGBS data and I want to extract the DMRs in a particular mutant.
I used bwameth to align the sequencing data to mm10 genome. After that, I sorted the alignment file and used methyldackel to see the methylation status. I used the metilene tool to extract DMRs but this doesn’t work. I get an error message like this
" null device 1 [WARNING] Tue Aug 9, 16:7:8, 2022 Input files need to be SORTED, i.e., use “bedtools sort -i file >file.sorted”. I sorted the methyldackel output files using Sortbed and ran metilene again. Same error message. I read from the manual for Metilene, that the input file has to be in a certain format and the input file can be prepared using bedtools unionbedg. But this tool is not present in usegalaxy.eu or .org. Can anyone help me to prepare the input files for metilene?
Thanks for asking! Yes I did, by removing the last two columns and first row (description) of the Methyldackel output. It then did the trick and Metilene took them with no problem! The weird thing is that I didn’t have this problem at all a year ago when I used the exact same files and process (using it for teaching a course).
I am using metilene in use galaxy.eu and ran into the same issue . The bed file output of Methyldackel is not getting executed. I have also tried sorting the bed files. Please can you suggest.
.