If I have reads of length x, inserrt size of y and genome size z, how do I convert from sequence to physical coverage?
•
written
by
Dan
1 answer
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.
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
Log in to answer this question.
-
bedtools: coordinates of read covered
written by Smriti •
-
Is my estimated genome coverage correct?
written by apcreyes29 •
-
How to generate sequence length distribution from Fasta file
written by 2822462298 •
-
coverage of the reads over coding sequence only for one gene
written by alireza346 •
-
How to calculate the joint entropy H(X,Y,Z) where X, Y, Z are position of MSA
written by sivanc7 •
-
Coverage for sample with many species
written by biobio •
-
Combining Dna Sequences Files Into One
written by nepgorkhey
-
Ggbio Issues: How To Align Tracks & Where Are Figure Legends?
written by sndrtj
-
How Can I Draw A Rod In Pymol?
written by yotama9 •
-
Finding Genome Coverage Using Random Reads
written by Smandape
maybe you miss n the number of reads?
n is implied by the coverage