View Single Post
Old 12-18-2014, 01:50 PM   #2
Brian Bushnell
Super Moderator
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707

BBDuk will do that. For example, if you want to remove phiX reads, and reads with average quality below 10: ref=phix.fasta k=31 maq=10 in=reads.fq out=good.fq outm=bad.fq

Passing reads go to out, and failing reads go to outm. I'm not sure whether you can do this with Trimmomatic; the manual seems to indicate that you can save pairs where one fails and the other doesn't, but does not seem to mention a way to catch all filtered reads.
Brian Bushnell is offline   Reply With Quote