Segmentation fault (core dumped) during bwa mem mapping

Hi,

I ran bwa mem with trimmed fastq files (ERR2593198) but I saw following error:

bwa mem CHO-PICR.fasta ../2.ngsShort/trimmed_ERR2593198_1.fastq ../2.ngsShort/trimmed_ERR2593198_2.fastq

[M::bwa_idx_load_from_disk] read 0 ALT contigs
@PG ID:bwa PN:bwa VN:0.7.17-r1188 CL:../downloads/bwa-0.7.17/bwa mem CHO-PICR.fasta ../2.ngsShort/trimmed_ERR2593198_1.fastq ../2.ngsShort/trimmed_ERR2593198_2.fastq
[M::process] read 92156 sequences (10000179 bp)…
Segmentation fault (core dumped)

To figure out what’s happening, I ran gdb and the report was like this:

bwa mem CHO-PICR.fasta ../2.ngsShort/trimmed_ERR2593198_1.fastq ../2.ngsShort/trimmed_ERR2593198_2.fastq

[Thread debugging using libthread_db enabled]

Using host libthread_db library “/lib/x86_64-linux-gnu/libthread_db.so.1”.

[M::bwa_idx_load_from_disk] read 0 ALT contigs

[New Thread 0x7ffff7c2a700 (LWP 18232)]

[New Thread 0x7ffff7429700 (LWP 18233)]

[M::process] read 92156 sequences (10000179 bp)…

[New Thread 0x7ffff6052700 (LWP 18234)]

Thread 4 “bwa” received signal SIGSEGV, Segmentation fault.

[Switching to Thread 0x7ffff6052700 (LWP 18234)]

bwt_2occ4 (bwt=bwt@entry=0x5555555ad6e0, k=18446744073709551615, l=18446744073709551615,
cntk=cntk@entry=0x7ffff6051a70, cntl=cntl@entry=0x7ffff6051a90) at bwt.c:203

203 memcpy(cntk, p, 4 * sizeof(bwtint_t));

(gdb) bt

.#0 bwt_2occ4 (bwt=bwt@entry=0x5555555ad6e0, k=18446744073709551615, l=18446744073709551615,
cntk=cntk@entry=0x7ffff6051a70, cntl=cntl@entry=0x7ffff6051a90) at bwt.c:203

.#1 0x000055555557a56b in bwt_extend (bwt=bwt@entry=0x5555555ad6e0,
ik=ik@entry=0x7ffff6051b50, ok=ok@entry=0x7ffff6051ba0, is_back=is_back@entry=0)
at bwt.c:266

.#2 0x000055555557a843 in bwt_smem1a (bwt=bwt@entry=0x5555555ad6e0, len=len@entry=96,
q=q@entry=0x7ffff0013900 “010301”, x=x@entry=0, min_intv=min_intv@entry=1,
max_intv=max_intv@entry=0, mem=0x7ffff1e1afd8, tmpvec=0x7ffff1e1aff0) at bwt.c:310

.#3 0x000055555557adf6 in bwt_smem1 (bwt=bwt@entry=0x5555555ad6e0, len=len@entry=96,
q=q@entry=0x7ffff0013900 “010301”, x=x@entry=0, min_intv=min_intv@entry=1,
mem=mem@entry=0x7ffff1e1afd8, tmpvec=0x7ffff1e1aff0) at bwt.c:355

.#4 0x00005555555819fd in mem_collect_intv (a=0x7ffff1e1afc0,
seq=0x7ffff0013900 “010301”, len=96, bwt=0x5555555ad6e0, opt=0x5555555ad3b0)
at bwamem.c:123

.#5 mem_chain (opt=opt@entry=0x5555555ad3b0, bwt=0x5555555ad6e0, bns=bns@entry=0x5555555adb50,
len=len@entry=96, seq=seq@entry=0x7ffff0013900 “010301”, buf=0x7ffff1e1afc0)
at bwamem.c:264

.#6 0x000055555558fac4 in mem_align1_core (opt=0x5555555ad3b0, bwt=<optimized out>,
bns=0x5555555adb50, pac=0x5555555add90 “”, l_seq=96, seq=0x7ffff0013900 “010301”,
buf=0x7ffff1e1afc0) at bwamem.c:1060

.#7 0x000055555558fe56 in worker1 (data=0x7ffff7c29cb0, i=<optimized out>, tid=<optimized out>)
at bwamem.c:1181

.#8 0x00005555555766ca in ktf_worker (data=0x7ffff7c29bf0) at kthread.c:42

.#9 0x00007ffff7e33609 in start_thread (arg=<optimized out>) at pthread_create.c:477

.#10 0x00007ffff7d4f293 in clone () at ../sysdeps/unix/sysv/linux/x86_64/clone.S:95


Can anyone please help me how to fix this error?

Thank you.

Read more here: Source link