This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Counts from STAR with stranded protocol

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
• 2,087 views • link
updated by swbarnes2 swbarnes2 • written by Jamie • Jamie

2 answers

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.

• 0 views • link
written by swbarnes2 swbarnes2

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

• 0 views • link
written by Jamie • Jamie

No. Read the documentation.

• 0 views • link
written by swbarnes2 swbarnes2

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

• 0 views • link
written by jkim • jkim

Log in to answer this question.

More posts like this