Trying to filter asv's based on length

I would like to filter out asv’s below a certain length, because they seem to be sequencing artifacts. I am trying to use the very elegant solution posted here by @thermokarst
Filtering out ASVs from DADA2 based on length

However, nothing seems to be happening at the terminal window when I run this. The command I used was:

qiime feature-table filter-seqs \
    --i-data rep-seqs.qza \
    --m-metadata-file rep-seqs.qza \
    --p-where 'length(sequence) > 350’ \
    --o-filtered-data 350-seqs.qza

Am I using this correctly?

Looks like you have a grave mark (or a smart quote) instead of a single quote after the number “350”. Updated, it would look like this:

qiime feature-table filter-seqs \
    --i-data rep-seqs.qza \
    --m-metadata-file rep-seqs.qza \
    --p-where 'length(sequence) > 350' \
    --o-filtered-data 350-seqs.qza
1 Like

Brilliant!
Thank you so much @thermokarst.
It works like a charm now!
:smiley::rainbow:

This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.