Hi,
I am working on a RNAseq data set of about 20 human tissue samples. However, some of the samples have a very low alignment rate (~30-40%) while the rest has a a decent alignment rate around 80%. I aligned the reads using STAR with default settings. The log file tells me that for those samples with low alignment rate many reads did not map due to multimapping. My first guess is that there could have been an issue with the rRNA depletion for those samples. How can I show that the rRNA depletion failed?
Take unmapped reads and align to rRNA sequences? Has anybody had similar issues before?
Any feedback is greatly appreciated!
•
written
by
lkmklsmn
1 answer
Yes, you can align against the 45S sequence, noting that this will (in my experience) somewhat underestimate the actual amount of rRNA in the sample.
As an aside, I should note that rRNA depletion only seems to work well if you do it on fresh samples. If you freeze them beforehand then the % of rRNA remaining is more likely to be high (this is purely observational).
Low RINs tend to correlate with degradation, so you might give salmon a try with whatever the position bias correction option is. Perhaps that'd help.
Log in to answer this question.
-
[RNASeq] Excluding samples that cluster differently in PCA plot
written by Texx •
-
Salmon low mapping rates
written by L.Chang •
-
High amount of intronic/intergenic reads in SMARTer stranded total bulk RNAseq
written by Mat •
-
anyway to save contaminated samples?
written by demoraesdiogo2017 •
-
RNA-seq pipeline for degraded RNA
written by oludhe •
-
What % of rRNA contamination is acceptable?
written by John
-
STAR output - high multi-mapping, low alignment (human)
written by anikb •
-
Is it normal to have low expression levels of house keeping genes in plant RNA-seq results
written by shuchen •
-
Ribosomal RNA in RNA-seq data
written by lkmklsmn
-
Tophat vs Tophat-fusion concordant alignment rate
written by IP
