SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Alignment PE reads with different length by TopHat Trudy Bioinformatics 4 10-17-2013 11:47 AM
segment length in tophat ickou RNA Sequencing 4 05-21-2013 12:30 PM
Limit of indel length with bwa nvteja Bioinformatics 1 12-09-2010 08:46 AM
Tophat --segment-length Pepe Bioinformatics 0 05-31-2010 10:33 PM
seed length in TopHat? joseph Bioinformatics 5 09-23-2009 09:23 PM

Reply
 
Thread Tools
Old 10-21-2009, 01:10 AM   #1
maria.b
Member
 
Location: Paris

Join Date: Sep 2009
Posts: 14
Default TopHat limit length

Dear all,

I'm trying to use Tophat with 8 million single-end reads of 75 bp. I set all options that concern minimum intron length to 30 (because the minimum intron length in my genome reference is 39), all other option are set to the default value:
-i 30
--min-coverage-intron 30
--min-segment-intron 30
--min-closure-intron 30

In this conditions it works until the step of "Joining segment hits", but it never passed this step. I'm still waiting since yesterday morning!!!!
Is it normal? How many times does it take?

If I test this conditions on 8 millions single-end read of 35bp it completly works in ~25 minutes.

thanks for your help

Maria
maria.b is offline   Reply With Quote
Old 10-21-2009, 08:32 AM   #2
Cole Trapnell
Senior Member
 
Location: Boston, MA

Join Date: Nov 2008
Posts: 212
Default

How are you running the program (i.e. command line options)? Which version are you using?

It's not normal: I routinely run more than 150 million 75bp reads, and that step takes no more than an hour (when our filesystem is running slowly).
Cole Trapnell is offline   Reply With Quote
Old 10-21-2009, 11:50 PM   #3
maria.b
Member
 
Location: Paris

Join Date: Sep 2009
Posts: 14
Default

Hi,

I use the last version of bowtie (0.11.3) and the last version of tophat (1.0.11) with the command line:

./tophat-1.0.11/bin/tophat -i 30 --min-coverage-intron 30 --min-segment-intron 30 --min-closure-intron 30 -o res-70/ index/trichoderma-ref 70a.fq

The same commande line works well on the 35bp read.
maria.b 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:01 PM.


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