SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
TopHat Error: Could not find Bowtie index files /bowtie-0.12.5/indexes/. rebrendi Bioinformatics 11 06-22-2016 10:55 AM
how to get coverage from TopHat output? joseph Bioinformatics 8 04-23-2012 05:44 PM
Maximum possible coverage depth junfeng Bioinformatics 1 01-18-2012 09:54 AM
TopHat 1.3.1 coverage islands Irina Pulyakhina Bioinformatics 0 08-15-2011 03:24 AM
maximum number of the reference input files for bowtie-build joseph Bioinformatics 2 09-04-2010 03:09 PM

Reply
 
Thread Tools
Old 03-17-2010, 11:13 AM   #1
Ender985
Member
 
Location: Spain

Join Date: Mar 2009
Posts: 12
Default Bowtie/TopHat hitting maximum coverage

Hey all,

I've been working with an RNA paired-end dataset of around 60 milion reads, in which a big part of the reads map to the same gene. However, it seems that there is a hard limit of around 65536 (=2^16) mapped reads per base, since this gene displays a weird saw-like pileup (going up steadily until around 65536 and then dropping to around 1000 to keep going up steadily), in contrast to the normal appearance of other less expressed genes that don't reach that coverage level, or even the very same gene in an alignment with a less deep coverage dataset:



(click for zoom) Blue track is the less deep dataset, grey track is this problematic pileup. The maximum value shown here is not 65000 reads, since I normalized the tracks before uploading them, it is the raw coverage.wig file TopHat produces that actually hits that value in the grey track. The effect is best seen in the left part of the 2nd exon.

I'm using the latest Bowtie version, 0.12.3.0, and TopHat version v1.0.12; even though when I first spotted this behavior I was using a previous Bowtie version, 0.11.3.0. This was run on a x86-64 Karmik Ubuntu.

I don't think there is an easy way around this, using comand line options and such, but I'm very interested to know. I can write a script that will deal with this problem, but probably this issue can be solved for the future versions of Bowtie/TopHat. I'm no IT guy but it seems like a simple int overflowing problem to me.

Last edited by Ender985; 03-17-2010 at 11:56 AM.
Ender985 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 01:55 AM.


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