Question: Good coverage value for de novo assembly with NanoPore Reads?
0
gravatar for marianaer98
17 months ago by
marianaer980 wrote:

Hello, I am trying to do a de-novo assembly of a turtle genome of about 2.15 gbases with Nanopore Minion reads. I have read about Canu as a good tool for this purpose, but it seems it yields good a good assembly with reads of at least 20x coverage, and the coverage of my reads is only about 2.5x. Is this value too low for de novo assembly with Nanopore reads? What would be a good coverage value?

Thank you everyone

ont sequencing assembly • 965 views
ADD COMMENTlink modified 16 months ago by colindaven2.3k • written 17 months ago by marianaer980
3
gravatar for colindaven
16 months ago by
colindaven2.3k
Hannover Medical School
colindaven2.3k wrote:

Yes, 2.5X is too low for de novo.

For structural variant detection vs a reference genome about 10X is recommended.

For de novo I would suggest 30-40X, whereby the longest reads are the most important.

Sorry, but 2.5 X is likely only maybe useful for somehow scaffolding Illumina contigs if you have them (difficult, poor results expected).

You can print the read distribution eg by using stats.sh from bbmap so we can have a look if you like.

The tool Miniasm or wtdbg2 https://github.com/ruanjue/wtdbg2 might give you a very very sketchy assembly but I doubt it.

ADD COMMENTlink written 16 months ago by colindaven2.3k

Hi,

For de novo I would suggest 30-40X, whereby the longest reads are the most important.

Do you have a source for this ? thanks.

ADD REPLYlink written 10 months ago by Picasa530

Pretty sure these figures come from work from the group around Adam Philippy who wrote Canu. Makes sense though, 2.5X is no way near enough to actually correct long reads from, and self corrected reads are needed before the assembly overlap step. They have published theoretical curves on effect of read coverage on NG50 , check the Canu and alternatively Flye and wtdgb2 papers.

ADD REPLYlink written 10 months ago by colindaven2.3k
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: 1491 users visited in the last hour