miga-base 1.2.17.0 → 1.2.17.2

data/utils/enveomics/Manifest/examples.json

@@ -0,0 +1,162 @@
+{
+  "_": "Input files and directories are included in the 'Tests' folder.",
+  "examples": [
+    {
+      "_": "== Examples of genome comparisons ==",
+      "task": "ogs.stats.rb",
+      "description": ["Statistics on the groups of orthology in the Primate",
+        "Lentivirus Group, including HIV-1, HIV-2, and SIV."],
+      "values": ["primate_lentivirus.ogs",null,null,null,null,null]
+    },
+    {
+      "task": "ani.rb",
+      "description": ["Average Nucleotide Identity (ANI) between two strains",
+        "of Mycoplasma genitalium (M2288 and M2321)."],
+      "values": ["Mgen_M2288.fna","Mgen_M2321.fna",null,null,null,null,null,
+        null,null,null,null,null,null,null,null,null,null,null,null,null,null,
+        null,null,null]
+    },
+    {
+      "task": "aai.rb",
+      "description": ["Average Amino acid Identity (AAI) between Mycoplasma",
+        "genitalium (Bacteria) and Nanoarchaeum equitans (Archaea)."],
+      "values": ["Mgen_M2288.faa","Nequ_Kin4M.faa",null,null,null,null,null,
+        null,null,null,null,null,null,null,null,null,null,null,null,null,null,
+        null,null,null]
+    },
+    {
+      "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", "hiv1-hiv2.rbm", null, null, null, null, null,
+        null, null,null, null
+      ]
+    },
+    {
+      "task": "ogs.mcl.rb",
+      "description": ["Groups of orthology in the Primate Letivirus Group,",
+         "including HIV-1, HIV-2, and SIV."],
+      "values": ["primate_lentivirus.ogs","primate_lentivirus.rbm",null,null,
+        null,null,null,null,null,null,null,null]
+    },
+    {
+      "task": "Table.df2dist.R",
+      "description": ["Transforms a list of AAI values between Xanthomonas",
+        "oryzae genomes into a distance matrix."],
+      "values": ["Xanthomonas_oryzae.aai.tsv",null,null,null,null,100.0,
+        "Xanthomonas_oryzae.aai-mat.tsv"]
+    },
+    {
+      "_": "== Recruitment plots",
+      "task": "BlastTab.catsbj.pl",
+      "description": ["Prepares recruitment plot files for a comparison",
+        "between a virome containing HIV and the HIV-1 genome."],
+      "values": [null,null,null,null,"hiv1.fna","hiv_mix-hiv1.blast.tsv"]
+    },
+    {
+      "task": "BlastTab.recplot2.R",
+      "description": ["Generates recruitment plots for a comparison",
+        "between a virome containing HIV and the HIV-1 genome."],
+      "values": ["hiv_mix-hiv1.blast.tsv",50,100,null,null,null,null,null,null,
+        null,null,null,"hiv_mix-hiv1.Rdata","hiv_mix-hiv1.pdf",null,null]
+    },
+    {
+      "_": "== Examples of functional annotations ==",
+      "task": "HMM.essential.rb",
+      "description": ["Typical single-copy bacterial genes present in",
+        "Mycoplasma genitalium."],
+      "values": ["Mgen_M2288.faa",null,null,null,null,null,null,true,null,null,
+        null,null,null,null,null,null,null,null,null]
+    },
+    {
+      "task": "HMM.essential.rb",
+      "description": ["Typical single-copy archaeal genes present in",
+        "Nanoarchaeum equitans."],
+      "values": ["Mgen_M2288.faa",null,null,null,null,null,null,null,true,null,
+        null,null,null,null,null,null,null,null,null]
+    },
+    {
+      "task": "Newick.autoprune.R",
+      "description": ["Prune an AlkB tree with 110 tips to get only distant",
+        "representatives (41)."],
+      "values": ["alkB.nwk",0.9,null,null,null,null,null,"alkB-pruned.nwk"]
+    },
+    {
+      "_": "== 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,
+        "a_mg.cds-go.blast-bm.tsv", 1, null, null
+      ]
+    },
+    {
+      "task": "BlastTab.sumPerHit.pl",
+      "description": ["Count the number of reads per gene in a mapping of a",
+        "metagenome to a metagenome-derived genes (from the 'A metagenome')."],
+      "values": [null,null,null,null,null,null,null,"a_mg.reads-cds.blast.tsv",
+        null,"a_mg.reads-cds.counts.tsv"]
+    },
+    {
+      "task": "BlastTab.sumPerHit.pl",
+      "description": ["Estimate the total abundance of Gene Ontology",
+        "annotations in the A metagenome, using metagenome-derived proteins,",
+        "and normalizing by the read counts of each protein."],
+      "values": ["a_mg.reads-cds.counts.tsv",null,null,null,null,true,null,
+        "a_mg.cds-go.blast.tsv",null,"a_mg.go.read-counts.tsv"]
+    },
+    {
+      "_": "== Examples of diversity ==",
+      "task": "Table.barplot.R",
+      "description": ["Barplot with the distribution of bacterial phyla in",
+        "four different sites, with taxa sorted by variance."],
+      "values": ["phyla_counts.tsv","250,100,75,200",null,null,null,null,null,
+        null,true,"var",2,null,null,"phyla_counts.pdf",10,null]
+    },
+    {
+      "task": "Chao1.pl",
+      "description": ["Phylum-richness estimated by the Chao1 index with 95%",
+        "confidence, using the distributions of bacterial phyla in four",
+        "different sites."],
+      "values": ["phyla_counts.tsv",null,1,null,null,true,null,
+         "phyla_chao1.tsv"]
+    },
+    {
+      "task": "AlphaDiversity.pl",
+      "description": ["Phylum-diversity estimated by the indices of Shannon",
+        "(H'), Inverse Simpson (1/Lambda), and true diversity of order 1 (1D),",
+        "using the distributions of bacterial phyla in four different sites."],
+      "values": ["phyla_counts.tsv",null,1,null,null,true,null,true,1,null,
+         "phyla_diversity.tsv"]
+    },
+    {
+      "_": "== Other miscelaneous examples ==",
+      "task": "CharTable.classify.rb",
+      "description": ["Classification of anthrax genomes based on can-SNPs, as",
+        "described in Van Ert 2007 (PLoS ONE 2(5):e461)."],
+      "values": ["anthrax-cansnp-data.tsv","anthrax-cansnp-key.tsv",
+        "anthrax-cansnp-classif.tsv","anthrax-cansnp-classif.nwk",null]
+    },
+    {
+      "task": "TRIBS.test.R",
+      "description": ["Test overclustering of Xanthomonas oryzae genomes",
+        "encoding for PilA using Transformed-space Resampling In Biased Sets",
+        "(TRIBS)."],
+      "values": ["Xanthomonas_oryzae.aai-mat.tsv","Xanthomonas_oryzae-PilA.txt",
+        5000,null,null,null,null,0,"Xanthomonas_oryzae-PilA.tribs.Rdata",100]
+    },
+    {
+      "task": "TRIBS.plot-test.R",
+      "description": ["Show the TRIBS-normalized distances between Xanthomonas",
+        "oryzae genomes (grey) and X. oryzae encoding for PilA (red)."],
+      "values": ["Xanthomonas_oryzae-PilA.tribs.Rdata",null,null,null,null,null,
+        null,null,"Xanthomonas_oryzae-PilA.tribs.pdf",null,null]
+    }
+  ]
+}

data/utils/enveomics/Manifest/tasks.json

@@ -0,0 +1,4 @@
+{
+  "_": "This file loads all the .json files inside 'Manifest/Tasks'.",
+  "_include": "Tasks/*.json"
+}

data/utils/enveomics/README.md

@@ -0,0 +1,42 @@
+# Enveomics Collection
+
+Scripts and reference libraries at [Kostas lab](http://enve-omics.gatech.edu).
+
+## Prerequisites
+
+The enveomics collection as a whole has very modest requirements, essentially a
+*nix system with `bash`, `perl`, `ruby`, and `R`. Some scripts may require
+additional libraries, or even external Software, but you'll be forewarned about
+these requirements in the documentation accompanying each script. If you prefer,
+you can also use the Graphical User Interface (GUI), that comes with additional
+tests to let you know if your system is ready to use any given script.
+
+## Graphical User Interface (GUI)
+
+The enveomics collection now has a graphical user interface! To learn more,
+please visit [enveomics-gui](https://github.com/lmrodriguezr/enveomics-gui).
+
+## License
+
+The files in this repository are licensed under the terms of the
+Artistic License 2.0, except when otherwise noted.
+
+You can find a copy of the license in [LICENSE.txt](LICENSE.txt) or at
+http://www.perlfoundation.org/artistic_license_2_0.
+
+## Documentation
+
+Most scripts in this repository are self-documented.  However,
+more extensive documentation (and some discussion) can be found at the
+[documentation website](http://enve-omics.ce.gatech.edu/enveomics/docs).
+Additional documentation for recruitment plots can be found
+[here](Docs/recplot2.md).
+
+## Citation
+
+If you use any of the utilitites in the Enveomics Collection in your research
+please cite:
+
+> Rodriguez-R LM & Konstantinidis KT (2016). The enveomics collection: a toolbox
+> for specialized analyses of microbial genomes and metagenomes.
+> [PeerJ Preprints 4:e1900v1](https://peerj.com/preprints/1900/).

data/utils/enveomics/Scripts/AAsubs.log2ratio.rb

@@ -0,0 +1,171 @@
+#!/usr/bin/env ruby
+#
+# @author  Luis M. Rodriguez-R
+# @update  Dec-21-2015
+# @license artistic license 2.0
+#
+
+$:.push File.expand_path(File.dirname(__FILE__) + "/lib")
+require "enveomics_rb/enveomics"
+
+o = {permutations: 1000, bootstraps: 1000, overwrite: false}
+OptionParser.new do |opt|
+   opt.banner = "
+   Estimates the log2-ratio of different amino acids in homologous sites using
+   an AAsubs file (see BlastPairwise.AAsubs.pl). It provides the point
+   estimation (.obs file), the bootstrap of the estimation (.boot file) and the
+   null model based on label-permutation (.null file).
+
+   Usage: #{$0} [options]".gsub(/^ +/,"")
+   opt.separator ""
+   opt.separator "Mandatory"
+   opt.on("-i", "--input FILE",
+      "Input file in AAsubs format (see BlastPairwise.AAsubs.pl)."
+      ){ |v| o[:file] = v}
+   opt.separator ""
+   opt.separator "Output files"
+   opt.on("-O", "--obs-file FILE",
+      "Output file with the log2-ratios per amino acid.",
+      "By default, '--input value'.obs."
+      ){ |v| o[:obs] = v }
+   opt.on("-B", "--bootstrap-file FILE",
+      "Output file with the bootstrap results of log2-ratios per amino acid.",
+      "By default, '--input value'.boot."
+      ){ |v| o[:boot] = v }
+   opt.on("-N", "--null-file FILE",
+      "Output file with the permutation results of log2-ratios per amino acid.",
+      "By default, '--input value'.null."
+      ){ |v| o[:null] = v }
+   opt.on("--overwrite",
+      "Overwrite existing files. By default, skip steps if the files already" +
+      " exist."){ |v| o[:overwrite] = v }
+   opt.separator ""
+   opt.separator "Parameters"
+   opt.on("-b", "--bootstraps INT",
+      "Number of bootstraps to run. By default: #{o[:bootstraps]}."
+      ){ |v| o[:bootstraps] = v.to_i }
+   opt.on("-p", "--permutations INT",
+      "Number of permutations to run. By default: #{o[:permutations]}."
+      ){ |v| o[:permutations] = v.to_i }
+   opt.on("-q", "--quiet", "Run quietly (no STDERR output)."){ o[:q] = TRUE }
+   opt.on("-h", "--help", "Display this screen.") do
+      puts opt
+      exit
+   end
+   opt.separator ""
+end.parse!
+
+# Initialize
+abort "--input is mandatory" if o[:file].nil?
+ALPHABET = %w(A C D E F G H I K L M N P Q R S T V W Y X)
+o[:obs] ||= "#{o[:file]}.obs"
+o[:boot] ||= "#{o[:file]}.boot"
+o[:null] ||= "#{o[:file]}.null"
+
+# Functions
+def dist_summary(a,b)
+   ALPHABET.map do |i|
+      Math.log(a[i].reduce(0,:+).to_f/b[i].reduce(0,:+), 10)
+   end
+end
+def empty_sample
+   Hash[ALPHABET.map{|k| [k, []]}]
+end
+
+# Initialize
+$stderr.puts "Initializing." unless o[:q]
+sample_A = empty_sample
+sample_B = empty_sample
+last_label = nil
+prot_index = -1
+
+# Read file
+$stderr.puts "Reading input file." unless o[:q]
+ifh = File.open(o[:file], "r")
+ifh.each do |l|
+   r = l.chomp.split /\t/
+   if r.first != last_label
+      prot_index +=1
+      last_label = r.first
+      ALPHABET.each do |a|
+         sample_A[a][prot_index] = 0
+         sample_B[a][prot_index] = 0
+      end
+   end
+   [1,2].each do |ds|
+      unless %w(- *).include? r[ds]
+	 abort "Unknown amino acid in line #{$.}: '#{r[ds]}'." unless
+	    ALPHABET.include? r[ds]
+	 sample_A[ r[ds] ][ prot_index ] += 1 if ds==1
+	 sample_B[ r[ds] ][ prot_index ] += 1 if ds==2
+      end
+   end
+end
+ifh.close
+$stderr.puts "  > Found #{prot_index+1} proteins." unless o[:q]
+$stderr.puts "  > Saving #{o[:obs]}" unless o[:q]
+sum = dist_summary(sample_A, sample_B)
+File.open(o[:obs], "w") do |fh|
+   fh.puts ["AA", "log10_AB"].join("\t")
+   ALPHABET.each do |i|
+      fh.puts [i, sum.shift].join("\t")
+   end
+end
+
+# Permutations
+if File.size? o[:null] and not o[:overwrite]
+   $stderr.puts "Skipping permutations." unless o[:q]
+else
+   $stderr.puts "Permutating." unless o[:q]
+   permut_sum = []
+   o[:permutations].times do |i|
+      permut_A = empty_sample
+      permut_B = empty_sample
+      (0 .. prot_index).each do |j|
+	 # Copy counts of the protein
+	 ALPHABET.each do |k|
+	    permut_A[k][j] = sample_A[k][j]
+	    permut_B[k][j] = sample_B[k][j]
+	 end
+	 # Swap labels at random
+	 permut_A,permut_B = permut_B,permut_A if rand(2)==1
+      end
+      permut_sum << dist_summary(permut_A, permut_B)
+   end
+   $stderr.puts "  > Performed #{o[:permutations]} permutations." unless o[:q]
+   $stderr.puts "  > Saving #{o[:null]}" unless o[:q]
+   File.open(o[:null], "w") do |fh|
+      fh.puts ALPHABET.join("\t")
+      permut_sum.each{ |s| fh.puts s.join("\t") }
+   end
+end
+
+# Bootstraps
+if File.size? o[:boot] and not o[:overwrite]
+   $stderr.puts "Skipping bootstraps." unless o[:q]
+else
+   $stderr.puts "Bootstrapping." unless o[:q]
+   boot_sum = []
+   o[:bootstraps].times do |i|
+      boot_A = empty_sample
+      boot_B = empty_sample
+      (0 .. prot_index).each do |j|
+	 # Sample randomly with replacement
+	 jr = rand(prot_index+1)
+	 # Copy counts of the protein
+	 ALPHABET.each do |k|
+	    boot_A[k][j] = sample_A[k][jr]
+	    boot_B[k][j] = sample_B[k][jr]
+	 end
+      end
+      boot_sum << dist_summary(boot_A, boot_B)
+   end
+   $stderr.puts "  > Performed #{o[:bootstraps]} bootstraps." unless o[:q]
+   $stderr.puts "  > Saving #{o[:boot]}" unless o[:q]
+   File.open(o[:boot], "w") do |fh|
+      fh.puts ALPHABET.join("\t")
+      boot_sum.each{ |s| fh.puts s.join("\t") }
+   end
+end
+
+$stderr.puts "Done. Yayyy!" unless o[:q]

data/utils/enveomics/Scripts/Aln.cat.rb

@@ -0,0 +1,221 @@
+#!/usr/bin/env ruby
+
+# @author  Luis M. Rodriguez-R
+# @license artistic license 2.0
+
+$VERSION = 1.0
+$:.push File.expand_path('../lib', __FILE__)
+require 'enveomics_rb/enveomics'
+
+o = {
+  q: false, missing: '-', model: 'AUTO', removeinvar: false, undefined: '-.Xx?'
+}
+
+OptionParser.new do |opt|
+  cmd = File.basename($0)
+  opt.banner = <<~BANNER
+
+    [Enveomics Collection: #{cmd} v#{$VERSION}]
+
+    Concatenates several multiple alignments in FastA format into a single
+    multiple alignment.  The IDs of the sequences (or the ID prefixes, if using
+    --ignore-after) must coincide across files.
+  
+    Usage: #{cmd} [options] aln1.fa aln2.fa ... > aln.fa
+
+  BANNER
+  opt.on(
+    '-c', '--coords FILE',
+    'Output file of coordinates in RAxML-compliant format'
+  ) { |v| o[:coords] = v }
+  opt.on(
+    '-i', '--ignore-after STRING',
+    'Remove everything in the IDs after the specified string'
+  ) { |v| o[:ignoreafter] = v }
+  opt.on(
+    '-I', '--remove-invariable', 'Remove invariable sites',
+    'Note: Invariable sites are defined as columns with only one state and',
+    'undefined characters.  Additional ambiguous characters may exist and',
+    'should be declared using --undefined'
+  ) { |v| o[:removeinvar] = v }
+  opt.on(
+    '-u', '--missing-char CHAR',
+    "Character denoting missing data. By default: '#{o[:missing]}'"
+  ) do |v|
+    if v.length != 1
+      abort "-missing-char can only be denoted by single characters: #{v}"
+    end
+    o[:missing] = v
+  end
+  opt.on(
+    '-m', '--model STRING',
+    'Name of the model to use if --coords is used. See RAxML docs;',
+    'supported values in v8+ include:',
+    '~ For DNA alignments:',
+    '  "DNA[F|X]", or "DNA[F|X]/3" (to estimate rates per codon position,',
+    '  particular notation for this script)',
+    '~ General protein alignments:',
+    '  "AUTO" (default in this script), "DAYHOFF" (1978), "DCMUT" (MBE 2005;',
+    '  22(2):193-199), "JTT" (Nat 1992;358:86-89), "VT" (JCompBiol 2000;',
+    '  7(6):761-776), "BLOSUM62" (PNAS 1992;89:10915), and "LG" (MBE 2008;',
+    '  25(7):1307-1320)',
+    '~ Specialized protein alignments:',
+    '  "MTREV" (mitochondrial, JME 1996;42(4):459-468), "WAG" (globular, MBE',
+    '  2001;18(5):691-699), "RTREV" (retrovirus, JME 2002;55(1):65-73),',
+    '  "CPREV" (chloroplast, JME 2000;50(4):348-358), and "MTMAM" (nuclear',
+    '  mammal proteins, JME 1998;46(4):409-418)'
+  ) { |v| o[:model] = v }
+  opt.on(
+    '--undefined STRING',
+    'All characters to be regarded as "undefined". It should include all',
+    'ambiguous and missing data chars.  Ignored unless --remove-invariable',
+    "By default: '#{o[:undefined]}'"
+  ) { |v| o[:undefined] = v }
+  opt.on('-q', '--quiet', 'Run quietly (no STDERR output)') { o[:q] = true }
+  opt.on('-V', '--version', 'Returns version') { puts $VERSION ; exit }
+  opt.on('-h', '--help', 'Display this screen') { puts opt ; exit }
+  opt.separator ''
+end.parse!
+files = ARGV
+abort 'Alignment files are mandatory' if files.nil? || files.empty?
+$QUIET = o[:q]
+
+# Read individual gene alignments and return them as a single hash with genome
+# IDs as keys and arrays of single-line strings as values
+#
+# IDs are trimmed after the first occurrence of +ignoreafter+, if defined
+def read_alignments(files, ignoreafter = nil)
+  aln = {}
+  files.each_with_index do |file, i|
+    key = nil
+    File.open(file, 'r').each do |ln|
+      ln.chomp!
+      if ln =~ /^>(\S+)/
+        key = $1
+        key.sub!(/#{ignoreafter}.*/, '') if ignoreafter
+        aln[key] ||= []
+        aln[key][i] = ''
+      else
+        if key.nil?
+          abort "Invalid FastA file: #{file}: Leading line not a def-line"
+        end
+        ln.gsub!(/\s/, '')
+        aln[key][i] += ln
+      end
+    end
+    abort "Empty alignment file: #{file}" if key.nil?
+  end
+  aln
+end
+
+# Remove invariable sites from the alignment hash +aln+, using +undefined+ as
+# a string including all characters representing undefined positions (e.g., X)
+#
+# Returns number of columns removed
+def remove_invariable(aln, undefined)
+  invs = 0
+  lengths = aln.values.first.map(&:length)
+  undef_chars = undefined.chars
+
+  lengths.each_with_index do |len, i|
+    (0 .. len - 1).each do |pos|
+      chr = nil
+      inv = true
+      aln.each_key do |key|
+        next if aln[key][i].nil?
+        chr = aln[key][i][pos] if chr.nil? || undefined.chars.include?(chr)
+        if chr != aln[key][i][pos] && !undef_chars.include?(aln[key][i][pos])
+          inv = false
+          break
+        end
+      end
+      if inv
+        aln.each_key { |key| aln[key][i][pos] = '!' unless aln[key][i].nil? }
+        lengths[i] -= 1
+        invs += 1
+      end
+    end
+    aln.each_key { |key| aln[key][i].gsub!('!', '') unless aln[key][i].nil? }
+  end
+  invs
+end
+
+# Concatenate the alignments hash +aln+ using the character +missing+ to
+# indicate missing alignments, and send each entry in the concatenated alignment
+# to +blk+ as two variables: key (name) and value (alignment string)
+#
+# Returns an array with the lengths of each individual alignment
+def concatenate(aln, missing, &blk)
+  say 'Concatenating'
+  lengths = aln.values.first.map(&:length)
+  aln.each_key do |key|
+    # Pad missing entries
+    lengths.each_with_index { |len, i| aln[key][i] ||= missing * len }
+
+    # Check length
+    obs_len = aln[key].map(&:length)
+    unless lengths == obs_len
+      abort "Inconsistent lengths in '#{key}'\nexp: #{lengths}\nobs: #{obs_len}"
+    end
+
+    # Pass entry to the block and remove from alignment hash
+    blk[key, aln[key].join('')]
+    aln.delete(key)
+  end
+  lengths
+end
+
+# Save the coordinates in +file+ based on +files+ paths (for the names), and
+# using +lengths+ individual alignment lengths
+#
+# The saved format is RAxML coords, including the +model+ for each alignment
+def save_coords(file, names, lengths, model)
+  File.open(file, 'w') do |fh|
+    s = 0
+    names.each_with_index do |name, i|
+      l = lengths[i]
+      next unless l > 0
+      name += "_#{i}" while names.count(name) > 1
+      if model =~ /(DNA.?)\/3/
+        fh.puts "#{$1}, #{name}codon1 = #{s + 1}-#{s + l}\\3"
+        fh.puts "#{$1}, #{name}codon2 = #{s + 2}-#{s + l}\\3"
+        fh.puts "#{$1}, #{name}codon3 = #{s + 3}-#{s + l}\\3"
+      else
+        fh.puts "#{model}, #{name} = #{s + 1}-#{s + l}"
+      end
+      s += l
+    end
+  end
+end
+
+# ------ MAIN ------
+begin
+  say 'Reading'
+  alignments = read_alignments(files, o[:ignoreafter])
+
+  if o[:removeinvar]
+    say 'Removing invariable sites'
+    inv = remove_invariable(alignments, o[:undefined])
+    say "  Removed #{inv} sites"
+  end
+
+  lengths = concatenate(alignments, o[:missing]) do |name, seq|
+    puts ">#{name}", seq.gsub(/(.{1,60})/, "\\1\n")
+  end
+  say "  #{lengths.inject(:+)} columns"
+
+  unless o[:coords].nil?
+    say 'Generating coordinates'
+    names = files.map do |i|
+      File.basename(i).gsub(/\..*/, '').gsub(/[^A-Za-z0-9_]/, '_')
+    end
+    save_coords(o[:coords], names, lengths, o[:model])
+  end
+
+  $stderr.puts 'Done' unless o[:q] 
+rescue => err
+  $stderr.puts "Exception: #{err}\n\n"
+  err.backtrace.each { |l| $stderr.puts l + "\n" }
+  err
+end
+

data/utils/enveomics/Scripts/Aln.convert.pl

@@ -0,0 +1,35 @@
+#!/usr/bin/env perl
+
+#
+# @author: Luis M. Rodriguez-R
+# @update: Mar-23-2015
+# @license: artistic license 2.0
+#
+
+use Bio::AlignIO;
+
+my($iformat,$oformat) = @ARGV;
+($iformat and $oformat) or die "
+Usage:
+   $0 in-format out-format < in_file > output_file
+
+   in-format	Input file's format.
+   out-format	Output file's format.
+   in_file	Input file.
+   out_file	Output file.
+
+Example:
+   # Re-format example.fa into Stockholm
+   $0 fasta stockholm < example.fa > example.stk
+
+Supported formats are:
+   bl2seq, clustalw, emboss, fasta, maf, mase, mega,
+   meme, metafasta, msf, nexus, pfam, phylip, po,
+   prodom, psi, selex, stockholm, XMFA, arp
+
+";
+
+$in  = Bio::AlignIO->new(-fh => \*STDIN, -format => $iformat);
+$out = Bio::AlignIO->new(-fh => \*STDOUT, -format => $oformat);
+while ( my $aln = $in->next_aln ) { $out->write_aln($aln) }
+