Coverage Calculation
1
6
Entering edit mode
9.9 years ago
Raygozak ★ 1.4k

HI, if i'm sending samples for sequencing at say 30X coverage of a given genome, what is the relation between the 30X,40X coverage metric and the minimum number of reads that i should expect to get back from the sequencing?

Thanks

next-gen coverage • 13k views
ADD COMMENT
9
Entering edit mode
9.9 years ago
Irsan ★ 7.6k

Let's say the read length is 100 bps/read and the genome you are sequencing is 3.000.000.000 bps (= human =3 Gbs) and the average coverage is 30x (calculated after alignment). The minimum number of bases that you should have in your raw data is

3.000.000.000 * 30= 90.000.000.000 bases.

Given you have a read length of 100 bps you should have

90.000.000.000 / 100 = 900.000.000 reads = 450.000.000 read pairs (if you did paired end sequencing)

However, not all your raw reads will map to the reference genome so the number will even be bigger, 450.000.000 is an estimate of the minimum

also see technote from illumina

When you are particularly interested in optimizing coverage for exome seq you might find this information useful

ADD COMMENT

Login before adding your answer.

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