Hi I want to ask a question about calculate the RPKM.
I see the previous people calculate the RPKM for a gene, by using all the reads included in intron and exons belong to a gene.
If I just want to get the RPKM for different isoforms. Can I just pick up the reads mapped to exons.
For example, Gene A has 4 exons:
Chr2 MSU_osa1r6 gene 4542759 4544980 . + . ID=13102.t00754;Name=unknown gene Chr2 MSU_osa1r6 CDS 4543031 4543177 . + 0 Parent=13102.m00974 Chr2 MSU_osa1r6 CDS 4543287 4543709 . + 0 Parent=13102.m00974 Chr2 MSU_osa1r6 CDS 4543836 4543952 . + 0 Parent=13102.m00974 Chr2 MSU_osa1r6 CDS 4544064 4544423 . + 0 Parent=13102.m00974
There are 4,011 reads that map to this gene (between positions 4542759 and 4544980).
There are 4 exons for this particular gene which contain a total of 1,043 base pairs.
So the RPKM for this particular gene is ((4,011 reads/1.043kb of exon)/31.8mill mapped reads) = 120.9RPKM
If I just calculate the RPKM for different isoforms, I just want to use the reads in exons, so it is about 3000 reads.
So my RPKM is ((3000/1043)/31.8mill mapped reads)=90RPKM.
So if I just want to compare the expression level of this isofoms across samples, is it ok to do that?
Thanks in advance.