Remove gene pair only if is_self is on (#747)

tanghaibao · Jan 29, 2025 · 3730993 · 3730993
1 parent bd6be19
commit 3730993
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/jcvi/compara/blastfilter.py b/src/jcvi/compara/blastfilter.py
@@ -54,7 +54,7 @@ def blastfilter_main(blast_file, p, opts):
     nwarnings = 0
     for b in blasts:
         query, subject = b.query, b.subject
-        if query == subject:
+        if is_self and query == subject:
             continue
 
         if ostrip: