SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat junctions.bed RockChalkJayhawk RNA Sequencing 7 12-12-2013 10:56 AM
Junctions in tophat jay2008 Bioinformatics 0 01-26-2012 06:23 PM
Missing Junctions in Tophat! ( after providing a known junctions & gene models files) avi Bioinformatics 2 08-02-2011 04:21 PM
tophat - no junctions crh Bioinformatics 2 07-25-2011 11:33 AM
tophat: junctions and reads Nicolas902 Bioinformatics 0 03-09-2010 02:01 AM

Reply
 
Thread Tools
Old 01-05-2011, 09:51 AM   #1
agent99
Member
 
Location: San Francisco

Join Date: Jul 2010
Posts: 10
Default junctions in tophat

When using TopHat v1.2.0, I am getting 5-fold more junctions identified when I give it an annotation file than when I do not give an annotation file.

12471 junctions - no annotation provided
58610 junctions - annotation provided, no novel junctions allowed
60048 junctions - annotation provided, novel junctions allowed

These are for 36 bp paired-end seqs from a mouse cell line.

Here's my (abbreviated) call with annotation:

tophat --GTF mm9.NCBIM37.59.fix.gtf --mate-inner-dist 104 --mate-std-dev 17 --solexa1.3-quals --segment-length 18 --keep-tmp --output-dir out mm9 left_sequence.txt right_sequence.txt &

Looks like TopHat has trouble finding junctions when it doesn't have annotation information. Has anyone else seen this? Thanks for input.
agent99 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 09:03 PM.


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