Question: I need to get the fold-coverage
0
gravatar for kamel
2.6 years ago by
kamel10
kamel10 wrote:

Hello by what tool I know the genome coverage (X) of assembly contigs. Thank you

assembly • 857 views
ADD COMMENTlink modified 2.6 years ago by dungnv1140 • written 2.6 years ago by kamel10

Without knowing more information, see this post.

ADD REPLYlink written 2.6 years ago by st.ph.n2.4k

Thank you for your reply. but I can't use the calculated theoretical C = R x L / G, because I have the reads length 32-250 (so it is not a single value) and I can't use bedtools because I want to know the fold-covrage from assembly of Novo result (not from mapping).

ADD REPLYlink written 2.6 years ago by kamel10

If you want to know the coverage of your denovo assembly you could try mapping the reads to your assembly and then calculating coverage as suggested by the following answers. A not-so-correct way is to see how much of data you have (in total bases) in your assembly and use this as the denominator to divide the total data you have (i.e in place of G in the above formula). This post might help - Counting Number Of Bases In A Fastq File

But you would obviously be over-expecting coverage since there would be missing data due to a fragmented assembly. So you definitely need to map back your reads to your contigs (keeping these as a reference) to understand how much of coverage you might have.

ADD REPLYlink written 2.6 years ago by Rohit1.3k
2
gravatar for dago
2.6 years ago by
dago2.5k
Germany
dago2.5k wrote:

Googe is full of answers:

Tools To Calculate Average Coverage For A Bam File?

http://bedtools.readthedocs.io/en/latest/content/tools/genomecov.html

http://seqanswers.com/forums/showthread.php?t=17725

ADD COMMENTlink written 2.6 years ago by dago2.5k

I want the fold-coverage from contigs (fatsa file) not from mapping (bam file)

ADD REPLYlink written 2.6 years ago by kamel10

Do you have a reference genome available? If not, what "reference" (would that be the de novo assembly you seem to refer to in the post above) are you planning to use to calculate fold coverage?

ADD REPLYlink written 2.6 years ago by genomax67k

I made a de novo assembly of a bacterial strain and I want to know their genome coverage from fasta file (contigs) ie how many times the sequenced reads but not according to a reference

ADD REPLYlink modified 2.6 years ago • written 2.6 years ago by kamel10
1

In that case your contigs are going to be the "reference". You would need to align your reads to that "reference" and then calculate coverage using one of the methods @dago noted above.

ADD REPLYlink written 2.6 years ago by genomax67k
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: 1707 users visited in the last hour