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

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.4k views
ADD COMMENT
2
Entering edit mode

maybe you miss n the number of reads?

ADD REPLY
0
Entering edit mode

n is implied by the coverage

ADD REPLY
1
Entering edit mode
10.7 years ago

N = number of reads, L=length of each read, G=genome size

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.

ADD COMMENT
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?

ADD REPLY
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

ADD REPLY

Login before adding your answer.

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