X-Git-Url: http://woldlab.caltech.edu/gitweb/?p=pysam.git;a=blobdiff_plain;f=samtools%2Fkprobaln.c.pysam.c;h=3404936c70689012234e956e754844f1b5d0618e;hp=b87dcf07c831dc946900028fc1a2b3a073b404b5;hb=ca46ef4ba4a883c57cea62d5bf1bc021f1185109;hpb=d02fe5283ed7a93a2f76a5d6dc6e37b40c11b9b1

diff --git a/samtools/kprobaln.c.pysam.c b/samtools/kprobaln.c.pysam.c
index b87dcf0..3404936 100644
--- a/samtools/kprobaln.c.pysam.c
+++ b/samtools/kprobaln.c.pysam.c
@@ -163,7 +163,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);