Commit 38670846 authored by Mikaël Salson's avatar Mikaël Salson
Browse files

Representative: We were always counting one kmer in the representative

parent eb7edce3
......@@ -131,10 +131,10 @@ void KmerRepresentativeComputer::compute() {
size_t length_run = 0;
size_t i = pos_required;
if (pos_required)
do {
while (i > 0 && isSufficienlyExpressed(counts[i-1].count, max)) {
i--;
length_run++;
} while (i > 0 && isSufficienlyExpressed(counts[i].count, max));
}
for (i = pos_required; i < counts.size(); i++) {
while (i < counts.size() &&
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment