- if (divs[k].bsMax >= (minBS-10)) {
- snpsLeft = getSNPS(divs[k].parentA.getAligned(), divs[k].querySeq.getAligned(), divs[k].parentB.getAligned(), divs[k].winLStart, divs[k].winLEnd);
- snpsRight = getSNPS(divs[k].parentA.getAligned(), divs[k].querySeq.getAligned(), divs[k].parentB.getAligned(), divs[k].winRStart, divs[k].winREnd);
-
- if (m->control_pressed) { delete q; delete leftParent; delete rightParent; return "no"; }
-
- numSNPSLeft = snpsLeft.size();
- numSNPSRight = snpsRight.size();
-
- //require at least 4 SNPs on each side of the break
- if ((numSNPSLeft >= 4) && (numSNPSRight >= 4)) {
- float BS_A, BS_B;
- bootstrapSNPS(snpsLeft, snpsRight, BS_A, BS_B, 1000);
+ if ((divs[k].bsMax >= (minBS-10)) && (iters < 1000)) {
+ bootstrapSNPS(snpsLeft, snpsRight, BS_A, BS_B, 1000);