Formula For Converting Sequence To Physical Coverage
1
0
Entering edit mode
10.1 years ago
Dan ▴ 520

If I have reads of length x, inserrt size of y and genome size z, how do I convert from sequence to physical coverage?

sequencing • 5.3k views
2
Entering edit mode

maybe you miss n the number of reads?

0
Entering edit mode

n is implied by the coverage

1
Entering edit mode
10.1 years ago

The "normal" definition of coverage would be C = N * L / G

Now if you want to extrapolate from that you could use the insert size for L but then you would need to also divide N by 2 as two reads form a pair. This would be more of a fragment coverage rather than actual base coverage.

0
Entering edit mode

Yes, that's what I call physical coverage. So I need to write two equations and solve to get the ratio? Is the ratio directly proportional to L1/L2?

0
Entering edit mode

instead of read length use fragment length, though it is not clear from your description what you mean by insert size? how do you know that? If it is by alignment then note that it could be wrong, hence extrapolating from that to physical coverage may also be incorrect