Hello,
I would like to use filter-seqs-length-by-taxon
to discard archea and bacteria to generate a database consists of only eukaryota. To achieve this, I have downloaded silva 138 with the command below.
qiime rescript get-silva-data --p-version 138 --p-target SSURef_NR99 --p-include-species-labels --p-download-sequences --o-silva-sequences silva-138-sequences.qza --o-silva-taxonomy silva-138-taxonomy.qza
The problem I face here is, filter-seqs-length-by-taxon
does not accept silva-138-sequences.qza due to it is a FeatureData[RNASequence]. Expected artifact type is FeatureData[Sequence].
I have skipped cull-seqs
step, because I have thought it would take less time if I discarded sequences first. I think that causes this error, seems like cull-seqs
changes the artifact type after culling.
If I'm not mistaken, filter-seqs-length-by-taxon
should accept FeatureData[RNASequence] artifacts, since it just filters by length not the type of sequence.
Could you inform me where my thought process goes wrong?
And also, thank you all for such a plugin. It helps a lot.