Runing MOABS on Bismark BAM files

[seqlad]

Hello,

I have BAM files generated with Bismark 0.19.0. I wish to detect differentially methylated Cytosines (CpG or non-CpG context). I have tried to run the MOABS mcall tool on those files (As suggested in the MOABS documentation) I get an error message:

program is parsing file according to BAM foramt
XR:Z or ZR:Z:, or ZS:Z: field not found, please enable -R option in bsmapXR:Z or ZR:Z:, or ZS:Z: field not found, please enable -R option in bsmap

Has anyone ever tried to run MOABS on bismark BAM files ?
Any idea as to how solve that problem would be greatly appreciated!
Also if you know some other tool that I could use to find DMCs from Bismark file please let me know!

Thanks