Control-FREEC - Segmentation fault (core dumped)
0
0
Entering edit mode
4.5 years ago
igor 12k

I am running Control-FREEC. I run into this error:

CG-content printed into ./GC_profile.targetedRegions.cnp
..using GC-content to normalize the control profile
..will remove all windows with read count in the control less than 50
Warning: control length is not equal to the sample length for chromosome 1
Warning: control length is not equal to the sample length for chromosome 2
Warning: control length is not equal to the sample length for chromosome 3
Warning: control length is not equal to the sample length for chromosome 4
Warning: control length is not equal to the sample length for chromosome 5
Warning: control length is not equal to the sample length for chromosome 6
Warning: control length is not equal to the sample length for chromosome 7
Warning: control length is not equal to the sample length for chromosome 8
Warning: control length is not equal to the sample length for chromosome 9
Warning: control length is not equal to the sample length for chromosome 10
Warning: control length is not equal to the sample length for chromosome 11
Warning: control length is not equal to the sample length for chromosome 12
Warning: control length is not equal to the sample length for chromosome 13
Warning: control length is not equal to the sample length for chromosome 14
Warning: control length is not equal to the sample length for chromosome 15
Warning: control length is not equal to the sample length for chromosome 16
Warning: control length is not equal to the sample length for chromosome 17
Warning: control length is not equal to the sample length for chromosome 18
Warning: control length is not equal to the sample length for chromosome 19
Warning: control length is not equal to the sample length for chromosome 20
Warning: control length is not equal to the sample length for chromosome 21
Warning: control length is not equal to the sample length for chromosome 22
Warning: control length is not equal to the sample length for chromosome X
Segmentation fault (core dumped)


I also tried compiling with a newer version of gcc (6.1 instead of 4.4), which produces a slightly different error message:

Warning: control length is not equal to the sample length for chromosome X
terminate called after throwing an instance of 'std::bad_alloc'
Aborted (core dumped)


Usually, a core dump happens when there is a memory issue, but it's only 45MB, so that can't be the reason.

If I repeat the analysis, but change from mateFile using BAM files to mateCopyNumberFile using the .cpn files from the failed analysis, the problem goes away:

CG-content printed into ./GC_profile.targetedRegions.cnp
..using GC-content to normalize the control profile
..will remove all windows with read count in the control less than 50
..will process the control file as well: removing all windows with read count in the control less than 50
..Set ploidy for the control genome equal to 2
..Running FREEC with ploidy set to 2


It's using the same .cpn files for doing that calculation as far as I can tell. Why does it get stuck depending on input? I tried multiple times with the same results.

This is based on the latest Control-FREEC 10.6. I used the same input files for older version 7.2 and that one works.

A similar problem was reported previously:

Unfortunately, those discussions don't actually have a solution posted. Has anyone run into this problem and solved it?

cnv cna freec control-freec • 2.9k views