samclip icon indicating copy to clipboard operation
samclip copied to clipboard

zero clipping at the beginning of the read

Open avilella opened this issue 5 years ago • 3 comments

I googled for a tool like this and I was glad to find it. Thanks @tseemann !

Could an option be added so that --start_max N limits the reads that are accepted if they have more than N clipped bases at the beginning of the read? E.g. for cases like:

4S96M - would be filtered out with --start_max 1 --max 5 96M4S - would not be filtered out with --start_max 1 --max 5

avilella avatar Jan 03 '19 10:01 avilella

@avilella am i correct in that you would like separate --max for "start" and "end" ? out of interest, what is the use-case for this behaviour?

tseemann avatar Sep 15 '19 23:09 tseemann

Yes, a max for start and one for end.

Thanks in advance

On Mon, 16 Sep 2019, 01:48 Torsten Seemann, [email protected] wrote:

@avilella https://github.com/avilella am i correct in that you would like separate --max for "start" and "end" ?

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/tseemann/samclip/issues/8?email_source=notifications&email_token=AABGSNZMSLZEHJW4GQVFWCDQJ3CTLA5CNFSM4GM7RXOKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD6X3TJY#issuecomment-531610023, or mute the thread https://github.com/notifications/unsubscribe-auth/AABGSNZCG65CBSVYXIN3TIDQJ3CTLANCNFSM4GM7RXOA .

avilella avatar Sep 16 '19 06:09 avilella

@avilella what does "start" mean in a SAM file? I think alignments are always begin..end with begin<end and revcom in a flag?

tseemann avatar Oct 16 '19 22:10 tseemann