+ //cout << q->getName() << endl << q->getAligned() << endl << endl;
+ //cout << leftParent.getName() << '\t' << leftParent.getAligned().length() << endl << endl;
+ //cout << rightParent.getName() << '\t' << rightParent.getAligned().length() << endl << endl;
+ //cout << q.getName() << '\t' << q.getAligned().length() << endl << endl;
+ //cout << rightParent->getName() << endl << rightParent->getAligned() << endl << endl;
+ //cout << " length = " << rightParent->getAligned().length() << endl;
+
+ map<int, int> spots; //map from spot in original sequence to spot in filtered sequence for query and both parents
+ vector<data_struct> divs = runBellerophon(q, leftParent, rightParent, spots);
+
+ if (m->control_pressed) { return "no"; }
+// cout << "examining:\t" << refSeqs[i]->getName() << '\t' << refSeqs[j]->getName() << endl;