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.hamilton • 0 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 COMMENT • link • modified 2.5 years ago by Jennifer Hillman Jackson ♦ 25k • written 2.5 years ago by emily.grace.hamilton • 0

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

ADD REPLY • link written 2.5 years ago by y.hoogstrate • 460
0
gravatar for Jennifer Hillman Jackson
2.5 years ago by
United States
Jennifer Hillman Jackson ♦ 25k 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 COMMENT • link written 2.5 years ago by Jennifer Hillman Jackson ♦ 25k
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