projects
/
samtools.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
| inline |
side by side
Mention ‘SAMtools’ in libbam-dev's description, to make it easier to find with apt...
[samtools.git]
/
kprobaln.c
diff --git
a/kprobaln.c
b/kprobaln.c
index 5201c1af75b9c24636ee67481b0461ab7713e7a1..894a2ae625f9e0b2b9d38337b8651d8f02bcce4a 100644
(file)
--- 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) Pr
1
+= -4.343 * log(p), p = 1.;
}
Pr1 += -4.343 * log(p * l_ref * l_query);
Pr = (int)(Pr1 + .499);