fixed minor bug where some phages were not reported

SushiLab · Jan 19, 2022 · bbbc9f6 · bbbc9f6
1 parent 7911534
commit bbbc9f6
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/mVIRs/mvirs.py b/mVIRs/mvirs.py
@@ -957,7 +957,7 @@ def extract_regions(clipped_file, opr_file, reference_fasta_file, output_fasta_f
         if minsize >= length or length >= maxsize:
             continue
         # 2:
-        if opr_cnt < minoprcount or hs_cnt < minhscount or (opr_cnt + hs_cnt) <= mincombcount:
+        if opr_cnt < minoprcount or hs_cnt < minhscount or (opr_cnt + hs_cnt) < mincombcount:
             continue
         # 3:
         scaffold_length = len(reference_header_2_sequence[scaffold])