RSeQC infer_experiment strange output
Entering edit mode
21 months ago


I tried the tool RSeQC on my RNAseq bacterial data. I know that my data is reverse stranded but I get the following output :

This is PairEnd Data
Fraction of reads failed to determine: 0.0693
Fraction of reads explained by "1++,1--,2+-,2-+": 0.3577
Fraction of reads explained by "1+-,1-+,2++,2--": 0.5730

With the doc of RSeQC we would have expected that we obtain like 0.99 for the "1+-,1-+,2++,2--" because it is reverse stranded, or 0.5/0.5 if it wasn't stranded. But we obtain a strange middle ground.

Do you have any idea of what is happening with the data ?

Thank you

RNA-Seq rseqc stranded • 550 views
Entering edit mode
21 months ago

map your reads to the genome of interest and visualize what the alignments look like.

What clarifies a lot of ambiguities is if you align only the data from pair 2 as single-end data.

All the reads from pair 2 should align in the sense (or antisense, depending on the library prep) direction.

Now you can visualize which do not and why.


Login before adding your answer.

Traffic: 1513 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6