SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat: segment-based junction search failed with err =-11 jdanderson Bioinformatics 15 10-14-2017 06:47 AM
Tophat2.0.8b segment-junc err=1 mattia Bioinformatics 8 01-08-2014 02:38 PM
tophat2 segment_juncs error: Error: segment-based junction search failed with err =-6 hulan0@gmail.com Bioinformatics 1 04-16-2012 06:37 AM
Tophat: segment join failed with err = -11 ytmnd85 Bioinformatics 4 03-15-2012 03:14 AM
TopHat ERROR: Segment join failed with err = 1 rebrendi Bioinformatics 10 12-17-2011 11:58 AM

Reply
 
Thread Tools
Old 05-07-2013, 01:25 AM   #1
mattia
Member
 
Location: Milano

Join Date: Aug 2011
Posts: 30
Question Tophat2.0.8b segment-junc err=1

Hi,
Last week I upgraded last version of Tophat; so now I'm using Tophat 2.0.8b with bowtie 1.0.0 (I need to align color space reads).
Before I used Tophat 1.4.1 and bowtie 0.12.8, and I never have any problems.

So I launch this command (the same of previous version except for --bowtie1) :
Quote:
tophat -C --bowtie1 --no-coverage-search --library-type fr-secondstrand -G $annotation_file -g 10 -p 12 -r 90 -o $output_path $reference_file $csfastaF3 $csfastaF5
with reference_file="/.../ENSEMBL_HG19/Homo_sapiens.GRCh37.71.dna.23chr_c"

So, now this message appers:
Quote:
[2013-05-06 16:04:37] Beginning TopHat run (v2.0.8b)
-----------------------------------------------
[2013-05-06 16:04:37] Checking for Bowtie
Bowtie version: 1.0.0.0
[2013-05-06 16:04:37] Checking for Samtools
Samtools version: 0.1.18.0
[2013-05-06 16:04:37] Checking for Bowtie index files
[2013-05-06 16:04:37] Checking for reference FASTA file
Warning: Could not find FASTA file /home/references/new_Human_annotations/ENSEMBL_HG19/Homo_sapiens.GRCh37.71.dna.23chr_c.fa
[2013-05-06 16:04:37] Reconstituting reference FASTA file from Bowtie index
Executing: ....

....
....


[2013-05-06 16:06:57] Reading known junctions from GTF file
[2013-05-06 16:07:11] Preparing reads

WARNING: read pairing issues detected (check prep_reads.log) !

left reads: min. length=75, max. length=75, 49812297 kept reads (187703 discarded)
right reads: min. length=35, max. length=35, 49679413 kept reads (320587 discarded)
[2013-05-06 16:20:43] Creating transcriptome data files..
[2013-05-06 16:21:59] Building Bowtie index from Homo_sapiens.GRCh37.71.23chr.fa
[2013-05-06 16:52:43] Mapping left_kept_reads to transcriptome Homo_sapiens.GRCh37.71.23chr with Bowtie
[2013-05-06 17:07:08] Mapping right_kept_reads to transcriptome Homo_sapiens.GRCh37.71.23chr with Bowtie
[2013-05-06 17:22:15] Resuming TopHat pipeline with unmapped reads
[2013-05-06 17:22:16] Mapping left_kept_reads.m2g_um to genome Homo_sapiens.GRCh37.71.dna.23chr_c with Bowtie
[2013-05-06 18:19:33] Mapping left_kept_reads.m2g_um_seg1 to genome Homo_sapiens.GRCh37.71.dna.23chr_c with Bowtie (1/3)
[2013-05-06 18:53:33] Mapping left_kept_reads.m2g_um_seg2 to genome Homo_sapiens.GRCh37.71.dna.23chr_c with Bowtie (2/3)
[2013-05-06 19:28:16] Mapping left_kept_reads.m2g_um_seg3 to genome Homo_sapiens.GRCh37.71.dna.23chr_c with Bowtie (3/3)
[2013-05-06 20:01:48] Mapping right_kept_reads.m2g_um to genome Homo_sapiens.GRCh37.71.dna.23chr_c with Bowtie
[2013-05-06 20:21:52] Mapping right_kept_reads.m2g_um_seg1 to genome Homo_sapiens.GRCh37.71.dna.23chr_c with Bowtie (1/1)
[2013-05-06 20:26:56] Searching for junctions via segment mapping
[FAILED]
Error: segment-based junction search failed with err =1
Error: could not get read# 11562180 from stream!
Any suggestions (maybe also from autors)?

Thank a lot.
mattia is offline   Reply With Quote
Old 05-08-2013, 12:12 AM   #2
Jeremy
Senior Member
 
Location: Pathum Thani, Thailand

Join Date: Nov 2009
Posts: 190
Default

This is a duplicate of http://seqanswers.com/forums/showthread.php?t=29981
Jeremy is offline   Reply With Quote
Reply

Tags
errors, segment_juncs, tophat 2

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 06:49 AM.


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