+//cout << "not equal " << Q << '\t' << A << '\t' << B << endl;
+
+ //ensure not neighboring a gap. change to 12/09 release of chimeraSlayer - not sure what this adds, but it eliminates alot of SNPS
+ if (
+ //did query loose a base here during filter??
+ ( i == 0 || abs (baseSpots[0][i] - baseSpots[0][i-1]) == 1) &&
+ ( i == query.length() || abs (baseSpots[0][i] - baseSpots[0][i+1]) == 1)
+ &&
+ //did parentA loose a base here during filter??
+ ( i == 0 || abs (baseSpots[1][i] - baseSpots[1][i-1]) == 1) &&
+ ( i == parentA.length() || abs (baseSpots[1][i] - baseSpots[1][i+1]) == 1)
+ &&
+ //did parentB loose a base here during filter??
+ ( i == 0 || abs (baseSpots[2][i] - baseSpots[2][i-1]) == 1) &&
+ ( i == parentB.length() || abs (baseSpots[2][i] - baseSpots[2][i+1]) == 1)
+ )
+ {