This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How can I know forwards reads are sense or anti-sense reads in paired end sequencing ?

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

3 answers

  • 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.

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.

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.

Log in to answer this question.