Tools for strand direction detection RNA-Seq
Hello everyone.
I'm having trouble trying to figure out how to determine the directionality of paired-end reads. I want to find whether the reads are unstranded, first stranded or second stranded. The tool suggestions I've got so far are : Salmon, RSeqC, and how_are_we_stranded_here. Does anyone have anyone experience using these? Which one do you recommend.
• 1,503 views
•
link
updated
by
iraun
•
written
by
pubsurfted •
•
written
by
pubsurfted •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Analize RNA-seq data using Hisat2
written by Maurice •
-
Transcript Directionality with Unstranded Single End Reads ?
written by Will •
-
HTSeq count stranded setting
written by ExtentHonest56 •
-
PolyA Tail Detection in RNA-seq data
written by pubsurfted •
-
Good QC Package (RSeQC or alternative) for Salmon mapping/quantification results?
written by skjw1029 •
-
How can I filter RNA-seq reads to remove reads where more than 50% of the read is low quality (Q<20…
written by mar77 •
-
What will be the impact of strandedness on differentially expressed genes?
written by bioinforesearchquestions
-
Remove Duplicates from Illumina AmpliSeq data
written by preeti.singh •
-
Meta analysis of RNA-seq data
written by bharata1803
-
Stranded Paired-End Data With Rsem
written by Ge •
Easiest is to run infer_experiment.py after mapping your reads against the reference genome.
Not everyone maps rna-seq reads to reference genome to get bam files (I rarely do this because you can't plug bam genome alignments into probabilistic rna-seq models which are required for accurate rna-seq quantification).
The question here is not whether everyone maps the reads against the genome or not. But if you want to use infer_experiment to infer the strand, this is a requirement.
The existing tools you listed are fine; any one of them works.