  
  
  #ifdef  _DEBUG
-@@ -764,6 +765,18 @@ Return Value:
+@@ -764,6 +766,19 @@ Return Value:
                  primaryResult->score = bestScore;
                  if (bestScore <= maxK) {
                      primaryResult->location = bestScoreGenomeLocation;

+			diffbufPair[2] = '\0';
+		}
+		ztzSeparator = "\tZT:Z:";
+		snprintf(ztzbuf, 1024, "%d,%s,%s,%d,%d,%d,%lld,%d,%d",
+				 -features.bestScore,
+				 diffbuf,
+				 diffbufPair,
+				 (int)(9999.0 * features.probabilityOfBestCandidate /
+				       (features.probabilityOfAllCandidates > 0 ? features.probabilityOfAllCandidates : 1)),
+				 features.popularSeedsSkipped,
+				 features.popularSeedsSkippedPair,
+				 features.countOfHashTableLookups,

         *spaceUsed = bamSize;
diff -rupN -x Xcode -x .git tmp/snap/SNAPLib/BaseAligner.cpp snap/SNAPLib/BaseAligner.cpp
--- tmp/snap/SNAPLib/BaseAligner.cpp	2016-01-18 14:51:10.000000000 -0500
+++ snap/SNAPLib/BaseAligner.cpp	2016-08-22 16:55:46.000000000 -0400
@@ -289,7 +289,7 @@ Return Value:
     primaryResult->score = UnusedScoreValue;
     primaryResult->status = NotFound;
