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

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.1k views
ADD COMMENT
1
Entering edit mode
8.1 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.
ADD COMMENT
0
Entering edit mode
8.1 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.

ADD COMMENT
0
Entering edit mode
8.1 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.

ADD COMMENT

Login before adding your answer.

Traffic: 2649 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6