shira
1
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
shira
3
Brilliant!
Thank you so much @thermokarst.
It works like a charm now!
system
(system)
Closed
4
This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.