Software Crash - Olympia oyster genome assembly with Masurca on Mox

Ah, the joys of bioinformatics. I just received an email from Mox indicating that [the Masurca assembly I started 11 DAYS AGO (!!)(2017/10/19/genome-assembly-olympia-oyster-illumina-pacbio-reads-using-masurca.html) crashed.


I’m probably not going to put much effort in to trying to figure out what went wrong, but here’s some log file snippets for reference. I’ll probably drop a line to the developers and see if they have any easy ways to address whatever caused the problems, but that’s about as much effort as I’m willing to put into troubleshooting this assembly.

Additionally, since this crashed, I’m not going to bother moving any of the files off of Mox. That means they will be deleted automatically by the system around Nov. 9th, 2017.

slurm-94620.out (tail)

[code lang=text] compute_psa 6601202 2632582819 Refining alignments Joining Generating assembly input files Coverage of the mega-reads less than 5 – using the super reads as well Coverage threshold for splitting unitigs is 138 minimum ovl 63 Running assembly /gscratch/srlab/programs/MaSuRCA-3.2.3/bin/ line 85: 24330 Aborted (core dumped) overlapStoreBuild -o $ASM_DIR/$ASM_PREFIX.ovlStore -M 65536 -g $ASM_DIR/$ASM_PREFIX.gkpStore $ASM_DIR/overlaps_dedup.ovb.gz > $ASM_DIR/overlapStore.rebuild.err 2>&1 Assembly stopped or failed, see [Mon Oct 30 23:19:37 PDT 2017] Assembly stopped or failed, see [/code] (tail)

[code lang=text] number of threads = 28 (OpenMP default)

ERROR: overlapStore ‘/gscratch/scrubbed/samwhite/20171019_masurca_oly_assembly/’ is incomplete; previous overlapStoreBuild probably crashed.

Failure message:

failed to unitig [/code]


[code lang=text] Scanning overlap files to count the number of overlaps. Found 277.972 million overlaps. Memory limit 65536MB supplied. Ill put 3246167525 IIDs (3435.97 million overlaps) into each of 1 buckets. bucketizing bucketizing DONE! overlaps skipped: 0 OBT - low quality 0 DUP - non-duplicate overlap 0 DUP - different library 0 DUP - dedup not requested terminate called after throwing an instance of std::bad_alloc what(): std::bad_alloc

Failed with Aborted

Backtrace (mangled):

overlapStoreBuild[0x40523a] /usr/lib64/[0x2af83b3c0100] /usr/lib64/[0x2af83c0395f7] /usr/lib64/[0x2af83c03ace8] /usr/lib64/[0x2af83b62d9d5] /usr/lib64/[0x2af83b62b946] /usr/lib64/[0x2af83b62b973] /usr/lib64/[0x2af83b62bb93] /usr/lib64/[0x2af83b62c12d] /usr/lib64/[0x2af83b62c1c9] overlapStoreBuild[0x402e10] /usr/lib64/[0x2af83c025b15] overlapStoreBuild[0x403089]

Backtrace (demangled):

[0] overlapStoreBuild() [0x40523a] [1] /usr/lib64/ + 0xf100 [0x2af83b3c0100] [2] /usr/lib64/ + 0x37 [0x2af83c0395f7] [3] /usr/lib64/ + 0x148 [0x2af83c03ace8] [4] /usr/lib64/ + 0x165 [0x2af83b62d9d5] [5] /usr/lib64/ + 0x5e946 [0x2af83b62b946] [6] /usr/lib64/ + 0x5e973 [0x2af83b62b973] [7] /usr/lib64/ + 0x5eb93 [0x2af83b62bb93] [8] /usr/lib64/ new(unsigned long) + 0x7d [0x2af83b62c12d] [9] /usr/lib64/ new + 0x9 [0x2af83b62c1c9] [10] overlapStoreBuild() [0x402e10] [11] /usr/lib64/ + 0xf5 [0x2af83c025b15] [12] overlapStoreBuild() [0x403089]

GDB: [/code]