why forward reads were much more than reverse reads In strand specific RNAseq data?
0
0
Entering edit mode
5.7 years ago
yuaner • 0

I aligned strand specific RNAseq data using tophat2 and set the parameter --fr-firststrand, After grepped forward strand by 'XS:A:+' and reverse strand by 'XS:A:-', I obtained about x4 forward reads more than reverse reads. I thought it was so strange, why forward reads were much more than reverse reads?

RNA-Seq • 1.7k views
ADD COMMENT
0
Entering edit mode

Does your organism have 4+ more genes on the plus strand? That'd be the likely explanation.

ADD REPLY
0
Entering edit mode

Or maybe, on average, the genes on plus strand are 4x more expressed?

ADD REPLY
0
Entering edit mode

Sure, but what could a biological explanation for that be?

ADD REPLY
0
Entering edit mode

Thanks for your reply. I checked the number of genes on plus and minus strands, but the ratio was about 1:1. And the organism of RNAseq data was rice.

ADD REPLY

Login before adding your answer.

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