Functional Metagenome & Pathway Analysis
For functional analysis, MOCAT2 was used to cluster high quality reads into reference gene catalogs which is then annotated to Evolutionary genealogy of genes: non-supervised orthologus group (eggnog) database46 . The gene abundance profile was calculated for each sample by mapping the raw reads back to the reference gene catalogue. Functional mapping and Analysis Pipeline (FMAP)47 , was used to compare the various time points for alteration in pathways, using the Kruskal-Wallis rank-sum test.