Question: Ribosomal RNA in RNA-seq data
gravatar for lkmklsmn
4.2 years ago by
United States
lkmklsmn930 wrote:

Hi biostars,
I am analyzing RNA-seq data (not for the first time) but for a subset of samples I am getting a very low alignment rate (~30%). To find out what is going on I have looked at some of the unaligned reads and they seem to come from ribosomal RNA.
Now my question: Are ribosomal RNA sequences not contained in the standard reference genome sequence (e.g. hg19)? Is that why the reads did not align or should I be able to map rRNA reads to hg19?


ribosomal rna-seq mapping rrna • 3.0k views
ADD COMMENTlink written 4.2 years ago by lkmklsmn930

Reads mapping to rRNA are generally multi mappers and they might have been excluded due to the parameter you selected for mapping. In other words, did you allow multi mapping reads or not ? You could also directly map the reads only to rRNA and see what fraction of reads map to rRNA and how much reads are left ?

ADD REPLYlink written 4.2 years ago by Chirag Nepal2.3k

Thanks. Yes the percentage of multi-mapping reads is very high in those affected samples. Given that only a subset of the samples are affected by this (Alignment rate<50%) my conclusion is that there must have been some error in the rRNA depletion step.

ADD REPLYlink written 4.2 years ago by lkmklsmn930
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1949 users visited in the last hour