working on slayer bug

[mothur.git] / distancedb.cpp

diff --git a/distancedb.cpp b/distancedb.cpp
index 8d0c6298578c4311a59407610ae28a87f48ee493..c1bf7e7fb08dc3fae5b3cf17f575580269677148 100644 (file)
--- a/distancedb.cpp
+++ b/distancedb.cpp
@@ -49,6 +49,7 @@ void DistanceDB::addSequence(Sequence seq) {
 vector<int> DistanceDB::findClosestSequences(Sequence* query, int numWanted){
        try {
                vector<int> topMatches;
+               Scores.clear();
                bool templateSameLength = true;
                string sequence = query->getAligned();
                vector<seqDist> dists; 
@@ -87,6 +88,7 @@ vector<int> DistanceDB::findClosestSequences(Sequence* query, int numWanted){
                                //fill topmatches with numwanted closest sequences indexes
                                for (int i = 0; i < numWanted; i++) {
                                        topMatches.push_back(dists[i].seq2);
+                                       Scores.push_back(dists[i].dist);
                                }
                        }else {
                                int bestIndex = 0;
@@ -103,6 +105,7 @@ vector<int> DistanceDB::findClosestSequences(Sequence* query, int numWanted){
                                }
                                searchScore = smallDist;
                                topMatches.push_back(bestIndex);
+                               Scores.push_back(smallDist);
                        }
                
                }else{