Hi Dhruv,

I regularly use CREST and RDP Classifier for generating abundance profiles for 16S RNA (at Phylum, Class, Order, Family and Genus level). You can use MetaPhlAn, TAXAassign (my own software), and MEGAN for WGS. To help people understand how to process metagenomic datasets, I have put up two tutorials (using an example dataset) on my website and discuss how you can go from Raw Reads to Abundance Tables + afterwards, how you can then perform multivariate statistical analysis using the abundance table generated from these classifiers against environmental parameters that you recorded. The links are given below, if you follow the same structure as I have you can simply copy/paste the commands from my tutorials.

• 16S rRNA Amplicons Processing Workflow for Illumina: http://userweb.eng.gla.ac.uk/umer.ijaz/bioinformatics/Illumina_workflow.html

• WGS Short-read based Workflow (MetaPhlAn pipelines): http://userweb.eng.gla.ac.uk/umer.ijaz/bioinformatics/Metagenomics_tutorial.html

• WGS Contig-based Workflow (CONCOCT): https://github.com/BinPro/CONCOCT

• Generating abundance tables from RDP and CREST: http://userweb.eng.gla.ac.uk/umer.ijaz/bioinformatics/oneliners.html?#PHYLOMAP

Best Wishes, Umer