v1.0.8

fixed an edge case when blast output contains a lot of hits with no taxid
GDKO · Dec 18, 2024 · 98c1029 · 98c1029
1 parent ed4e038
commit 98c1029
Show file tree

Hide file tree

Showing 4 changed files with 5 additions and 5 deletions.
diff --git a/aux_scripts/calculate_ai.py b/aux_scripts/calculate_ai.py
@@ -101,7 +101,7 @@ def main():
                         if gene not in best_hit_ingroup.keys():
                             best_hit_ingroup[gene] = {}
                             best_hit_ingroup[gene]["hit"] = hit
-                            best_hit_ingroup[gene]["pos"] = str(num_hits[gene])
+                            best_hit_ingroup[gene]["pos"] = str(num_hits[gene] + skip)
                             best_hit_ingroup[gene]["iden"] = iden
                             best_hit_ingroup[gene]["evalue"] = str(evalue)
                             best_hit_ingroup[gene]["bitscore"] = bitscore
@@ -111,7 +111,7 @@ def main():
                         if gene not in best_hit_donor.keys():
                             best_hit_donor[gene] = {}
                             best_hit_donor[gene]["hit"] = hit
-                            best_hit_donor[gene]["pos"] = str(num_hits[gene])
+                            best_hit_donor[gene]["pos"] = str(num_hits[gene] + skip)
                             best_hit_donor[gene]["iden"] = iden
                             best_hit_donor[gene]["evalue"] = str(evalue)
                             best_hit_donor[gene]["bitscore"] = bitscore

diff --git a/avp b/avp
@@ -1,5 +1,5 @@
 #!/usr/bin/env python3
-#v1.0.7
+#v1.0.8
 
 from depot.interface import main
 

diff --git a/depot/interface.py b/depot/interface.py
@@ -23,7 +23,7 @@
 from docopt import docopt
 
 def main():
-    args = docopt(__doc__,version='1.0.7', options_first=True)
+    args = docopt(__doc__,version='1.0.8', options_first=True)
     if args['<module>'] == 'prepare':
         import depot.prepare as prepare
         prepare.main()

diff --git a/depot/prepare.py b/depot/prepare.py
@@ -193,7 +193,7 @@ def main():
                 if len(query_id) > 255:
                     sys.exit("[x] Query names are too long")
                 if query_id in queries_info.keys(): # Queries that pass the initial selection
-                    if(len(query_dict_set[query_id]) <= queries_info[query_id]["pos"]):
+                    if(len(query_dict_set[query_id]) < queries_info[query_id]["pos"]):
                         query_hit_id = L_hitqline[1]
                         if "@" in query_id:
                             sys.exit("@ symbol is not allowed: " + query_id)