A program that computes functional enrichment across groups of genomes..
This program computes functional enrichment across groups of genomes and returns a functional-enrichment-txt file.
Please also see anvi-display-functions which can both calculate functional enrichment, AND give you an interactive interface to display the distribution of functions.
You can use this program by combining genomes described through external-genomes, internal-genomes, and/or stored in a genomes-storage-db. In addition to sources for your genomes, you will need to provide a groups-txt file to declare which genome belongs to which group for enrichment analysis to consider.
How does it work?
Aggregate functions from all sources. Gene calls in each genome are tallied according to their functional annotations from the given annotation source.
Quantify the distribution of functions in each group of genomes. This information is then used by
anvi-script-enrichment-statsto fit a GLM to determine (1) the level that a particular functional annotation is unique to a single group and (2) the percent of genomes it appears in in each group. This produces a functional-enrichment-txt file.
You can use it with a single source of genomes:
You can get a tab-delimited matrix describing the occurrence (counts) of each function within each genome using the
Edit this file to update this information.
Are you aware of resources that may help users better understand the utility of this program? Please feel free to edit this file on GitHub. If you are not sure how to do that, find the
__resources__ tag in this file to see an example.