fixed ties issue in maligner getBlastSeqs

[mothur.git] / blastdb.cpp

diff --git a/blastdb.cpp b/blastdb.cpp
index 776955588abb474ae8f1123f4cb93ad8bbb73883..026c4ca3054f0d9db16eb53ec1327115f67239d4 100644 (file)
--- a/blastdb.cpp
+++ b/blastdb.cpp
@@ -116,6 +116,7 @@ vector<int> BlastDB::findClosestMegaBlast(Sequence* seq, int n) {
                vector<int> topMatches;
                
                ofstream queryFile;
+
                m->openOutputFile((queryFileName+seq->getName()), queryFile);
                queryFile << '>' << seq->getName() << endl;
                queryFile << seq->getUnaligned() << endl;
@@ -145,16 +146,17 @@ vector<int> BlastDB::findClosestMegaBlast(Sequence* seq, int n) {
                        
                        m->gobble(m8FileHandle);
                        topMatches.push_back(templateAccession);
+                       megaScores.push_back(searchScore);
 //cout << templateAccession << endl;
                }
                m8FileHandle.close();
                remove((queryFileName+seq->getName()).c_str());
                remove((blastFileName+seq->getName()).c_str());
-//cout << "\n\n" ;             
+//cout << "\n" ;               
                return topMatches;
        }
        catch(exception& e) {
-               m->errorOut(e, "BlastDB", "findClosest");
+               m->errorOut(e, "BlastDB", "findClosestMegaBlast");
                exit(1);
        }
 }