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?
Does your organism have 4+ more genes on the plus strand? That'd be the likely explanation.
Or maybe, on average, the genes on plus strand are 4x more expressed?
Sure, but what could a biological explanation for that be?
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.