estimate genome coverage
2
0
Entering edit mode
5.2 years ago
Ric ▴ 360

Hi, I ran KmerGenie and I got the best Kmer and the estimated genome size. How do I get the estimated genome coverage?

Thank you in advance.

Best wishes,

Mic

masurca kmergenie Assembly • 2.9k views
ADD COMMENT
1
Entering edit mode
5.2 years ago
Medhat 9.0k

you can use bedtools like this calculate coverage


Or GATK

ADD COMMENT
1
Entering edit mode
5.2 years ago
dbrowne.up ▴ 70

Or more broadly and easily... C = R x L / G

Where C is physical coverage, R is the total number of reads, L is (average) read length, and G is the genome size.

So for example, if you have a 100 Mbp genome and an Illumina library with 10 M reads of length 150 bp, the coverage of the genome is 15X (10,000,000 * 150 bp / 100,000,000 bp).

ADD COMMENT

Login before adding your answer.

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