This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Filtering FASTQ Files

I am trying to filter paired end reads from a FASTQ file based on ones that match to a specific genome. I want to keep any reads that align to the genome. So far all examples of this I have found are to remove the matches, does anyone have any examples of how I can do the reverse?

Thank you.

fastq bwa

1 answer

You can use bbsplit.sh from BBMap suite with the genome of your interest to bin the reads or use bbduk.sh from same suite.

Log in to answer this question.