1. a more detailed protocol can be found here:
https://github.com/nioo-knaw/orochi/blob/0.1/Snakefile#L500
2. Please find in the zenodo link the detals of the snakemake and the extraction of 16S reads from the metagenome.
https://zenodo.org/record/3405565#.YGWXxNztZhE
Also see below the description for the sup material:
4.Taxonomic analysis of 16S rRNA reads extracted from metagenome Ribosomal RNA reads were extracted from the quality and contamination filtered reads using a kmer strategy in BBDuk from the BBMap tool suite (54). These rRNA reads were aligned to SILVA 119 (55) reference set provided by Qiime (56) using the usearch_global algorithm implemented in VSEARCH (57). The alignment in usearch format (uc) was converted to a BIOM file with the biom 2.1.5 package (58) and the from-uc flag. All steps where implemented in a Snakemake workflow (59). The OTU table was filtered using QIIME (1.9.1) custom scripts (60).
I have asked Mattias de Hollander to give more details if possible. Also if you have more questions, do not hesitate in contacting me.
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this
article to respond.