This task uses PICRUSt2 : Phylogenetic Investigation of Communities by Reconstruction of Unobserved States(paper can be found here).It wraps a number of tools to generate functional predictions based on 16S rRNA gene sequencing data.
The input files should be a FASTA of amplicon sequences variants (ASVs; i.e. your representative sequences, not your raw reads, which is ASV-sequences.fasta generated using Q2FeatureInference task) and a qza table of the abundance of each ASV across each sample which is table.qza (generated using the same previous task ).
The key output files are:
- EC_metagenome_out : Folder containing unstratified EC number metagenome predictions.
(pred_metagenome_unstrat.tsv.gz), sequence table normalized by predicted 16S copy number abundances.
(seqtab_norm.tsv.gz), and the per-sample NSTI values weighted by the abundance of each ASV (weighted_nsti.tsv.gz).
- KO_metagenome_out : As EC_metagenome_out above, but for Kegg orthology metagenomes.
- pathways_out : Folder containing predicted pathway abundances and coverages per-sample, based on predicted EC number abundances.
