Make the option min-score and min-length work in centrifuge-kreport w…

…hen using lca.
DaehwanKimLab · May 10, 2020 · 7ba4885 · 7ba4885
1 parent 6cc874e
commit 7ba4885
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/centrifuge-kreport b/centrifuge-kreport
@@ -97,11 +97,12 @@ if ($is_cnts_table) {
     my $queryLength = $cols[ $headerMap{ "queryLength" } ] ; 
     my $numMatches = $cols[ $headerMap{ "numMatches" } ] ; 
 
+    next if defined $min_length && $hitLength < $min_length;
+    next if defined $min_score && $score < $min_score;
+
     $taxID = 1 if ( !isTaxIDInTree( $taxID ) ) ;
 
     if ($no_lca) {
-      next if defined $min_length && $hitLength < $min_length;
-      next if defined $min_score && $score < $min_score;
       $taxo_counts{$taxID} += 1/$numMatches;
       $seq_count += 1/$numMatches;
     } else {