Question: TOPHAT error: segment-based junction search failed with err =-6 what(): std::bad_alloc
0
gravatar for emily.grace.hamilton
2.5 years ago by
emily.grace.hamilton0 wrote:

Anyone know what this error is about or how to troubleshoot? Most samples work fine, but some throw this error.

[2016-05-25 15:33:36] Beginning TopHat run (v2.0.12)

[2016-05-25 15:33:36] Checking for Bowtie
Bowtie version: 2.1.0.0
[2016-05-25 15:33:36] Checking for Samtools
Samtools version: 0.1.12.a
[2016-05-25 15:33:36] Checking for Bowtie index files (transcriptome)..
[2016-05-25 15:33:36] Checking for Bowtie index files (genome)..
[2016-05-25 15:33:37] Checking for reference FASTA file
[2016-05-25 15:33:37] Generating SAM header for genome
[2016-05-25 15:33:41] Reading known junctions from GTF file
[2016-05-25 15:33:56] Pre-filtering multi-mapped left reads
[2016-05-25 15:33:56] Mapping HS1272_1 to genome genome with Bowtie2
[2016-05-25 18:50:17] Pre-filtering multi-mapped right reads
[2016-05-25 18:50:17] Mapping HS1272_2 to genome genome with Bowtie2
[2016-05-25 21:52:03] Preparing reads
left reads: min. length=50, max. length=50, 184399753 kept reads (5480263 discarded)
right reads: min. length=50, max. length=50, 183890626 kept reads (5989390 discarded)
[2016-05-25 23:29:45] Using pre-built transcriptome data..
[2016-05-25 23:30:12] Mapping left_kept_reads to transcriptome hg19_GencodeCompV19 with Bowtie2
[2016-05-26 02:54:18] Mapping right_kept_reads to transcriptome hg19_GencodeCompV19 with Bowtie2
[2016-05-26 05:59:29] Resuming TopHat pipeline with unmapped reads
[2016-05-26 06:42:55] Mapping left_kept_reads.m2g_um_seg1 to genome genome with Bowtie2 (1/2)
[2016-05-26 07:33:29] Mapping left_kept_reads.m2g_um_seg2 to genome genome with Bowtie2 (2/2)
[2016-05-26 08:54:12] Mapping right_kept_reads.m2g_um_seg1 to genome genome with Bowtie2 (1/2)
[2016-05-26 09:43:46] Mapping right_kept_reads.m2g_um_seg2 to genome genome with Bowtie2 (2/2)
[2016-05-26 10:20:18] Searching for junctions via segment mapping
[FAILED]
Error: segment-based junction search failed with err =-6
what(): std::bad_alloc

rna-seq tophat • 1.7k views
ADD COMMENTlink modified 2.5 years ago by Jennifer Hillman Jackson25k • written 2.5 years ago by emily.grace.hamilton0

This exception is usually thrown if the computer needs more RAM than available (http://www.cplusplus.com/reference/new/bad_alloc/).

ADD REPLYlink written 2.5 years ago by y.hoogstrate460
0
gravatar for Jennifer Hillman Jackson
2.5 years ago by
United States
Jennifer Hillman Jackson25k wrote:

Hello,

The error is coming from the underlying tool itself, the same as if run line-command. Double check the inputs carefully for format and content - both the fastq sequences and the custom reference genome.

Help: https://wiki.galaxyproject.org/Support

Best, Jen, Galaxy team

ADD COMMENTlink written 2.5 years ago by Jennifer Hillman Jackson25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 176 users visited in the last hour