Commit 8d259859 authored by Mikaël Salson's avatar Mikaël Salson

filter.cpp: Threshold for limit number of kmers must be soft

parent bb37a769
......@@ -114,7 +114,7 @@ BioReader FilterWithACAutomaton::filterBioReaderWithACAutomaton(
if(!element.first.isGeneric()){
continue;
}
if(nbKmers == kmer_threshold && nb_kmers_limit == element.second){
if(nbKmers == kmer_threshold && nb_kmers_limit <= element.second){
// We have reached our limit of number of genes recovered but we
// continue taking sequences are they have a similar number of
// matching k-mers.
......
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