SeekSubclones::getSegmentsInfo SegFault at Line 83

[deathstar550]

Line 83 in SeekSubclones.cpp:

if (breakpoints.back()<length) {breakpoints.push_back(length);}

If breakpoints is empty, breakpoints.back() throws a segfault. This triggered on some of my data (out of about 100+ samples, two of the samples had one chromosome where this occurred).

The above needs to be rewritten to include handling of the special case when breakpoints is empty:

        if (breakpoints.empty()) {
            breakpoints.push_back(length);
        } else if (breakpoints.back()<length) {
            breakpoints.push_back(length);
        }

This solved the issue for me.