bacterial-annotator 0.9.3 → 0.9.4

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: e6b26cd37618fedd557b8d9f7fe2c27a3b95dee5404d7751a8fb10ddcb63edc8
-  data.tar.gz: 2129f8372c3a48654cf3bdf6b75eccce3148b45f0e15191f3a555ff0769f9bd4
+  metadata.gz: 6873b45085524788e9ea2c82b29c5e3cfb6998c6cd7ff1f16d56dc64378c910e
+  data.tar.gz: 7b02f34f8ffef086ade18526db05b5fffdeabee339a4808ec09181a7cb108b53
 SHA512:
-  metadata.gz: 2400660df26bab507b5f9a3310105514d415a4950a06fdb7a3cbc9a4de2b4268d77c3725a13507e19c92fcdcc688f2bda4184aeadb88afd4a444addce1117fb9
-  data.tar.gz: 9caef8af954e6de43ed3eb2cd5a729065bafd2451aa5c73e00e60b400cb849666a53fb1e39d549f4b33925edbfea33f209a966da7f64d91ac1aa3d8b75e5f4d2
+  metadata.gz: 4b01f90fe7bc53e79a4245fbd64830a8fa967447e4a609715feef1ec21a68583531b9d1120631234204215bfa3bf6cfdd0b311704f030cf4cdf08d670fcabcf8
+  data.tar.gz: 15f9a395c37249bf937052d1917c89a42122a720bb1beddd1c4b2bddef3883111cb7e2e883da7757ccd9170eb3675499bbb70ac8dcfc62d4e722b9434549c92c

data/lib/bacterial-annotator/sequence-synteny.rb

@@ -1,10 +1,10 @@
 # -*- coding: utf-8 -*-
 # author:  	maxime déraspe
 # email:	maximilien1er@gmail.com
-# review:  	
+# review:
 # date:    	15-02-24
 # version: 	0.0.1
-# licence:  	
+# licence:
 
 require 'json'
 require 'zlib'
@@ -111,7 +111,7 @@ class SequenceSynteny
         system("#{@root}/diamond.linux makedb --db #{subject_file} --in #{subject_file} > /dev/null 2>&1")
         db_file = subject_file
       end
-      system("#{@root}/diamond.linux blastp --db #{db_file} -q #{query_file} -o #{@outdir}/#{@name}.diamond.tsv -f 6 > /dev/null 2>&1")
+      system("#{@root}/diamond.linux blastp --masking none --db #{db_file} -q #{query_file} -o #{@outdir}/#{@name}.diamond.tsv -f 6 > /dev/null 2>&1")
     else
       # system("#{@root}/glsearch36.linux -b 3 -E 1e-25 -m 8 #{@subject_file} #{@query_file} > #{@outdir}/#{@name}.fasta36.tsv")
     end

data/lib/bacterial-comparator.rb

@@ -45,6 +45,9 @@ class BacterialComparator
       min_pid = min_pid/100
     end
 
+    @min_cov = min_cov
+    @min_pid = min_pid
+
     @aln_opt = options[:align].downcase
     @run_phylo = 0
     if options[:phylogeny] == 1
@@ -592,8 +595,8 @@ class BacterialComparator
                                            query_prot_file,
                                            ref_prot_file,
                                            "Prot-Ref",
-                                           0.80,
-                                           0.80,
+                                           @min_cov,
+                                           @min_cov,
                                            "prot")
 
     print "# Running alignment with Reference Genome CDS (diamond).."
@@ -613,6 +616,7 @@ class BacterialComparator
       next if ! syn_val.has_key? :homology
       next if syn_val[:homology][:assert_cutoff].inject(:+) < 3
       next if ref_annotated.has_key? syn_val[:homology][:hits][0] and ref_annotated[syn_val[:homology][:hits][0]][:partial] == 0
+
       ref_annotated[syn_val[:homology][:hits][0]] = {
         key: prot,
         pId: syn_val[:homology][:pId],
@@ -622,6 +626,7 @@ class BacterialComparator
         length: syn_val[:homology][:length][0],
         partial: (syn_val[:partial] ? 1 : 0)
       }
+
       # ref_annotated[syn_val[:homology][:hits][0]] = {
       #   key: prot,
       #   pId: syn_val[:homology][:pId],
@@ -631,6 +636,7 @@ class BacterialComparator
       #   length: syn_val[:homology][:length][0],
       #   partial: (syn_val[:partial] ? 1 : 0)
       # }
+
     end
 
     # print ref_annotated
@@ -642,13 +648,21 @@ class BacterialComparator
       coverage_query = ""
       query_length = ""
       pId = ""
+
       if ref_annotated[ref_v[:protId]] != nil
-        gene = ref_annotated[ref_v[:protId]][:key]
-        coverage_ref = ref_annotated[ref_v[:protId]][:cov_subject]
-        query_length = query_lengths[ref_annotated[ref_v[:protId]][:key]]
-        coverage_query = ref_annotated[ref_v[:protId]][:cov_query]
-        pId = ref_annotated[ref_v[:protId]][:pId]
-        partial = ref_annotated[ref_v[:protId]][:partial]
+
+        if ref_annotated[ref_v[:protId]][:pId] >= @min_pid and
+          ref_annotated[ref_v[:protId]][:cov_query] >= @min_cov and
+          ref_annotated[ref_v[:protId]][:cov_subject] >= @min_cov
+
+          gene = ref_annotated[ref_v[:protId]][:key]
+          coverage_ref = ref_annotated[ref_v[:protId]][:cov_subject]
+          query_length = query_lengths[ref_annotated[ref_v[:protId]][:key]]
+          coverage_query = ref_annotated[ref_v[:protId]][:cov_query]
+          pId = ref_annotated[ref_v[:protId]][:pId]
+          partial = ref_annotated[ref_v[:protId]][:partial]
+        end
+
       end
 
       _locus_tag = ref_v[:locustag] || ""

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: bacterial-annotator
 version: !ruby/object:Gem::Version
-  version: 0.9.3
+  version: 0.9.4
 platform: ruby
 authors:
 - Maxime Deraspe
 autorequire:
 bindir: bin
 cert_chain: []
-date: 2022-12-19 00:00:00.000000000 Z
+date: 2023-07-17 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bio