How can I know forwards reads are sense or anti-sense reads in paired end sequencing ?
3
0
Entering edit mode
8.4 years ago
jack ▴ 960

I have Illumina paired-end sequencing read. I want to do quantification and I need to specify to the software, which read(forward, backward) is sense and which read is anti-sense.

How can I know that, forward reads are sense and the backwards are antisense or vice versa?

software alignment error Assembly RNA-Seq • 4.2k views
1
Entering edit mode
8.4 years ago
• It's dependent on how the libraries were made, so just ask whomever made them (having said that, read#2 is commonly determines the strand).
• Just look at a couple genes in IGV. If you color reads by the orientation of read#1 then it should be quickly apparent which (if either) read is determining the strand.
0
Entering edit mode
8.4 years ago
Prakki Rama ★ 2.6k

Quick BLAST must do this. Take subset of reads. Blast it again few reference genes and see how your forward reads and reverse reads are aligning by looking at the subject start and subject end. Check if most of your reads are aligning (-) orientation or in (+) orientation.

0
Entering edit mode
8.4 years ago
Pascal ▴ 250

You could also use infer_experiment.py from the RSeQC package. It will tell you to which strand your read map (and if it is the second or first read) or if your library is not strand specific.