tracking down chimera.slayer issue

[mothur.git] / maligner.cpp

diff --git a/maligner.cpp b/maligner.cpp
index 5ea41fcd867d291fe928f096cb7d90e03d4d19d4..7d7145f4ca6bd9fac7a316745ace4432e4c26fe0 100644 (file)
--- a/maligner.cpp
+++ b/maligner.cpp
@@ -733,11 +733,13 @@ vector<Sequence*> Maligner::getBlastSeqs(Sequence* q, int num) {
 //cout << q->getName() << " merged results size = " << mergedResults.size() << '\t' << "numwanted = " << numWanted <<  endl;           
                for (int i = 0; i < numWanted; i++) {
 //cout << db[mergedResults[i]]->getName()  << '\t' << mergedResults[i] << endl;        
-                       //if (db[mergedResults[i]]->getName() != q->getName()) { 
+                       
+                       if (db[mergedResults[i]]->getName() != q->getName()) { 
                                Sequence* temp = new Sequence(db[mergedResults[i]]->getName(), db[mergedResults[i]]->getAligned());
                                refResults.push_back(temp);
                                indexes.push_back(mergedResults[i]);
-                       //}
+                       }
+                       
 //cout << mergedResults[i] << endl;
                }
 //cout << "done " << q->getName()  << endl;