Tophat 1.4.1 paired end alignment error -11 healp

mediator · 02-22-2012, 08:07 AM

Hi All,
I was running Tophat1.4.1 to do a pair-end alignment but got err=-11. Here is the command I used:

module load bowtie/0.12.7
module load tophat/1.4.1
tophat -m 2 -p 16 -r 210 -o <output dir> <hg19_index> /data/PE_read_1.fastq /data/PE_read_2.fastq

Output:

Quote:

[Tue Feb 21 10:51:53 2012] Beginning TopHat run (v1.4.1)
-----------------------------------------------
[Tue Feb 21 10:51:53 2012] Preparing output location <output dir>
[Tue Feb 21 10:51:53 2012] Checking for Bowtie index files
[Tue Feb 21 10:51:53 2012] Checking for reference FASTA file
[Tue Feb 21 10:51:53 2012] Checking for Bowtie
Bowtie version: 0.12.7.0
[Tue Feb 21 10:51:53 2012] Checking for Samtools
Samtools Version: 0.1.16
[Tue Feb 21 10:51:53 2012] Generating SAM header for <hg19>
format: fastq
quality scale: phred33 (default)
[Tue Feb 21 10:52:19 2012] Preparing reads
left reads: min. length=50, count=121496435
right reads: min. length=50, count=121483750
[Tue Feb 21 12:08:01 2012] Mapping left_kept_reads against genome with Bowtie
[Tue Feb 21 14:28:21 2012] Processing bowtie hits
[Tue Feb 21 16:06:02 2012] Mapping left_kept_reads_seg1 against genome with Bowtie (1/2)
[Tue Feb 21 16:52:36 2012] Mapping left_kept_reads_seg2 against genome with Bowtie (2/2)
[Tue Feb 21 17:54:56 2012] Mapping right_kept_reads against genome with Bowtie
[Tue Feb 21 20:20:10 2012] Processing bowtie hits
[Tue Feb 21 22:00:24 2012] Mapping right_kept_reads_seg1 against genome with Bowtie (1/2)
[Tue Feb 21 22:58:31 2012] Mapping right_kept_reads_seg2 against genome with Bowtie (2/2)
[Tue Feb 21 23:58:39 2012] Searching for junctions via segment mapping
[FAILED]
Error: segment-based junction search failed with err =-11
Loading left segment hits...

Any suggestions? Thanks!

pbluescript · 02-23-2012, 08:02 AM

I would try dropping the -m 2 option. The issue might be that by allowing 2 mismatches in an anchor region only 8 bp long, you are getting a lot of potential matches that are making junction mapping difficult.

mediator · 02-23-2012, 08:09 AM

@pbluescript
Will try that, thanks!

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
tophat paired end read alignment	idyll_ty	Bioinformatics	0	02-18-2012 10:58 PM
Tophat paired end alignment help: comma or blank to separate two pairs?	mediator	Bioinformatics	3	01-31-2012 09:44 AM
Tophat paired end alignment.	adeslat	SOLiD	5	08-07-2011 04:05 AM
Tophat - Error using Paired End	jling	Bioinformatics	0	02-10-2010 11:26 AM
how paired end alignment works?	totalnew	Bioinformatics	8	04-27-2009 01:46 PM

02-23-2012, 08:02 AM	#2
pbluescript Senior Member Location: Boston Join Date: Nov 2009 Posts: 224	I would try dropping the -m 2 option. The issue might be that by allowing 2 mismatches in an anchor region only 8 bp long, you are getting a lot of potential matches that are making junction mapping difficult.

02-23-2012, 08:09 AM	#3
mediator Member Location: New England Join Date: Nov 2010 Posts: 27	@pbluescript Will try that, thanks!