Practices for mapping reads to De Novo metagenomic assembly
0
0
Entering edit mode
2.9 years ago
doron • 0

Hi all,

I have performed de novo assembly for several metagenomic samples (using MegaHIT). The overall statistics (e.g. N50) seem OK. I am now interested in mapping the original reads to the contigs (to estimate coverage, which I need for binning).

However, it is not entirely clear what parameters should I use for the mapping, in terms of number of allowed indels, mismatches, longest indel etc. Are there any known practices? I was not able to find any.

Note: Olson et. al discuss the various sources of assembly errors, but I wasn't able to make it useful.

Thanks!

assembly metagenomic de novo reads megahit • 746 views
ADD COMMENT
0
Entering edit mode

Hi There,

I think the best thing would be to just use default parameters in mapping tool such as bowtie2 and then look at the alignment stats. You can use tools such as weeSAM that can produce alignment stats.

ADD REPLY
0
Entering edit mode

Thanks! However, I'm confident that different parameters will yield different alignment stats... I doubt that same default parameters may suit all alignment objectives. Does my concern make sense?

ADD REPLY
0
Entering edit mode

If you are too paranoid then you can use --sensitive (I think this is the default on bowtie2) or --very-sensitive parameters with bowtie2. However, the default parameters would still help you get an idea of the number of mapped reads, read depth and other assembly stats. Alternatively, you could use other tools such as QUAST to assess the assembly quality.

ADD REPLY

Login before adding your answer.

Traffic: 2269 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