Dear all,
I run CeleraAssembler using 60 G bp next-sequencing data, and I come into and error when building overlap store after 1-overlapper. I am wondering if I have run out of memory. Are there any recommendations for changing the parameters? I have also attach my spec_file below. This is part of my log file:
----------------------------------------START Mon Dec 6 20:30:41 2021 /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/gatekeeper -o /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.gkpStore.BUILDING /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/Tbrevi-01F_notCombined.frg /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/Tbrevi-01F_extendedFrags.frg > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.gkpStore.err 2>&1 ----------------------------------------END Mon Dec 6 20:57:44 2021 (1623 seconds) numFrags = 341784124 No need to run meryl for OBT (OBT is disabled). ----------------------------------------START Mon Dec 6 20:57:44 2021 /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/meryl -B -C -v -m 22 -memory 320000 -threads 60 -c 0 -L 2 -s /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.gkpStore:chain -o /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./0-mercounts/Tbrevi-01F-C-ms22-cm0 > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./0-mercounts/meryl.err 2>&1 ----------------------------------------END Mon Dec 6 22:44:57 2021 (6433 seconds) ----------------------------------------START Mon Dec 6 22:44:57 2021 /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/estimate-mer-threshold -m /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./0-mercounts/Tbrevi-01F-C-ms22-cm0 > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./0-mercounts/Tbrevi-01F-C-ms22-cm0.estMerThresh.out 2> /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./0-mercounts/Tbrevi-01F-C-ms22-cm0.estMerThresh.err ----------------------------------------END Mon Dec 6 22:44:57 2021 (0 seconds) ----------------------------------------START Mon Dec 6 22:44:57 2021 /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/meryl -Dt -n 5134 -s /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./0-mercounts/Tbrevi-01F-C-ms22-cm0 > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./0-mercounts/Tbrevi-01F.nmers.ovl.fasta 2> /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./0-mercounts/Tbrevi-01F.nmers.ovl.fasta.err ----------------------------------------END Mon Dec 6 22:45:28 2021 (31 seconds) Reset OVL mer threshold from to 5134. ----------------------------------------START Mon Dec 6 22:45:29 2021 /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/overlap_partition \ -g /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.gkpStore \ -bl 1000000000 \ -bs 0 \ -rs 24000000 \ -rl 0 \ -o /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper \ > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/overlap_partition.err 2>&1 ----------------------------------------END Mon Dec 6 22:45:46 2021 (17 seconds) Created 462 overlap jobs. Last batch '001', last job '000462'. ----------------------------------------START CONCURRENT Mon Dec 6 22:45:46 2021 /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/overlap.sh 1 > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/000001.out 2>&1 /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/overlap.sh 2 > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/000002.out 2>&1 ...omitting about 400 lines here /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/overlap.sh 461 > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/000461.out 2>&1 /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/overlap.sh 462 > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/000462.out 2>&1 ----------------------------------------END CONCURRENT Tue Dec 7 09:04:31 2021 (37125 seconds) ----------------------------------------START Tue Dec 7 09:04:36 2021 find -L /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper \( -name \*ovb.gz -or -name \*ovb \) -print > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.ovlStore.list ----------------------------------------END Tue Dec 7 09:04:41 2021 (5 seconds) ----------------------------------------START Tue Dec 7 09:04:41 2021 /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/overlapStoreBuild -o /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.ovlStore.BUILDING -g /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.gkpStore -M 500000 -L /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.ovlStore.list > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.ovlStore.err 2>&1 sh: line 1: 77876 Aborted /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/overlapStoreBuild -o /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.ovlStore.BUILDING -g /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.gkpStore -M 500000 -L /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.ovlStore.list > /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.ovlStore.err 2>&1 ----------------------------------------END Tue Dec 7 16:05:22 2021 (25241 seconds) ERROR: Failed with signal ABRT (6) ================================================================================ runCA failed. ---------------------------------------- Stack trace: at /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/runCA line 1628, <J> line 462. main::caFailure("failed to create the overlap store", "/public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-"...) called at /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/runCA line 3696 main::createOverlapStore() called at /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/runCA line 6255 ---------------------------------------- Last few lines of the relevant log file (/public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./Tbrevi-01F.ovlStore.err): bucketizing /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/001/000349.ovb.gz bucketizing /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/001/000127.ovb.gz bucketizing /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/001/000099.ovb.gz bucketizing /public1/home/sc61333/project/gaoyuanqiu/Tbrevi-01F/./1-overlapper/001/000347.ovb.gz 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): /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/overlapStoreBuild(_Z17AS_UTL_catchCrashiP9siginfo_tPv+0x2a)[0x40945a] /lib64/libpthread.so.0(+0xf630)[0x2b71ec0ec630] /lib64/libc.so.6(gsignal+0x37)[0x2b71ec32f387] /lib64/libc.so.6(abort+0x148)[0x2b71ec330a78] /lib64/libstdc++.so.6(_ZN9__gnu_cxx27__verbose_terminate_handlerEv+0x165)[0x2b71eb6f7a95] /lib64/libstdc++.so.6(+0x5ea06)[0x2b71eb6f5a06] /lib64/libstdc++.so.6(+0x5ea33)[0x2b71eb6f5a33] /lib64/libstdc++.so.6(+0x5ec53)[0x2b71eb6f5c53] /lib64/libstdc++.so.6(_Znwm+0x7d)[0x2b71eb6f61ed] /lib64/libstdc++.so.6(_Znam+0x9)[0x2b71eb6f6289] /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/overlapStoreBuild(main+0x1450)[0x4067a0] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b71ec31b555] /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/overlapStoreBuild[0x406a8d] Backtrace (demangled): [0] /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/overlapStoreBuild::AS_UTL_catchCrash(int, siginfo_t*, void*) + 0x2a [0x40945a] [1] /lib64/libpthread.so.0::(null) + 0xf630 [0x2b71ec0ec630] [2] /lib64/libc.so.6::(null) + 0x37 [0x2b71ec32f387] [3] /lib64/libc.so.6::(null) + 0x148 [0x2b71ec330a78] [4] /lib64/libstdc++.so.6::__gnu_cxx::__verbose_terminate_handler() + 0x165 [0x2b71eb6f7a95] [5] /lib64/libstdc++.so.6::(null) + 0x5ea06 [0x2b71eb6f5a06] [6] /lib64/libstdc++.so.6::(null) + 0x5ea33 [0x2b71eb6f5a33] [7] /lib64/libstdc++.so.6::(null) + 0x5ec53 [0x2b71eb6f5c53] [8] /lib64/libstdc++.so.6::operator new(unsigned long) + 0x7d [0x2b71eb6f61ed] [9] /lib64/libstdc++.so.6::operator new[](unsigned long) + 0x9 [0x2b71eb6f6289] [10] /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/overlapStoreBuild::(null) + 0x1450 [0x4067a0] [11] /lib64/libc.so.6::(null) + 0xf5 [0x2b71ec31b555] [12] /public1/home/sc61333/software/wgs-8.3rc2/Linux-amd64/bin/overlapStoreBuild() [0x406a8d] GDB: ---------------------------------------- Failure message: failed to create the overlap store
This is my spec_file for running runAC:
ovlConcurrency=7 ovlThreads=12 cnsConcurrency =60 merSize=22 merylMemory=320000 merylThreads=60 # merThreshold=5000 doOBT=0 overlapper=ovl ovlStoreMemory = 500000 ovlRefBlockSize = 24000000 ovlHashBits = 26 ovlHashBlockLength = 1000000000 doFragmentCorrection=0 unitigger=bogart batMemory=320 batThreads = 60 doExtendClearRanges=0 doToggle = 0 gkpFixInsertSizes=0 Tbrevi-01F_notCombined.frg Tbrevi-01F_extendedFrags.frg