X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=distancedb.cpp;h=ca6ffe8ba61217ac0b66b353e8f33461f4916e52;hb=050220fe7822cc660615972a0054cf4a83eefbe4;hp=5c49a936e9293c629e05ea5565f334a43e5cc37c;hpb=1d18c13c5d71e92ee850a0d8d7cda4f0ef530bf7;p=mothur.git

diff --git a/distancedb.cpp b/distancedb.cpp
index 5c49a93..ca6ffe8 100644
--- a/distancedb.cpp
+++ b/distancedb.cpp
@@ -48,6 +48,8 @@ vector<int> DistanceDB::findClosestSequences(Sequence* query, int numWanted){
 		bool templateSameLength = true;
 		string sequence = query->getAligned();
 		vector<seqDist> dists;
+		
+		searchScore = -1.0;
 	
 		if (numWanted > data.size()) { m->mothurOut("numwanted is larger than the number of template sequences, using "+ toString(data.size()) + "."); m->mothurOutEndLine(); numWanted = data.size(); }
 		
@@ -66,6 +68,9 @@ vector<int> DistanceDB::findClosestSequences(Sequence* query, int numWanted){
 			
 			sort(dists.begin(), dists.end(), compareSequenceDistance);  //sorts by distance lowest to highest
 			
+			//save distance of best match
+			searchScore = dists[0].dist;
+			
 			//fill topmatches with numwanted closest sequences indexes
 			for (int i = 0; i < numWanted; i++) {
 				topMatches.push_back(dists[i].seq2);