This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Multiple columns in RNAseq count matrix?

Hi, in GSE131592, each file corresponds to the counting matrix of one sample. In each file there are three columns. Does anyone know how to get the count of this sample?

The picture shows GSM3790428 as an example. No reply from the dataset author yet.

example

rnaseq
• 3,097 views • link
updated by swbarnes2 swbarnes2 • written by Jeff-Gui • Jeff-Gui

1 answer

That is htseq count data. The three columns represent the number of reads assigned to that gene assuming the library is 1) unstranded 2) stranded in the forward direction, 3) stranded in the reverse direction.

• 0 views • link
written by swbarnes2 swbarnes2

Thanks. Without checking raw reads to see how exactly the library is generated, is it reasonable to simply take the second column because of the much fewer ambiguous count than the first?

• 0 views • link
written by Jeff-Gui • Jeff-Gui

I'd look over more of the columns, and if the general trend is that columns 1 and 2 are nearly the same, that means the prep is stranded forward, and you should use column 2.

• 0 views • link
written by swbarnes2 swbarnes2

Log in to answer this question.

More posts like this