Why do ribo-seq sequencing data have only a small proportion of uniquely mapped reads?
2
0
Entering edit mode
6.7 years ago
Cacau ▴ 520

I downloaded several published ribo-seq data and mapped them with STAR. However, only 15% ~ 50% of them are uniquely mapped reads. The reads are short, the length being about 35bp long after trimming. Is this the reason? Or might be because many rRNA reads are included?

ribo-seq RNA-Seq • 3.4k views
ADD COMMENT
2
Entering edit mode
6.7 years ago
mark.ziemann ★ 1.9k

You're correct in that it could just be short length or rRNA. It is relatively easy to check the following:

  1. Screen rRNA content http://genomespot.blogspot.com.au/2015/08/screen-for-rrna-contamination-in-rna.html
  2. Simulate cDNA reads of the different lengths and determine the proportion that map uniquely http://genomespot.blogspot.com.au/2014/10/rna-seq-aligner-accuracy-tested-with.html
ADD COMMENT
0
Entering edit mode

That's very helpful! Thank you.

ADD REPLY
2
Entering edit mode
6.7 years ago
mark.ziemann ★ 1.9k

You're correct in that it could just be short length or rRNA. It is relatively easy to check the following:

  1. Screen rRNA content http://genomespot.blogspot.com.au/2015/08/screen-for-rrna-contamination-in-rna.html
  2. Simulate cDNA reads of the different lengths and determine the proportion that map uniquely http://genomespot.blogspot.com.au/2014/10/rna-seq-aligner-accuracy-tested-with.html
ADD COMMENT

Login before adding your answer.

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