Taxonomic profiles were generated using mOTUs (v2.0). Furthermore, reads were assembled into contigs, which were used for gene predictions using Prokka , and subsequently annotated using the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologous categories. The KEGG orthologous (KO) genes were identified, and their coverages based on number of reads mapping to individual KOs were used for subsequently analyses.