seqtk icon indicating copy to clipboard operation
seqtk copied to clipboard

trimfq to drop seq under minlen and output cutpoints

Open a1an77 opened this issue 8 years ago • 1 comments

Hi, I needed to change the behaviour of the trimming function to drop reads shorter than minlen after trimming instead of providing low quality bases and I added an option to output the cut points (-1 if discarded) in a separate file. If you like the changes to be included I can open a pull request. Also added a Q option to set the threshold with a Phred score

a1an77 avatar Dec 19 '17 11:12 a1an77

I think you can just seq trimfq ... | seqtk seq -L minlen ...

Might also be able to use seqtk seq -q PHRED .... to do the quality trimming.

tseemann avatar May 27 '18 07:05 tseemann