SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Question:Most convenient way to get annotation information for a chromosome location? serenaliao Epigenetics 2 02-14-2013 12:36 PM
Information about RNA-Seq Coverage oiluig87 Bioinformatics 0 11-08-2012 01:23 AM
Question: Information about the Caltech RNA-Seq data on ENCODE and its corresponding tmy1018 Bioinformatics 1 08-15-2012 01:04 PM
samtools question - redefine mate information after applying filters caswater Bioinformatics 0 03-07-2012 08:37 PM

Reply
 
Thread Tools
Old 07-28-2015, 09:46 AM   #1
ml439
Junior Member
 
Location: Cambridge

Join Date: May 2015
Posts: 9
Default question about RNA-Seq and strandness information

Hi,
I'm fairly new to RNA-Seq data analysis and I'm a bit confused about something and I'd like to get your views on it. I have run an RNA-Seq analysis for a set of samples that i have been told they were prepared as second strand. So I've run tophat2 (didn't include strandness there, left if as default), then I counted reads with HTSeq (using -s reverse), and done DESeq and DEXSeq. I was happy with my data, since I could detect our positive control splicing events, until I had to run another data set, and something went wrong.
I've contacted the people that had prepared the libraries and run the sequencing for us and it turns out that the first data set had to be run using fr-firststrand.
Now I'm very confused, why does my analysis show differential splicing in positive control splicing events, if it seems that i have run the wrong commands? is the fact that I run the mapping as default the explanation? but then how about counting with HTSeq? surely that should be counting the wrong strand and throwing away a lot of useful data?

I would appreciate if anyone could help, i'm extremely confused here.

Thank you very much for your time and help,

Miriam
ml439 is offline   Reply With Quote
Reply

Tags
first strand, htseq, second strand, tophat2

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 06:03 AM.


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