Error thrown during split-read assembly

[eeshaanrehani]

Hello,

I'm running Delly on a large WGS paired tumor-normal sample (from SEQC2) and ran into this error during the split-read assembly step of delly call:

terminate called after throwing an instance of 'std::length_error'
what():  basic_string::_M_create
./delly_no_small.sh: line 17: 1626942 Aborted                 (core dumped) delly call -n 1000 -o seqc2.temp.delly.bcf -g ${REF} ${SEQC2_T} ${SEQC2_N}

The specific BAM files I am using are tumor and normal, and I'm using the GDC's GRCh38.d1.vd1 reference.

Interestingly, the issue goes away when I remove the -n 1000 from the command. I'm using Delly v1.2.6, installed from bioconda.

Thanks!

Eeshaan

[sourav8888]

Hi Tobias,

I do have the same issue with delly run. Although we are working on Plasmodium (haploid).
Please let us know is there any way to resolve this issue.

Best,
Sourav.