SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

Reply
 
Thread Tools
Old 02-22-2012, 08:07 AM   #1
mediator
Member
 
Location: New England

Join Date: Nov 2010
Posts: 27
Default Tophat 1.4.1 paired end alignment error -11 healp

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!
mediator is offline   Reply With Quote
Old 02-23-2012, 08:02 AM   #2
pbluescript
Senior Member
 
Location: Boston

Join Date: Nov 2009
Posts: 224
Default

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.
pbluescript is offline   Reply With Quote
Old 02-23-2012, 08:09 AM   #3
mediator
Member
 
Location: New England

Join Date: Nov 2010
Posts: 27
Default

@pbluescript
Will try that, thanks!
mediator 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 10:59 PM.


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