gfa 0.5.0 → 0.6.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 9eaaaf8bcac372e7d6b63d96aee94493faa51d78f5cc2b7edb525c3fbd7efaa2
-  data.tar.gz: 950e9ffa9cd07cdea8ef5cc80b6c162bf009ff3bef5ab757dcda51c443925043
+  metadata.gz: dc98fa1d3479be3fee554dab64876d6cd42795aec08a6b60be6ccb148c57c679
+  data.tar.gz: 28ef77ca42b374e58c983f474d44801060dadb2df912230e431cf40ee8f1386c
 SHA512:
-  metadata.gz: 92e53b9d6c09b1814f3feda86b67dcea6192bd935e3c04dcca7b6ee9c0e7eed919c711941985313c2918f02772f3598883f037638e6e1b7cc9cf8897fec80caf
-  data.tar.gz: 4c14ce2bd91582b7b81bc6123443bd5b0cffe8c83550fb85dc48c8da030d35dafaf393506d537c5df7550b25b87ab10003a938c407bd6210983e5f8841552317
+  metadata.gz: 1a7226dbe5813cc9e7b0cde0e756b0b41395d45e690fd3df0e6ad266f51a209d55d72fdaeaab5aaebd3fec7ce033b54bb04510903cc3eadef230bd30485f037e
+  data.tar.gz: 32f081a9cf35219f05cb21654011a32a1645f13cbf1d7dda320ebef0b7d893523943cd71f8c9f8a310f1021b6adfaa99b1f863d498be1bbb77bf0a5e04779efd

data/bin/gfa-mean-depth

@@ -0,0 +1,37 @@
+#!/usr/bin/env ruby
+
+# @package MiGA
+# @license Artistic-2.0
+
+$LOAD_PATH.push File.expand_path('../../lib', __FILE__)
+$LOAD_PATH.push File.expand_path('../../lib', File.realpath(__FILE__))
+
+require 'gfa'
+
+input, threads = ARGV
+
+unless input
+  $stderr.puts <<~HELP
+    Calculate the average sequencing depth of all segments in the GFA
+    weighted by the segment lengths
+
+    gfa-mean-depth <input> [<threads>]
+
+    <input>      Input GFA file to read
+    <threads>    If passed, parallelize process with these many threads
+  HELP
+  exit(1)
+end
+
+$stderr.puts "Loading GFA: #{input}"
+gfa = GFA.load_parallel(input, (threads || 1).to_i)
+
+$stderr.puts 'Calculating average depth'
+n = gfa.total_length
+avg =
+  gfa.segments.set.map do |segment|
+    raise "Some segments are missing depth data" unless segment.DP
+    segment.DP.value * segment.length / n
+  end.inject(:+)
+puts avg
+

data/bin/gfa-merge

@@ -0,0 +1,34 @@
+#!/usr/bin/env ruby
+
+# @package MiGA
+# @license Artistic-2.0
+
+$LOAD_PATH.push File.expand_path('../../lib', __FILE__)
+$LOAD_PATH.push File.expand_path('../../lib', File.realpath(__FILE__))
+
+require 'gfa'
+
+output = ARGV.shift
+input  = ARGV
+
+if input.empty?
+  $stderr.puts <<~HELP
+    Combine several GFAs into a single GFA. Requires uniqueness of element names
+
+    gfa-merge <output> <input...>
+
+    <output>     Output GFA file to be created
+    <input...>   List of input GFA files to read
+  HELP
+  exit(1)
+end
+
+gfa = GFA.new
+input.each do |i|
+  $stderr.puts "Merging GFA: #{i}"
+  gfa.merge! GFA.load(i)
+end
+
+$stderr.puts "Saving GFA: #{output}"
+gfa.save(output)
+

data/bin/gfa-paths-to-fasta

@@ -0,0 +1,29 @@
+#!/usr/bin/env ruby
+
+# @package MiGA
+# @license Artistic-2.0
+
+$LOAD_PATH.push File.expand_path('../../lib', __FILE__)
+$LOAD_PATH.push File.expand_path('../../lib', File.realpath(__FILE__))
+
+require 'gfa'
+
+input, output, threads = ARGV
+
+unless output
+  $stderr.puts <<~HELP
+    Extract the sequences of the paths from a GFA to FastA file
+
+    gfa-merge <input> <output> [<threads>]
+
+    <input>      Input GFA file to read
+    <output>     Output FastA file to be created
+    <threads>    If passed, parallelize process with these many threads
+  HELP
+  exit(1)
+end
+
+$stderr.puts "Loading GFA: #{input}"
+gfa = GFA.load_parallel(input, (threads || 1).to_i)
+
+

data/lib/gfa/common.rb

@@ -67,4 +67,17 @@ class GFA
   def total_length
     segments.total_length
   end
+
+  ##
+  # Adds the entrie of +gfa+ to itself
+  def merge!(gfa)
+    records.each { |k, v| v.merge!(gfa.records[k]) }
+    self
+  end
+
+  ##
+  # Creates a new GFA based on itself and appends all entries in +gfa+
+  def merge(gfa)
+    GFA.new(opts).merge!(self).merge!(gfa)
+  end
 end

data/lib/gfa/version.rb

@@ -1,5 +1,5 @@
 class GFA
-  VERSION       = '0.5.0'
+  VERSION       = '0.6.0'
   VERSION_ARRAY = VERSION.split(/\./).map { |x| x.to_i } # :nodoc:
   VERSION_MAJOR = VERSION_ARRAY[0] # :nodoc:
   VERSION_MINOR = VERSION_ARRAY[1] # :nodoc:

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: gfa
 version: !ruby/object:Gem::Version
-  version: 0.5.0
+  version: 0.6.0
 platform: ruby
 authors:
 - Luis M. Rodriguez-R
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2023-03-26 00:00:00.000000000 Z
+date: 2023-10-05 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rgl
@@ -65,6 +65,9 @@ files:
 - Rakefile
 - bin/gfa-add-gaf
 - bin/gfa-greedy-modules
+- bin/gfa-mean-depth
+- bin/gfa-merge
+- bin/gfa-paths-to-fasta
 - bin/gfa-subgraph
 - lib/gfa.rb
 - lib/gfa/common.rb