This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Rpkm In Masked Sequences

Hello.

I'm confused about calculating RPKM when the target sequence is masked.

In this post, they discuss about the use of "N", in

rpkm=10^9*C/NL, where C is the reads number of the transcript, L is the length of the transcript and N is the total reads count of the sample.

http://seqanswers.com/forums/showthread.php?t=7630

And what about L? Does it change, when the sequence is masked?

Thank you.

rpkm
• 2,514 views • link
written by Leandro Lima Leandro Lima

1 answer

If you mapped the reads to the masked sequence, then I would use the masked length to calculate the RPKM.

The rationale is that you are effectively removing reads that would have mapped to that masked region, so you should also remove the length of the masked region when normalizing by length.

• 128 views • link
written by Damian Kao Damian Kao

Log in to answer this question.

More posts like this