How to visualize MAGs
2
0
Entering edit mode
12 months ago

I have recently analyzed my metagenomic data using the metagenome-atlas. I have a set of MAGs of varying completeness/contamination and from different lineages and am looking for a way to visualize the coverage.

I have previously used BRIG before but not for coverage mapping, but I am aware that it does have a coverage mapping module. However, I am confused as to what my reference data would be.

I have also had some R packages recommended to me by colleagues (ggbio, circlize) but have no experience using these.

Any suggestions on what to use? Unfortunately I only have the atlas outputs because I do not have access to a server which can operate metagenome-atlas anymore.

next-gen metagenome R • 544 views
ADD COMMENT
0
Entering edit mode
12 months ago
Mensur Dlakic ★ 15k

If you don't have coverage information, that can be obtained on a per-contig basis by mapping the sequencing reads onto your metagenomic assembly.

Assuming you have some kind of embedding information for your bins - say coordinates from t-SNE or UMAP embedding - it becomes simply a matter of creating a scatter plot and coloring each dot by coverage values.

enter image description here

ADD COMMENT
0
Entering edit mode

Looks great - do you have any tool suggestions to do this? I currently don't have coverage stats.

ADD REPLY
0
Entering edit mode

There are many tools to calculate contig coverage. I will recommend calculate_read_coverage.py from Autometa.

Making a plot like I did above is fairly trivial once you have [X, Y, coverage] columns. It is a scatter plot of [X, Y] colored by coverage values, which can be done using many tools such as pandas.plot.

ADD REPLY
0
Entering edit mode
12 months ago

If you want to analyze the coverage of you MAG trought the metagenomes you should use AnviĀ“o :) .... and if you want to see very quickly you can use anvi snakemate pipeline for metagenome analysis

https://merenlab.org/2018/07/09/anvio-snakemake-workflows/#metagenomics-workflow

ADD COMMENT

Login before adding your answer.

Traffic: 2121 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6