How to get actual alignment % from STAR?
1
0
Entering edit mode
18 months ago
John ▴ 220

Hi

Is this formula right to calculate actual alignment percentage of my RNA seq data? I want to find out the reads mapped only to annotated genes, my sample is from Ribosome bound RNA, so I would not have introns, I want to exclude rRNA alignments from the percentage alignment.

   Alignment % = Total reads - (N_unmapped - N_multimapping - N_noFeature -N_ambiguous)

I would get these four value from ReadsPerGene.out.tab file. Thanks

RNA-Seq R • 429 views
ADD COMMENT
0
Entering edit mode
18 months ago
padwalmk ▴ 100

Hi, You can use MultiQC which gives everything in a very good format.

ADD COMMENT

Login before adding your answer.

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