This is a test version of Biostars. For the public version, visit https://www.biostars.org.
filter_fasta.py not removing sequences from fastq based on read IDs

I'm attempting to use filter_fasta.py in macqiime to remove sequences from a fastq based on a .txt file of read IDs. 

filter_fasta.py -f my_reads.fastq -o filtered_reads.fq -s read_ids.txt -n

The input file has 366000 reads. The output file is 366000 reads. Nothing is being removed. I ensured the read IDs are actually represented in the fastq with grep. My read ID .txt file has only one ID per line. Any ideas what could be wrong?

sequencing sequence fastq qiime
• 2,494 views • link
updated by lakhujanivijay lakhujanivijay • written by fjs5035 • fjs5035

Hi fjs5035

I just added a hyperlink to the script you are referring to.

Additionally, could you add the outputs of following commands, it will you get quick answers

  1. output few lines from your id file

    head -n 10 read_ids.txt

  2. output few read headers

    grep "^@" my_reads.fastq | head -n 10

• 1 views • link
written by lakhujanivijay lakhujanivijay

There is also a very useful tool in bbMap for your requirement. In case your issue persists, give a try with filterbyname.sh in bbMap suite of tools available at link.

• 1 views • link
written by Jeffin Rockey Jeffin Rockey

Just a hunch, does your read_ids.txt contain the "@" at the start of the fastq identifiers?

• 1 views • link
written by cschu181 cschu181

It does. I take it from your comment that that shouldn't happen.

• 1 views • link
written by fjs5035 • fjs5035

0 answers

No answers yet.

Log in to answer this question.

More posts like this