Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

modbam reads and hardclipping #75

Open
jts opened this issue May 30, 2024 · 0 comments
Open

modbam reads and hardclipping #75

jts opened this issue May 30, 2024 · 0 comments

Comments

@jts
Copy link

jts commented May 30, 2024

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant