miga-base 1.3.22.0 → 1.3.22.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 18af6df383be9cb5a5ecba8ed914ec5c1878a70c8e8cb3825a590049b44a06b5
-  data.tar.gz: a79ed870803ce5e4c7eb8e3380c79e43157de8090790d7d0ca9ace51fc6b12c7
+  metadata.gz: 5e4cafb2b821ae172aa9537d5c71c1cafddd6ed73adfa993da1267876dddaedf
+  data.tar.gz: 9fe748cc53152b70fe03a3d85bcf56b517fdae772c7eb1a96329394253f5ca78
 SHA512:
-  metadata.gz: bce0394b7a0f750df581b767992d5eb54f842b40aa06bbfa5dbd8b7d0d6e8f8ae813930432756369e6da07a9b0a95d1c6067147d5e6ae2fae7df697f46a0b41d
-  data.tar.gz: 8f5c085855514811ae1b781d4e58ce3af1b82034d0d92e10b3a558fd2130d7312fba7eb548852e67083eff310435bbb44bde074c342971acca140f9be80c88af
+  metadata.gz: 64c9508befc2bc9c7bf35041ab8ab2dc90aca74e38d80fc401dd6135e9a703135456a8cde5b88e61918007b6249b9e923b279ade15ade5c63ad2799d27b0350f
+  data.tar.gz: 3a74e3c63ca8332c16d5d3583890f6f4e47d8b0afd4cd178046877164083a73d918b37c04a8a8086c6c31d1b55efdf5bc913aa890869220ae5f2e68f8a976a1d

data/lib/miga/project/base.rb

@@ -154,7 +154,14 @@ module MiGA::Project::Base
     },
     aai_save_rbm: {
       desc: 'Should RBMs be saved for OGS analysis?',
-      default: proc { |project| project.clade? },
+      default: proc do |project|
+        project.clade? && project.option(:aai_p) != 'no'
+      end,
+      in: [true, false]
+    },
+    run_ogs: {
+      desc: 'Should orthologous groups be estimated?',
+      default: proc { |project| project.option(:aai_save_rbm) },
       in: [true, false]
     },
     ogs_identity: {

data/lib/miga/project/result.rb

@@ -31,9 +31,11 @@ module MiGA::Project::Result
   ##
   # Is this +task+ to be bypassed?
   def ignore_task?(task)
-    return true if metadata["run_#{task}"] == false
+    opt = "run_#{task}"
+    return true if metadata[opt] == false
+    return true if option?(opt) && option(opt) == false
 
-    !clade? && @@INCLADE_TASKS.include?(task) && metadata["run_#{task}"] != true
+    !clade? && @@INCLADE_TASKS.include?(task) && metadata[opt] != true
   end
 
   ##
@@ -72,20 +74,21 @@ module MiGA::Project::Result
   ##
   # Add result type +:clade_finding+ at +base+ (no +_opts+ supported).
   def add_result_clade_finding(base, _opts)
+    r = nil
     if result_files_exist?(base, %w[.empty])
       r = MiGA::Result.new("#{base}.json")
       r.add_file(:empty, 'miga-project.empty')
-      return r
-    end
-    return nil unless result_files_exist?(base, %w[.proposed-clades])
-    unless clade? ||
-           result_files_exist?(
-             base, %w[.pdf .classif .medoids .class.tsv .class.nwk]
-           )
-      return nil
+    else
+      return nil unless result_files_exist?(base, %w[.proposed-clades])
+      unless clade? ||
+             result_files_exist?(
+               base, %w[.pdf .classif .medoids .class.tsv .class.nwk]
+             )
+        return nil
+      end
+      r = add_result_iter_clades(base)
     end
 
-    r = add_result_iter_clades(base)
     r.add_file(:aai_dist_rds, 'miga-project.dist.rds')
     r.add_file(:aai_dist_rda, 'miga-project.dist.rda')
     r.add_file(:aai_tree,     'miga-project.aai.nwk')

data/lib/miga/version.rb

@@ -12,7 +12,7 @@ module MiGA
   # - String indicating release status:
   #   - rc* release candidate, not released as gem
   #   - [0-9]+ stable release, released as gem
-  VERSION = [1.3, 22, 0].freeze
+  VERSION = [1.3, 22, 2].freeze
 
   ##
   # Nickname for the current major.minor version.
@@ -20,7 +20,7 @@ module MiGA
 
   ##
   # Date of the current gem relese.
-  VERSION_DATE = Date.new(2025, 2, 17)
+  VERSION_DATE = Date.new(2025, 2, 20)
 
   ##
   # References of MiGA

data/utils/distance/commands.rb

@@ -115,37 +115,37 @@ module MiGA::DistanceRunner::Commands
     f1 = tmp_file('largecontigs.fa')
     return unless File.size?(f1)
 
-    # Run FastANI
-    empty = true
-    File.open(f2 = tmp_file, 'w') do |fh|
-      targets.each do |target|
-        target_asm = target&.result(:assembly)&.file_path(:largecontigs)
-        if target_asm
-          fh.puts target_asm
-          empty = false
-        end
+    # Select targets with assemblies
+    target_asm =
+      targets.map do |target|
+        target&.result(:assembly)&.file_path(:largecontigs)
+      end.compact
+
+    # Run FastANI in batches of up to 100 genomes
+    until target_asm.empty?
+      File.open(f2 = tmp_file, 'w') do |fh|
+        target_asm.shift(100).each { |i| fh.puts i }
       end
-    end
-    return if empty
-    run_cmd <<~CMD
-              fastANI -q "#{f1}" --rl "#{f2}" -t #{opts[:thr]} \
-              -o "#{f3 = tmp_file}"
-            CMD
+      run_cmd <<~CMD
+                fastANI -q "#{f1}" --rl "#{f2}" -t #{opts[:thr]} \
+                -o "#{f3 = tmp_file}"
+              CMD
 
-    # Retrieve resulting data and save to DB
-    data = {}
-    File.open(f3, 'r') do |fh|
-      fh.each do |ln|
-        row = ln.chomp.split("\t")
-        n2 = File.basename(row[1], '.gz')
-        n2 = File.basename(n2, '.LargeContigs.fna')
-        data[n2] = [row[2].to_f, 0.0, row[3].to_i, row[4].to_i]
+      # Retrieve resulting data and save to DB
+      data = {}
+      File.open(f3, 'r') do |fh|
+        fh.each do |ln|
+          row = ln.chomp.split("\t")
+          n2 = File.basename(row[1], '.gz')
+          n2 = File.basename(n2, '.LargeContigs.fna')
+          data[n2] = [row[2].to_f, 0.0, row[3].to_i, row[4].to_i]
+        end
       end
-    end
-    batch_data_to_db(:ani, data)
+      batch_data_to_db(:ani, data)
 
-    # Cleanup
-    [f2, f3].each { |i| File.unlink(i) }
+      # Cleanup
+      [f2, f3].each { |i| File.unlink(i) }
+    end
   end
 
   ##

data/utils/subclades.R

@@ -226,8 +226,9 @@ generate_empty_files <- function (out_base) {
   plot(1, t = "n", axes = F)
   legend("center", "No data", bty = "n")
   dev.off()
-  file.create(paste(out_base, ".1.classif", sep = ""))
-  file.create(paste(out_base, ".1.medoids", sep = ""))
+  file.create(paste0(out_base, ".1.classif"))
+  file.create(paste0(out_base, ".1.medoids"))
+  file.create(paste0(out_base, ".empty"))
 }
 
 write_text_report <- function (out_base, ani.d, ani.medoids, ani.types) {

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: miga-base
 version: !ruby/object:Gem::Version
-  version: 1.3.22.0
+  version: 1.3.22.2
 platform: ruby
 authors:
 - Luis M. Rodriguez-R
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2025-02-17 00:00:00.000000000 Z
+date: 2025-02-20 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: daemons