Hi Mike,
Thanks for you reply! Learned a lot! Wii try!
I guess the last one also applies to qiime rescript get-silva-data, like this:
## get-silva-data
qiime rescript get-silva-data \
--p-version '138.1' \
--p-target 'SSURef_NR99' \
--p-include-species-labels \
--o-silva-sequences silva-138-99-seqs.qza \
--o-silva-taxonomy silva-138-99-tax.qza
## Dereplicate
qiime rescript dereplicate \
--i-sequences silva-138-99-seqs.qza \
--i-taxa silva-138-99-tax.qza \
--p-mode 'uniq' \
--p-threads 8 \
--o-dereplicated-sequences silva-138-99-seqs-derep.qza \
--o-dereplicated-taxa silva-138-99-tax-derep.qza
## extract-reads
qiime feature-classifier extract-reads \
--i-sequences silva-138-99-seqs-derep.qza \
--p-f-primer GAGTTTGATYMTGGCTCAG \ #28F-YM
--p-r-primer GCTGCCTCCCGTAGGAGT \ #388R
--p-n-jobs 8 \
--o-reads ilva-138-99-seqs-segments.qza
Many thanks!
Birong