wgbs_tools icon indicating copy to clipboard operation
wgbs_tools copied to clipboard

modbam reads and hardclipping

Open jts opened this issue 9 months ago • 0 comments

Hi,

Thanks for providing this tool, it is very useful.

I have been testing it on nanopore data and noticed a potential issue in the modbam parser. In the current implementation of MM tags reads that are hardclipped by the aligner are unreliable as it breaks the correspondence between SEQ and the MM tag (if interested there's a discussion here). As far as I can tell ont.cpp doesn't check for this condition so it is possible the .pat output is incorrect for hard clipped records.

Jared

jts avatar May 30 '24 17:05 jts