Fix Filter2 bug of HH-suite in MMseqs2 soedinglab/hh-suite#182

soedinglab · Jan 16, 2020 · 9bcff28 · 9bcff28
1 parent 01db79d
commit 9bcff28
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/alignment/MsaFilter.cpp b/src/alignment/MsaFilter.cpp
@@ -395,7 +395,7 @@ void MsaFilter::filter(const int N_in, const int L, const int coverage, const in
 //            // DEBUG
 //            printf("%20.20s with %20.20s:  diff=%i  diff_min_frac*cov_kj=%f  diff_suff=%i  nres=%i  cov_kj=%i\n",sname[k],sname[j],diff,diff_min_frac*cov_kj,diff_suff,nres[k],cov_kj);
 //            printf("%s\n%s\n\n",seq[k],seq[j]);
-                if (diff < diff_suff && float(diff) <= diff_min_frac * cov_kj)
+                if (diff < diff_suff && float(diff) <= diff_min_frac * cov_kj && cov_kj > 0)
                     break;  //dissimilarity < acceptace threshold? Reject!
 
             }