ZeroDivisionError about bsmapz

[praxsubba]

I had the same issue. Just download the methratio.py script from GitHub and replace it with the one you used previously. The issue arises because of the samtools -X option in the old script. Hope that helps!

[praxsubba]

You're welcome @NaKent. You should check your methratio.py output files to see if they have any data at all. If they do, then, it should work. I recommend following the methylkit tutorial: https://bioconductor.org/packages/devel/bioc/vignettes/methylKit/inst/doc/methylKit.html

[NaKent]

Thank you very much!! The problem was solved.
If you don't mind, I want to ask unrelated question. (Sorry, I'm stranger to this field.)
I tried to run DMR with MethylKit, but the format of methylation ratio file on this BSMAPz seems to be different from original BSMAP's one.
So, THIS PAGE is no longer helpful for me. How can I do DMR with the result of methratio.py?

[NaKent]

Thank you for your help.
After I learned the page and awk, I achieved analyzing the data with R :)