SUMMARY

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.

INPUT

Functional profile object, output of the "Run Fama Read Profiling" app.

OUTPUT

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.

Additional resources

• More info about Fama
• Example narrative

Team members who implemented App in KBase: Alexey Kazakov.For questions, please contact us.