Is there a simple command that returns read counts to keep track of what's happening while going through the QIIME2 pipeline for 16s rRNA data?

User Support
3

Hi @Nicholas_Bokulich,
Thanks for the response! That's what I figured but wanted check just in case since generating the feature-table for the classifier reference reads was taking a (little) bit of time (each time I do it) and all I really needed was the read count and sequence lengths information.

Interesting note about q2-quality-control. I will check it out. It looks like there are mock communities that I can test against. Awesome!

So far I have,
qiime feature-classifier extract-reads
--i-sequences SILVA-138-SSURef-Full-Seqs.qza
--p-f-primer GTGYCAGCMGCCGCGGTAA
--p-r-primer GGACTACNVGGGTWTCTAAT
--o-reads ref_seqs_SILVA_138

I was playing around with the idea of adding --p-identity 0.9 (due to this) although I don't see any reason to.

I would really appreciate your general thoughts about it. Thanks again for all the great help and support!

1 Like