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

Clipped reads not removed #11

Open
annecmg opened this issue Jun 9, 2020 · 0 comments
Open

Clipped reads not removed #11

annecmg opened this issue Jun 9, 2020 · 0 comments

Comments

@annecmg
Copy link

annecmg commented Jun 9, 2020

I was happy to find this tool which does exactly what I want: Remove reads with many clipped bases, but not in the beginning and end of contigs.
However, I find that too many reads remain in the mapping, I run the debug mode and here is one example of a read that should be BAD:

GOOD CHROM=GC00000017_11_p1:1..357 POS=149..299 CIGAR=89M62S L=0 R=0 | HL=0 SL=0 SR=62 HR=0 max=5)

The read has only 89M bases so I think POS should be 149..237 and the end of the contig (357) is not reached yet.

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