SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Running tophat on a cluster NM_010117 Bioinformatics 14 10-31-2016 08:16 AM
tophat Error running running 'prep_reads' victoryhe Bioinformatics 2 10-17-2011 04:53 AM
Running Tophat for a few samples dnusol Bioinformatics 0 01-28-2011 12:59 AM
a tophat running problem clariet Bioinformatics 1 04-12-2010 05:05 AM
32bp? yamayaya Illumina/Solexa 5 02-03-2009 04:39 AM

Reply
 
Thread Tools
Old 11-19-2009, 07:12 AM   #1
statsteam
Member
 
Location: Californica

Join Date: Sep 2009
Posts: 18
Default Running TopHat with 32bp or 36bp reads

Hi all,

I am trying to run TopHat with 32bp or 36bp reads. I am setting the length of segments to the half of the original length (16bp for 32bp reads and 18bp for 36bp reads). Besides it takes too long to map such short reads to the genome, I am not so sure about the reliability of the results because there are going to be simply too many mappable hits on the genome.

Any thoughts on this?

Thank you,
Statsteam
statsteam is offline   Reply With Quote
Old 11-20-2009, 06:55 PM   #2
HTS
Member
 
Location: Toronto

Join Date: Nov 2009
Posts: 24
Default

Further breaking down reads as short as 32 or 36bp doesn't seem to be a sensible thing to me. I wouldn't use segments shorter than 25bp myself. Doesn't the default setting more or less work for you in your case? I used the -a 5 option with 32bp reads and things seem to work fine for me. To further increase sensitivity, you could supply known junctions via a GFF3 file or tune down the -F parameter.

BTW, I saw in your other post that you have worked with Burge's cerebellum data. Where did you get it? By following the accession numbers provided in their Nature paper, I get all but the cerebellum data of six individuals. Thanks!

-- Leo
HTS is offline   Reply With Quote
Reply

Tags
mapping, tophat

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 07:00 AM.


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