SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Bias adjustment in cufflinks: stranded and unstranded protocol question... adeonari Bioinformatics 0 03-21-2013 03:53 AM
stranded or non-stranded protocol? enelkinsan RNA Sequencing 7 12-18-2012 12:42 PM
Unstranded Library RNAseq > Cufflinks --GTF carmeyeii Bioinformatics 0 09-07-2012 08:23 AM
Small RNA library annotation tool? vebaev Bioinformatics 3 08-04-2011 11:28 PM
Library prep from single-stranded DNA? genlyai Sample Prep / Library Generation 3 11-12-2010 12:26 PM

Reply
 
Thread Tools
Old 04-21-2013, 10:33 PM   #1
MichalO
Member
 
Location: CH

Join Date: Jan 2011
Posts: 10
Default tool for checking the stranded vs unstranded library

Assuming that we have by chance a FASTQ from a library that we do not know if was done with stranded or unstranded protocol - how to check this fact?

The library was most likely done as paired end with Tru Seq - whole transcript RNAseq.

A tool for counting antisense matches probably would do - what would you suggest? Meci!
MichalO is offline   Reply With Quote
Old 04-22-2013, 01:47 AM   #2
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994
Default

What about simply _looking_ at your data instead of pulling through some script? You shouldn't go ahead with an analysis without having made at least two or three spot check with a genome browser (e.g., IGV), and that tells you not only about strandedness.
Simon Anders is offline   Reply With Quote
Old 04-22-2013, 04:35 AM   #3
mknut
Member
 
Location: UK

Join Date: Jul 2012
Posts: 23
Default

Here is the answer (with pretty pictures). I agree with simon - checking manually a few reads will not only tell you about the strandedness, but will also give you an idea about the insert length.
mknut is offline   Reply With Quote
Old 04-29-2013, 04:47 AM   #4
MichalO
Member
 
Location: CH

Join Date: Jan 2011
Posts: 10
Default

Thanks, MK!

Yup, Simon - looking in-depth is the first choice for me too
But we see some strange, antisense-looking thingies in IGV in some of the places- and some not very relaxed people call for a "metrics"
MichalO is offline   Reply With Quote
Old 04-29-2013, 06:25 AM   #5
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994
Default

And these "not very relaxed" persons think that a good way of dealing with such strange looking worrying features is to condense them into a single number in the hope that that makes them go ahead? That actually rather sounds too relaxed to me. ;-)

If you must, you could count everything with htseq-count, once with '-s yes' and once with '-s reverse' and make a scatter plot of counts per gene on the sense versus on the antisense strand.

BTW, don't forget that there are lots of antisense transcripts in many species.
Simon Anders is offline   Reply With Quote
Old 04-29-2013, 08:49 AM   #6
kcchan
Senior Member
 
Location: USA

Join Date: Jul 2012
Posts: 182
Default

If you want a very quick and dirty check for strandedness, you can look for A/T and G/C skew in your reads. Stranded libraries will have a visible skew while unstranded libraries should have an even A/T and G/C skew. Of course, this isn't fool proof, so definitely go more in depth as others have suggested if you want to make sure.
kcchan is offline   Reply With Quote
Reply

Tags
library, quality control, stranded, unstranded

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:57 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO