Counts from STAR with stranded protocol
2
1
Entering edit mode
2.1 years ago
Jamie ▴ 10

Hi,

I have some RNA-seq fastq files which I aligned to GRCh38 with STAR. The RNA-seq protocol was stranded. I used --quantMode GeneCounts option in STAR. An example of one of the count files can be seen below. I'm looking to run DESeq2 and wanted to check if the right method is to use the counts from the third column?

Single sample gene count from STAR

thanks

RNA-seq DEG STAR DESeq2 • 1.0k views
ADD COMMENT
0
Entering edit mode
2.1 years ago

Use the far right column, not the nearly empty column you've circled. It looks like this library prep returns reads in the reverse direction.

ADD COMMENT
0
Entering edit mode

Ah okay, so the column with lower count values is defined as the reverse strand?

ADD REPLY
0
Entering edit mode

No. Read the documentation.

ADD REPLY
0
Entering edit mode
2.1 years ago
jkim ▴ 170

The column you referred to is counts for the forward stranded. Looks like yours is reverse stranded sequencing data. This is a screen-capture of STAR manual.

enter image description here

ADD COMMENT

Login before adding your answer.

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