X-Git-Url: http://woldlab.caltech.edu/gitweb/?p=samtools.git;a=blobdiff_plain;f=kprobaln.c;h=894a2ae625f9e0b2b9d38337b8651d8f02bcce4a;hp=5201c1af75b9c24636ee67481b0461ab7713e7a1;hb=d59366b97f0e6a0cbb0dfcdc5679f0ed88c8fe1e;hpb=016b50ab60e879a0b8f81cb76ce11ea360a03d4a

diff --git a/kprobaln.c b/kprobaln.c
index 5201c1a..894a2ae 100644
--- a/kprobaln.c
+++ b/kprobaln.c
@@ -161,7 +161,7 @@ int kpa_glocal(const uint8_t *_ref, int l_ref, const uint8_t *_query, int l_quer
 		double p = 1., Pr1 = 0.;
 		for (i = 0; i <= l_query + 1; ++i) {
 			p *= s[i];
-			if (p < 1e-100) Pr += -4.343 * log(p), p = 1.;
+			if (p < 1e-100) Pr1 += -4.343 * log(p), p = 1.;
 		}
 		Pr1 += -4.343 * log(p * l_ref * l_query);
 		Pr = (int)(Pr1 + .499);