small change in hit rate calculation in cyclecheck modul

soedinglab · Feb 6, 2020 · a66c02c · a66c02c
1 parent e54e2e3
commit a66c02c
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/assembler/cyclecheck.cpp b/src/assembler/cyclecheck.cpp
@@ -171,7 +171,8 @@ int cyclecheck(int argc, const char **argv, const Command& command) {
                         unsigned int diagbandHits = 0;
 
                         for (size_t i = lower; i <= upper; i++) {
-                            diagbandHits += diagHits[i];
+                            if (diagHits[i] <= diagHits[d])
+                                diagbandHits += diagHits[i];
                         }
 
                         float diagbandHitRate = static_cast<float>(diagbandHits) / (diaglen - kmerSize + 1);