This is a test version of Biostars. For the public version, visit https://www.biostars.org.
What Does "Length Per Kilobase" Mean ?

Let's say if I have to find "length per KB" of a read.

How should I go about it ?

I can find the "length" by subtracting the END coordinate with the read with the START coordinate.

Let's say Read 1 spans from 300 to 330. So, it's length is 30.

Does "Length per KB" mean I divide 30 by 1000 i.e 0.030 ?

length
• 4,644 views • link
updated by Biostar Biostar • written by Tarun Mohite • Tarun Mohite

2 answers

The "length per kilobase" as you pointed is the sequence size in kilobases units. But it's applied to gene sizes which spans for 1-2kb in average in the human genome, no reads that are really short. I think you are trying to understand RPKM see http://seqanswers.com/forums/showthread.php?t=586

• 0 views • link
written by JC JC

JC, thanks a lot ! RPKM is exactly what I am trying to understand. :)

• 0 views • link
written by Tarun Mohite • Tarun Mohite

Hello JC,

According to my understanding of the paper, the simplest way I can write a wordy formula for RPKM is :

RPKM = (sum of all reads mapping to all exons of the gene) / ((sum of lengths of all exons of the gene / 1000) * (sum of all reads mapping to all genes in the sample / 1000000))

Is that correct ?

However, if I have to find the RPKM for each exon of each gene, would I modify the formula to :

RPKM = (total #reads mapping to the exon) / ((length of exon / 1000) * (sum of all reads mapping to all genes in the sample / 1000000))

Please correct if I am going wrong.

Thanks.

• 2 views • link
written by Tarun Mohite • Tarun Mohite

Log in to answer this question.

More posts like this