SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
DGE using stranded and unstranded Illumina libraries Sentinel156 RNA Sequencing 1 03-22-2016 08:43 AM
CuffDiff Outputs richhogg RNA Sequencing 0 04-29-2015 02:57 AM
Cuffdiff outputs no significant AS genes in RNA-seq analysis. gkandoi RNA Sequencing 0 01-26-2015 08:48 PM
Cuffdiff outputs mostly FAIL agc Bioinformatics 0 01-11-2011 04:45 AM

Reply
 
Thread Tools
Old 12-13-2017, 03:02 PM   #1
capricy
Senior Member
 
Location: 63130

Join Date: Apr 2012
Posts: 125
Default unexplained cuffdiff outputs for stranded libraries

Dear Community,

I recently had sequencing data from stranded library (second strand). The count output from STAR indicated that most of the reads mapped to the second strand. I then ran cuffdiff using the bam files from STAR with -library-type parameter specified.

I happened to find that the either -library-type fr-firststrand or -library-type fr-secondstrand give the fpkms. How should I explain this?

For example:

0610009B22Rik, has 131 reads mapped to the second strand, and only 1 read mapped to the first strand (from STAR count file)

however, it has fpkm of 10.05 when I used "-library-type fr-firststrand", but fpkm of 2.32 when I used "-library-type fr-secondstrand". My expectation was that the latter should be much higher than the former...

Could anyone provide some inputs?

Thanks a lot!

C.

Last edited by capricy; 12-13-2017 at 05:41 PM.
capricy is offline   Reply With Quote
Reply

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:48 PM.


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