Question: RSeQC infer_experiment strange output
0
gravatar for nonolescargot
13 days ago by
nonolescargot0 wrote:

Hello

I tried the tool RSeQC infer_experiment.py 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

stranded rna-seq rseqc • 75 views
ADD COMMENTlink modified 13 days ago by Istvan Albert ♦♦ 82k • written 13 days ago by nonolescargot0
0
gravatar for Istvan Albert
13 days ago by
Istvan Albert ♦♦ 82k
University Park, USA
Istvan Albert ♦♦ 82k wrote:

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.

ADD COMMENTlink modified 13 days ago • written 13 days ago by Istvan Albert ♦♦ 82k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1542 users visited in the last hour