From d0cebdade6cfd0dbae6ba25d0f94508d0ac158d6 Mon Sep 17 00:00:00 2001 From: huhlim Date: Tue, 14 Jan 2020 17:09:33 -0500 Subject: [PATCH] accept sequence during filter2 when there is no overlap (cov_kj=0) --- src/hhalignment.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/hhalignment.cpp b/src/hhalignment.cpp index 2f2da73b..f45ff7b9 100755 --- a/src/hhalignment.cpp +++ b/src/hhalignment.cpp @@ -1916,7 +1916,7 @@ int Alignment::Filter2(char keep[], int coverage, int qid, float qsc, } //dissimilarity < acceptace threshold? Reject! - if (diff < diff_suff && float(diff) <= diff_min_frac * cov_kj) + if (diff < diff_suff && float(diff) < diff_min_frac * cov_kj) break; }