chimera.slayer

diff --git a/blastdb.cpp b/blastdb.cpp
index 8bf79c552c97e905f77011f2791a8f7f92a82172..aadff822bc9efb930f3a8f4c969fc3e5bc08e3fb 100644 (file)
--- a/blastdb.cpp
+++ b/blastdb.cpp
@@ -147,7 +147,7 @@ vector<int> BlastDB::findClosestMegaBlast(Sequence* seq, int n, int minPerID) {
                        //while (!m8FileHandle.eof())   {       char c = m8FileHandle.get(); if (c == 10 || c == 13){   break;  }else{ cout << c; }     } //
                                //cout << endl;
                        m->gobble(m8FileHandle);
-                       if (score >= minPerID) { //this follows broads CS, but probably should be searchScore...
+                       if (searchScore >= minPerID) { 
                                topMatches.push_back(templateAccession);
                                Scores.push_back(searchScore);
                        }