Question: How to visualize MAGs
0
gravatar for ellenphante
9 days ago by
ellenphante0 wrote:

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.

metagenome next-gen R • 82 views
ADD COMMENTlink modified 8 days ago by gamamiguelangel0 • written 9 days ago by ellenphante0
0
gravatar for Mensur Dlakic
9 days ago by
Mensur Dlakic7.2k
USA
Mensur Dlakic7.2k wrote:

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 COMMENTlink written 9 days ago by Mensur Dlakic7.2k

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

ADD REPLYlink written 8 days ago by ellenphante0

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 REPLYlink written 7 days ago by Mensur Dlakic7.2k
0
gravatar for gamamiguelangel
8 days ago by
gamamiguelangel0 wrote:

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 COMMENTlink written 8 days ago by gamamiguelangel0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1252 users visited in the last hour