miga-base 1.0.2.0 → 1.0.5.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 42180f18b6cf955c87a7bd4e5d798bba23f529e5a0571409186e6be566067907
-  data.tar.gz: 85d1f429054d5b773938018bd91523dbf0d5f040e266e8bb3bf263eea48b08d7
+  metadata.gz: 9867b3df06290fb7eb25e10b39dd439c9db52d33c9beb3657e76b1f8da570c47
+  data.tar.gz: a53428a16e8d34b48e0e889c54169f287ad34473ad935885b27d7eb74bec090e
 SHA512:
-  metadata.gz: 0d238933adb3c6799b66c5d5e36854c49dd24cc86a330a45c10531fe341a70b5e684a9f6a35e899f582ed64d33494e4300fa4c37236696187b0e2df7da195ad2
-  data.tar.gz: 31ab6d81294dee95801cdc3c38c515493e65d9b727659f856af474212e5ef5dd44d0c733f92d432e1b5d946bbf45580b70a6a815d4163c894ba967bfef343a0c
+  metadata.gz: 53305fda5e76e31d2fae77480de9c6c66b32b3780b27b5f028fba2f95bcebe6bde23635dc6f456244533c5859d732c48c1e595a9b3ddba4a4b8a6cca679ab6df
+  data.tar.gz: 822cd0ded5cb9508406e313abca2fdfb36899dd444e20bbc8839d28fafea2d647b4eeafc87e16e30e6918b0d3bf106ed0a46cfa4482ed6c53caa93f21a98c496

data/lib/miga/cli/action/daemon.rb

@@ -74,6 +74,11 @@ class MiGA::Cli::Action::Daemon < MiGA::Cli::Action
 
   def perform
     cli.operation or raise 'Please specify a daemon operation'
+
+    # Cleanup environment
+    %w[PROJECT RUNTYPE CORES DATASET].each { |i| ENV[i] = nil }
+
+    # Configure and run daemon
     p = cli.load_project
     d = MiGA::Daemon.new(p, cli[:json])
     dopts = %i[latency maxjobs nodelist ppn shutdown_when_done]

data/lib/miga/cli/action/env.rb

@@ -0,0 +1,26 @@
+# @package MiGA
+# @license Artistic-2.0
+
+require 'miga/cli/action'
+
+class MiGA::Cli::Action::Env < MiGA::Cli::Action
+  def parse_cli
+    cli.parse { |_| }
+  end
+
+  def perform
+    puts <<~BASH
+      MIGA="#{MiGA::MiGA.root_path}"
+      MIGA_HOME=${MIGA_HOME:-"$HOME"}
+      . "$MIGA_HOME/.miga_rc"
+      # Ensure MiGA & submodules are first in PATH
+      export PATH="$MIGA/bin:$PATH"
+      for util in enveomics/Scripts FastAAI/FastAAI multitrim ; do
+        export PATH="$MIGA/utils/$util:$PATH"
+      done
+    BASH
+  end
+
+  def empty_action
+  end
+end

data/lib/miga/cli/action/run.rb

@@ -34,7 +34,7 @@ class MiGA::Cli::Action::Run < MiGA::Cli::Action
     if cli[:env]
       cli[:project] ||= ENV['PROJECT']
       cli[:dataset] ||= ENV['DATASET']
-      cli[:thr] ||= ENV['CORES'].to_i unless ENV['CORES'].nil?
+      cli[:thr] = ENV['CORES'].to_i unless ENV['CORES'].nil?
       cli[:result] = File.basename(cli[:result].to_s, '.bash').to_sym
     end
     %i[project dataset result].each do |i|

data/lib/miga/cli/action/wf.rb

@@ -183,6 +183,7 @@ module MiGA::Cli::Action::Wf
     cmd += ['--json', cli[:daemon_json]] unless cli[:daemon_json].nil?
     cmd += ['--max-jobs', cli[:jobs]] unless cli[:jobs].nil?
     cmd += ['--ppn', cli[:threads]] unless cli[:threads].nil?
+    cmd += ['--debug', MiGA::MiGA.debug_trace? ? '2' : '1'] if MiGA::MiGA.debug?
     cwd = Dir.pwd
     call_cli(cmd)
     Dir.chdir(cwd)

data/lib/miga/cli/base.rb

@@ -41,6 +41,7 @@ module MiGA::Cli::Base
     lair: 'Control groups of daemons for several MiGA projects',
     date: 'Return the current date in standard MiGA format',
     console: 'Open an IRB console with MiGA',
+    env: 'Shell code to load MiGA environment',
     # Taxonomy
     tax_set: 'Register taxonomic information for datasets',
     tax_test: 'Return test of taxonomic distributions for query datasets',

data/lib/miga/common/base.rb

@@ -29,10 +29,18 @@ class MiGA::MiGA
     ##
     # Send debug message
     def DEBUG(*args)
-      $stderr.puts(*args) if @@DEBUG
+      $stderr.puts(*args) if debug?
       $stderr.puts(
         caller.map { |v| v.gsub(/^/, '     ') }.join("\n")
-      ) if @@DEBUG_TRACE
+      ) if debug_trace?
+    end
+
+    def debug?
+      @@DEBUG ||= false
+    end
+
+    def debug_trace?
+      @@DEBUG_TRACE ||= false
     end
   end
 end

data/lib/miga/daemon.rb

@@ -329,6 +329,7 @@ class MiGA::Daemon < MiGA::MiGA
   def launch_job(job, hostk = nil)
     # Execute job
     job[:cmd] = job_cmd(job)
+    MiGA::MiGA.DEBUG "CMD: #{job[:cmd]}"
     case runopts(:type)
     when 'ssh'
       # Remote job

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.0, 2, 0].freeze
+  VERSION = [1.0, 5, 1].freeze
 
   ##
   # Nickname for the current major.minor version.
@@ -20,7 +20,7 @@ module MiGA
 
   ##
   # Date of the current gem release.
-  VERSION_DATE = Date.new(2021, 4, 13)
+  VERSION_DATE = Date.new(2021, 8, 26)
 
   ##
   # References of MiGA

data/scripts/miga.bash

@@ -1,21 +1,21 @@
 #!/bin/bash
 
+###
 # Setup environment
 set -e
-MIGA_HOME=${MIGA_HOME:-"$HOME"}
+eval "$("$MIGA/bin/miga" env)"
 SCRIPT=${SCRIPT:-$(basename "$0" .bash)}
-# shellcheck source=/dev/null
-. "$MIGA_HOME/.miga_rc"
-
-# Ensure MiGA & submodules are first in PATH
-export PATH="$MIGA/bin:$PATH"
-for util in enveomics/Scripts FastAAI/FastAAI multitrim ; do
-  export PATH="$MIGA/utils/$util:$PATH"
-done
 
+###
 # Ancillary functions
+
+# Evaluates if the first passed argument is an existing file
 function exists { [[ -e "$1" ]] ; }
+
+# Evaluates if the first passed argument is a function
 function fx_exists { [[ $(type -t "$1") == "function" ]] ; }
+
+# Initiate a project-wide run
 function miga_start_project_step {
   local dir="$1"
   local dir_r="${dir}.running"
@@ -24,6 +24,8 @@ function miga_start_project_step {
   cd "$dir_r"
   miga date > "miga-project.start"
 }
+
+# Finalize a project-wide run
 function miga_end_project_step {
   local dir="$1"
   local dir_r="${dir}.running"

data/utils/enveomics/Manifest/Tasks/blasttab.json

@@ -589,7 +589,17 @@
         "|",
         { "arg": "task" },
         "--blast",
-        "/dev/stdin",
+        "-",
+        {
+          "name": "Output",
+          "opt": "--out",
+          "arg": "out_file",
+          "mandatory": true,
+          "description": [
+            "Output filtered BLAST file.",
+            "Supports compression with .gz extension."
+          ]
+        },
         {
           "opt": "--top",
           "arg": "integer",
@@ -599,19 +609,13 @@
         {
           "opt": "--sort-by",
           "arg": "select",
-          "values": ["bitscore", "evalue", "identity", "length"],
+          "values": ["bitscore", "evalue", "identity", "length", "no"],
           "default": "bitscore",
           "description": "Parameter used to detect the 'best' hits."
         },
         {
           "opt": "--quiet",
           "description": "Run quietly."
-        },
-        ">",
-        {
-          "arg": "out_file",
-          "mandatory": true,
-          "description": "Output (filtered) Tabular BLAST."
         }
       ]
     },

data/utils/enveomics/Manifest/Tasks/sequence-identity.json

@@ -536,8 +536,10 @@
     },
     {
       "task": "rbm.rb",
-      "description": ["Finds the reciprocal best matches between two sets of",
-        "sequences."],
+      "description": [
+        "Finds the reciprocal best matches between two sets of",
+        "sequences."
+      ],
       "help_arg": "--help",
       "cite":[
         ["Camacho et al, 2009, BMC Bioinf (BLAST+)",
@@ -564,11 +566,22 @@
           "mandatory": true,
           "description": "FastA file containing the genome 2."
         },
+        {
+          "name": "Output",
+          "opt": "--out",
+          "arg": "out_file",
+          "mandatory": true,
+          "description": [
+            "Reciprocal Best Matches in BLAST tabular format.",
+            "Supports compression with .gz extension."
+          ]
+        },
         {
           "name": "Nucleotides",
           "opt": "--nucl",
-          "description": ["Sequences are assumed to be nucleotides (proteins",
-            "by default)."]
+          "description": [
+            "Sequences are assumed to be nucleotides (proteins by default)."
+          ]
         },
         {
           "name": "Length",
@@ -582,8 +595,10 @@
           "opt": "--fract",
           "arg": "float",
           "default": 0.0,
-          "description": ["Minimum alignment length (as a fraction of the",
-            "query). If set, requires BLAST+ or Diamond (see Program)."]
+          "description": [
+            "Minimum alignment length (as a fraction of the query).",
+            "If set, requires BLAST+ or Diamond (see Program)."
+          ]
         },
         {
           "name": "Identity",
@@ -603,8 +618,9 @@
           "name": "Executables",
           "opt": "--bin",
           "arg": "in_dir",
-          "description": ["Directory containing the binaries of the search",
-            "program."]
+          "description": [
+            "Directory containing the binaries of the search program."
+          ]
         },
         {
           "opt": "--program",
@@ -612,9 +628,11 @@
           "values": ["blast+", "blast", "blat", "diamond"],
           "default": "blast+",
           "description": "Search program to be used. Default: blast+.",
-          "note": ["Make sure that you have installed the search program you",
+          "note": [
+            "Make sure that you have installed the search program you",
             "want to use. If you have downloaded the program, but it's not",
-            "installed, please use the Executables option above."]
+            "installed, please use the Executables option above."
+          ]
         },
         {
           "opt": "--threads",
@@ -625,12 +643,6 @@
         {
           "opt": "--quiet",
           "description": "Run quietly (no STDERR output)."
-        },
-        ">",
-        {
-          "arg": "out_file",
-          "mandatory": true,
-          "description": "Reciprocal Best Matches in BLAST tabular format."
         }
       ]
     }

data/utils/enveomics/Manifest/examples.json

@@ -26,10 +26,14 @@
     },
     {
       "task": "rbm.rb",
-      "description": ["Reciprocal Best Matches between the proteomes of the",
-        "two major HIV types (HIV-1 and HIV-2)."],
-      "values": ["hiv1.faa","hiv2.faa",null,null,null,null,null,null,null,null,
-        null,null,"hiv1-hiv2.rbm"]
+      "description": [
+        "Reciprocal Best Matches between the proteomes of the",
+        "two major HIV types (HIV-1 and HIV-2)."
+      ],
+      "values": [
+        "hiv1.faa", "hiv2.faa", "hiv1-hiv2.rbm", null, null, null, null, null,
+        null, null,null, null
+      ]
     },
     {
       "task": "ogs.mcl.rb",
@@ -83,10 +87,14 @@
     {
       "_": "== Examples of BLAST statistics and manipulation",
       "task": "BlastTab.topHits_sorted.rb",
-      "description": ["Extract the best match of metagenome-derived proteins",
-        "(from the 'A metagenome') against a Gene Ontology collection."],
-      "values": ["sort","a_mg.cds-go.blast.tsv",null,null,null,null,1,null,null,
-        null,"a_mg.cds-go.blast-bm.tsv"]
+      "description": [
+        "Extract the best match of metagenome-derived proteins",
+        "(from the 'A metagenome') against a Gene Ontology collection."
+      ],
+      "values": [
+        "sort","a_mg.cds-go.blast.tsv", null, null, null, null,
+        "a_mg.cds-go.blast-bm.tsv", 1, null, null
+      ]
     },
     {
       "task": "BlastTab.sumPerHit.pl",

data/utils/enveomics/Scripts/BlastTab.topHits_sorted.rb

@@ -1,101 +1,123 @@
 #!/usr/bin/env ruby
 
-#
-# @author:  Luis M. Rodriguez-R <lmrodriguezr at gmail dot com>
-# @update:  Dec-11-2015
-# @license: artistic license 2.0
-#
+# frozen_string_literal: true
 
-require "optparse"
+$:.push File.expand_path('../lib', __FILE__)
+require 'enveomics_rb/enveomics'
+require 'enveomics_rb/match'
+$VERSION = 1.0
 
-$opts = {n:5, sortby:"bitscore", q:false}
-$cols = {"bitscore"=>11, "evalue"=>10, "identity"=>2, "length"=>3}
-ARGV << "-h" if ARGV.size==0
+o = { n: 5, sortby: :bitscore, out: '-' }
 OptionParser.new do |opts|
-   opts.banner = "Reports the top-N best hits of a BLAST, pre-sorted by query."
-   opts.separator ""
-   opts.separator "Mandatory"
-   opts.on("-i", "--blast FILE",
-      "Path to the BLAST file."){ |v| $opts[:blast]=v }
-   opts.separator ""
-   opts.separator "Optional"
-   opts.on("-n", "--top INTEGER",
-      "Maximum number of hits to report for each query.",
-      "By default: #{$opts[:n]}"){ |v| $opts[:n]=v.to_i }
-   opts.on("-s", "--sort-by STRING",
-      "Parameter used to detect the 'best' hits.",
-      "Any of: bitscore (default), evalue, identity, length."
-      ){ |v| $opts[:sortby]=v }
-   opts.on("-q", "--quiet", "Run quietly."){ $opts[:q]=true }
-   opts.on("-h", "--help", "Display this screen") do
-      puts opts
-      exit
-   end
-   opts.separator ""
+  opts.version = $VERSION
+  Enveomics.opt_banner(
+    opts, 'Reports the top-N best hits of a BLAST, pre-sorted by query',
+    "#{File.basename($0)} -i in.tsv -o out.tsv [options]"
+  )
+  
+  opts.separator 'Mandatory'
+  opts.on(
+    '-i', '--blast FILE',
+    'Path to the BLAST file',
+    'Supports compression with .gz extension, use - for STDIN'
+  ) { |v| o[:in] = v }
+  opts.on(
+    '-o', '--out FILE',
+    'Output filtered BLAST file',
+    'Supports compression with .gz extension, use - for STDOUT (default)'
+  ) { |v| o[:out] = v }
+  opts.separator ''
+  opts.separator 'Filter Options'
+  opts.on(
+    '-n', '--top INTEGER', Integer,
+    'Maximum number of hits to report for each query',
+    "By default: #{o[:n]}"
+  ) { |v| o[:n] = v }
+  opts.on(
+    '-s', '--sort-by STRING',
+    'Parameter used to detect the "best" hits',
+    'Any of: bitscore (default), evalue, identity, length, no (pick first)'
+  ) { |v| o[:sortby] = v.to_sym }
+  opts.separator ''
+  opts.separator 'Other Options'
+  opts.on('-q', '--quiet', 'Run quietly (no STDERR output)') { $QUIET = true }
+  opts.on('-h', '--help', 'Display this screen') { puts opts; exit }
+  opts.separator ''
 end.parse!
 
-abort "-i/--blast is mandatory." if $opts[:blast].nil?
-abort "Unrecognized value for -s/--sortby: #{$opts[:sortby]}." if
-   $cols[ $opts[:sortby] ].nil?
-
-class Hit
-   attr_reader :blast_line
-   def initialize(blast_line)
-      @blast_line = blast_line.chomp.split(/\t/)
-   end
-   def col(i)
-      @blast_line[i]
-   end
-   def <=>(other)
-      ans = self.col( $cols[ $opts[:sortby] ] ).to_f <=> other.col( $cols[ $opts[:sortby] ] ).to_f
-      ans = ans * -1 unless $opts[:sortby] == "evalue"
-      return ans
-   end
-   def to_s
-      @blast_line.join("\t")
-   end
+raise Enveomics::OptionError.new('-i is mandatory') if o[:in].nil?
+unless o[:sortby] == :no || Enveomics::Match.column(o[:sortby])
+  raise Enveomics::OptionError.new("Unrecognized value for -s: #{o[:sortby]}")
 end
 
-class HitSet
-   attr_reader :query, :hits
-   def initialize
-      @hits = []
-      @query = nil
-   end
-   def <<(hit)
-      @query = hit.col(0) if @query.nil?
-      raise "Inconsistent query, expecting #{self.query}" unless
-	 self.query == hit.col(0)
-      @hits << hit
-   end
-   def empty?
-      self.hits.length == 0
-   end
-   def filter!
-      @hits.sort!
-      @hits.slice!($opts[:n], @hits.length)
-   end
-   def to_s
-      @hits.join("\n")
-   end
+class Enveomics::Match
+  attr_accessor :sortby
+
+  def <=>(other)
+    ans = send(sortby) <=> other.send(sortby)
+    sortby == :evalue ? ans : ans * -1
+  end
 end
 
-$stderr.puts "Parsing BLAST." unless $opts[:q]
-fh = File.open $opts[:blast], "r"
-hs = HitSet.new
-while ln=fh.gets
-   hit = Hit.new( ln )
-   if hs.query != hit.col(0)
-      hs.filter!
-      puts hs unless hs.empty?
-      hs = HitSet.new
-      $stderr.print "Parsing line #{$.}... \r" unless $opts[:q]
-   end
-   hs << hit
+class Enveomics::MatchSet
+  attr_reader :query, :hits, :sortby
+
+  def initialize(sortby)
+    @hits = []
+    @query = nil
+    @sortby = sortby
+  end
+
+  def <<(hit)
+    @query ||= hit.qry
+    unless query == hit.qry
+      raise "Inconsistent query, expecting #{query}"
+    end
+
+    @hits << hit.tap { |i| i.sortby = sortby }
+  end
+
+  def empty?
+    hits.empty?
+  end
+
+  def filter!(n)
+    @hits.sort! unless sortby == :no
+    @hits.slice!(n, @hits.length)
+  end
+
+  def to_s
+    hits.join("\n")
+  end
 end
-$stderr.print "Parsed #{$.} lines.   \n" unless $opts[:q]
-fh.close
 
-hs.filter!
-puts hs unless hs.empty?
+begin
+  ifh = reader(o[:in])
+  ofh = writer(o[:out])
+
+  say 'Parsing BLAST'
+  hs = Enveomics::MatchSet.new(o[:sortby])
+  lno = 0
+  ifh.each do |ln|
+    lno += 1
+    hit = Enveomics::Match.new(ln)
+    if hs.query != hit.qry
+      hs.filter! o[:n]
+      ofh.puts hs unless hs.empty?
+      hs = Enveomics::MatchSet.new(o[:sortby])
+      say_inline("Parsing line #{lno}... \r")
+    end
+    hs << hit
+  end
+  say("Parsed #{lno} lines    ")
+  ifh.close
+
+  hs.filter! o[:n]
+  ofh.puts hs unless hs.empty?
+  ofh.close
+rescue => err
+  $stderr.puts "Exception: #{err}\n\n"
+  err.backtrace.reverse.each { |l| $stderr.puts "DEBUG: %s\n" % l }
+  err
+end
 

data/utils/enveomics/Scripts/aai.rb

@@ -316,6 +316,7 @@ Dir.mktmpdir do |dir|
   end
   res = File.open(o[:res], "w") unless o[:res].nil?
   rbm = File.open(o[:rbm], "w") unless o[:rbm].nil?
+  sqlite_db.execute('BEGIN TRANSACTION') unless o[:sqlite3].nil?
   [1,2].each do |i|
     qry_seen = []
     q = "#{dir}/seq#{i}.fa"
@@ -413,6 +414,7 @@ Dir.mktmpdir do |dir|
       o[:sqlite3].nil?
     puts id2/n2 if o[:auto]
   end
+  sqlite_db.execute('COMMIT') unless o[:sqlite3].nil?
   res.close unless o[:res].nil?
   fo.close unless o[:out].nil?
 end

data/utils/enveomics/Scripts/lib/enveomics_rb/match.rb

@@ -4,44 +4,69 @@ module Enveomics
   # A simple object representing a sequence match from a search engine
   # supporting tabular BLAST output
   class Match
-    attr :row
+    class << self
+      def column_types
+        {
+          qseqid: String,    sseqid: String,    pident: Float,
+          length: Integer,   mismatch: Integer, gapopen: Integer,
+          q_start: Integer,  q_end: Integer,    s_start: Integer,
+          s_end: Integer,    evalue: Float,     bitscore: Float,
+          # Non-standard (but frequently used in Enveomics Collection):
+          qry_len: Integer,  sbj_len: Integer
+        }
+      end
 
-    ##
-    # Initialize Enveomics::Match object from a tabular blast line String +ln+
-    def initialize(ln)
-      @row = ln.chomp.split("\t")
-    end
+      def column_type(sym)
+        column_types[colname(sym)]
+      end
 
-    def qry
-      row[0]
-    end
+      def to_column_type(sym, value)
+        case column_type(sym).to_s
+        when 'String' ; value.to_s
+        when 'Float'  ; value.to_f
+        when 'Integer'; value.to_i
+        end
+      end
 
-    def sbj
-      row[1]
-    end
+      def columns
+        column_types.keys
+      end
 
-    def id
-      @id ||= row[2].to_f
-    end
+      def column(sym)
+        columns.index(colname(sym))
+      end
 
-    def len
-      @len ||= row[3].to_i
-    end
+      def colsynonyms
+        {
+          qry: :qseqid, sbj: :sseqid,
+          id: :pident, len: :length, score: :bitscore
+        }
+      end
 
-    def evalue
-      @evalue ||= row[9].to_f
-    end
+      def colnames
+        columns + colsynonyms.keys
+      end
 
-    def score
-      @score ||= row[10].to_f
+      def colname(sym)
+        s = sym.to_sym
+        column_types[s] ? s : colsynonyms[s]
+      end
     end
 
-    def qry_len
-      @qry_len ||= row[12].to_i
+    ####--- Instance Level ---###
+
+    attr :row
+
+    ##
+    # Initialize Enveomics::Match object from a tabular blast line String +ln+
+    def initialize(ln)
+      @row = ln.chomp.split("\t")
     end
 
-    def sbj_len
-      @sbj_len ||= row[13].to_i
+    colnames.each do |sym|
+      define_method sym do
+        self.class.to_column_type(sym, row[self.class.column(sym)])
+      end
     end
 
     def qry_fract

data/utils/enveomics/Scripts/lib/enveomics_rb/utils.rb

@@ -18,10 +18,11 @@ def use(gems, mandatory = true)
 end
 
 def say(*msg)
-  return if $QUIET ||= false
+  $stderr.puts('[%s] %s' % [Time.now, msg.join('')]) unless $QUIET ||= false
+end
 
-  o = '[%s] %s' % [Time.now, msg.join('')]
-  $stderr.puts(o)
+def say_inline(*msg)
+  $stderr.print('[%s] %s' % [Time.now, msg.join('')]) unless $QUIET ||= false
 end
 
 ##

data/utils/enveomics/Scripts/rbm.rb

@@ -2,18 +2,19 @@
 
 # frozen_string_literal: true
 
-$VERSION = 1.0
+$VERSION = 1.01
 $:.push File.expand_path('../lib', __FILE__)
 require 'enveomics_rb/rbm'
 require 'tmpdir'
 
 bms_dummy = Enveomics::RBM.new('1', '2').bms1
-o = { q: false }
+o = { q: false, out: '-' }
 %i[thr len id fract score bin program nucl].each do |k|
   o[k] = bms_dummy.opt(k)
 end
 
 OptionParser.new do |opts|
+  opts.version = $VERSION
   cmd = File.basename($0)
   opts.banner = <<~BANNER
 
@@ -34,6 +35,11 @@ OptionParser.new do |opts|
     '-2', '--seq2 FILE',
     'Path to the FastA file containing the set 2'
   ) { |v| o[:seq2] = v }
+  opts.on(
+    '-o', '--out FILE',
+    'Reciprocal Best Matches in BLAST tabular format.',
+    'Supports compression with .gz extension, use - for STDOUT (default)'
+  ) { |v| o[:out] = v }
   opts.separator ''
   opts.separator 'Search Options'
   opts.on(
@@ -80,7 +86,7 @@ OptionParser.new do |opts|
   ) { |v| o[:thr] = v }
   opts.separator ''
   opts.separator 'Other Options'
-  opts.on('-q', '--quiet', 'Run quietly (no STDERR output)') { o[:q] = true }
+  opts.on('-q', '--quiet', 'Run quietly (no STDERR output)') { $QUIET = true }
   opts.on('-h', '--help', 'Display this screen') { puts opts ; exit }
   opts.separator ''
 end.parse!
@@ -90,10 +96,12 @@ raise Enveomics::OptionError.new('-2 is mandatory') if o[:seq2].nil?
 raise Enveomics::OptionError.new(
   'Argument -f/--fract requires -p blast+ or -p diamond'
 ) if o[:fract] > 0.0 && !%i[blast+ diamond].include?(o[:program])
-$QUIET = o[:q]
 
 rbm = Enveomics::RBM.new(o[:seq1], o[:seq2], o)
-rbm.each { |bm| puts bm.to_s }
+ofh = writer(o[:out])
+rbm.each { |bm| ofh.puts bm.to_s }
+ofh.close
+
 say('Forward Best Matches: ', rbm.bms1.count)
 say('Reverse Best Matches: ', rbm.bms2.count)
 say('Reciprocal Best Matches: ', rbm.count)

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: miga-base
 version: !ruby/object:Gem::Version
-  version: 1.0.2.0
+  version: 1.0.5.1
 platform: ruby
 authors:
 - Luis M. Rodriguez-R
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2021-04-13 00:00:00.000000000 Z
+date: 2021-08-26 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: daemons
@@ -141,6 +141,7 @@ files:
 - lib/miga/cli/action/doctor.rb
 - lib/miga/cli/action/doctor/base.rb
 - lib/miga/cli/action/edit.rb
+- lib/miga/cli/action/env.rb
 - lib/miga/cli/action/files.rb
 - lib/miga/cli/action/find.rb
 - lib/miga/cli/action/generic.rb