This is a test version of Biostars. For the public version, visit https://www.biostars.org.
strandness of paired-end data

I have used infer_experiment.py from RSeQC and got this output:

This is PairEnd Data
Fraction of reads failed to determine: 0.0659
Fraction of reads explained by "1++,1--,2+-,2-+": 0.0055
Fraction of reads explained by "1+-,1-+,2++,2--": 0.9285

Is it shows that my RNA-Seq data is strand specific and strandness of read1 is opposite with that of gene model, while strandness of read2 is consistent to the strand of reference gene model?

Is that means that its reverse stranded and I have to use s -2 option or s -1 option in featureCounts?

rna-seq featurecounts
• 3,916 views • link
updated by DareDevil DareDevil • written by waqaskhokhar999 waqaskhokhar999

This is not a "Forum" this is a Question type post. I have corrected it for you this time.

• 1 views • link
written by Joe Joe

1 answer

For dummies, if the output is:

Fraction of reads explained by "1++,1--,2+-,2-+": 0.0055
Fraction of reads explained by "1+-,1-+,2++,2--": 0.9285

You can use -s 2 option for featurecounts

or if the output is:

Fraction of reads explained by "1++,1--,2+-,2-+": 0.9285
Fraction of reads explained by "1+-,1-+,2++,2--": 0.0055

You can use -s 1 option for featurecounts

or if the output is:

Fraction of reads explained by "1++,1--,2+-,2-+": 0.4555
Fraction of reads explained by "1+-,1-+,2++,2--": 0.4585

You can use -s 0 option for featurecounts

• 0 views • link
written by DareDevil DareDevil

Log in to answer this question.

More posts like this