This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Have you ever seen perfect (100%) mapping rate (RNASeQC) on human RNASeq (hg38) using STAR?

Hi,

I've been out of the RNA-Seq world for a while, so I'm trying to get up to speed on the latest expectations, protocols, QC, etc.

We've run four test samples (human blood) for our specific protocol (hemoglobin AND rRNA depletion, obviously no poly-A selection). I used STAR to align to the hg38 genome (annotated using gencode v25). I then ran FastQC, RNASeQC, and RSeQC (via QuaCRS) and am digging through the various QC metrics. I was immediately surprised to find the mapping rate is 1 (100%) for all four samples, according to both RSeqC and RNASeQC. The "Mapped Unique Rate" ranges from 70%-80% (according to RNASeQC).

It seems very strange that every read maps to the human genome, even from blood, but this is the first time I've run a blood RNASeq. Based on various google searches, I've only found posts concerned about low mapping rates.

I really appreciate any thoughts!

rna-seq hg38 star alignment
• 3,387 views • link
updated by h.mon h.mon • written by mark.ebbert • mark.ebbert

Could you tell us which STAR version and your exact command-line?

• 0 views • link
written by h.mon h.mon

1 answer

Search the STAR manual for --outSAMunmapped, and see an explanation here.

edit: check the final.log files, it should state the number of unmapped reads.

• 0 views • link
written by h.mon h.mon

uggh, you're right. I was looking for much more complicated answers than that.

• 0 views • link
written by mark.ebbert • mark.ebbert

Log in to answer this question.

More posts like this