- //require at least 3 SNPs on each side of the break
- if ((numSNPSLeft >= 3) && (numSNPSRight >= 3)) {
-
- //removed in 12/09 version of chimeraSlayer
- //int winSizeLeft = divs[k].winLEnd - divs[k].winLStart + 1;
- //int winSizeRight = divs[k].winREnd - divs[k].winRStart + 1;
-
- //float snpRateLeft = numSNPSLeft / (float) winSizeLeft;
- //float snpRateRight = numSNPSRight / (float) winSizeRight;
- //float logR = log(snpRateLeft / snpRateRight) / log(2.0);
-
- // do not accept excess snp ratio on either side of the break
- //if (abs(logR) < 1 ) {
+// cout << numSNPSLeft << '\t' << numSNPSRight << endl;
+ //require at least 4 SNPs on each side of the break
+ if ((numSNPSLeft >= 4) && (numSNPSRight >= 4)) {