Have you ever seen perfect (100%) mapping rate (RNASeQC) on human RNASeq (hg38) using STAR?
1
0
Entering edit mode
5.6 years ago
mark.ebbert ▴ 20

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 • 2.1k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY
1
Entering edit mode
5.6 years ago
h.mon 34k

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.

ADD COMMENT
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

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