View a functional profile generated by Fama
Functional profile captures genetic potential of microbial community for biological processes of interest, just like taxonomic profile captures microbial diversity. This app generates visualization of a functional profile generated by the "Run Fama Read Profiling" app. Depending on the type and size of the input read libraries, Fama Read Profiling generates either normalized scores or raw read counts.
For comparison between functions and between samples, raw read counts are normalized by library size, target gene size and predicted average genome size in the sample. For normalization by average genome size, Fama employs MicrobeCensus tool. The normalization metric for single-read libraries is ERPKG (number of reads per kb of effective gene length per genome-equivalent):
ERPKG = (reads mapped to to gene) / (effective gene length in kb) / (genome equivalents),
where effective gene length = (actual gene length) + (read length) - 2 * (minimal alignment length) + 1,<\br> genome equivalents = (number of reads in library) / (average genome size)
The normalization metric for paired-end read libraries is EFPKG(number of fragments per kb of effective gene length per genome-equivalent), which is calculated similarly to ERPKG, but with fragment count instead of read count, and some minor differences in calculation of effective gene length.
If calculation of normalized scores is not possible (because of small number of reads etc.), only read counts and fragment counts are reported for single read and paired-end read libraries, respectively.
Functional profile object, output of the "Run Fama Read Profiling" app.
HTML page with two tabs: "Profile Statistics" and "Profile Heatmap". The "Profile Statistics" tab shows aggregating statistics for rows and columns. The "Profile Heatmap" tab contains interactive 2D heatmap of scores with samples in columns and functions in rows.
- Kazakov A, Novichkov P. Fama: a computational tool for comparative analysis of shotgun metagenomic data. Great Lakes Bioinformatics conference (poster presentation). 2019. , https://iseq.lbl.gov/mydocs/fama_glbio2019_poster.pdf
- Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nature Methods. 2015;12: 59-60. doi: 10.1038/nmeth.3176. Publication about third-party program used by Fama. , https://pubmed.ncbi.nlm.nih.gov/25402007/
- Nayfach S, Pollard KS. Average genome size estimation improves comparative metagenomics and sheds light on the functional ecology of the human microbiome. Genome Biology. 2015;16: 51. doi: 10.1186/s13059-015-0611-7. Publication about third-party program used by Fama. , https://www.ncbi.nlm.nih.gov/pmc/articles/pmid/25853934/
- Ondov B, Bergman NH et al. Interactive metagenomic visualization in a Web browser. BMC Bioinformatics. 2011;12: 385. doi: 10.1186/1471-2105-12-385. Publication about third-party program used by Fama. , https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3190407/
Module Commit: d9db15ea217e3be2aab65c356564a6d345b4f410