gfa 0.6.3 → 0.7.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 6e858670d104b7602c27ef71b2942f355027eb3e34dfc9c9aa5d9aa6a3a8ebbc
-  data.tar.gz: 4a905a283972c8d0704c90feb81a90becc73af7dffb328cf108145c0e0103023
+  metadata.gz: a7a646a57c98f25bda2120b3e04783099433d7aa7c8f0658a1bfdcd07e0dea35
+  data.tar.gz: c8aec87b21fd9b1d7972228fa7118cf71b40f58a5dd20e385e6e80ef86b2ef92
 SHA512:
-  metadata.gz: cda7d12b5f2c6bd18e991d8591a2788d85a2dd29c99bfcc0e0d20c40978199e81099875453d9f34a188e4f723177acc09c2990d25062ef7082556a96f97b56ac
-  data.tar.gz: 18e5ff044ffd143e95fbc8525756215437528468b094136f622a30cfde7a51b8d25ca4c1e8a26300a3027ec79fafdbdda343af3cebf9e1c2813467c983b67873
+  metadata.gz: 3f4c54a483e4cd6baaf66fbddda2e7eb41a13cd14bd492c5132a2d2610bdd2f283d7a92c261dca10f54fb9e03611ca7f7c05036fa783c361d4b8e752e9fd6b53
+  data.tar.gz: 2b17c35013f3c9dbdb99a20e060c2dde13732c139f0551f3ce5b65e317ad5d1918b84447570542b11b67f5d65387729606de914b18994c420e5fe8b41b9bf5af

data/bin/gfa-paths-to-fasta

@@ -8,29 +8,42 @@ $LOAD_PATH.push File.expand_path('../../lib', File.realpath(__FILE__))
 
 require 'gfa'
 
-input, output, threads = ARGV
+input, output, threads, orphan = ARGV
 
 unless output
   $stderr.puts <<~HELP
     Extract the sequences of the paths from a GFA to FastA file
 
-    gfa-merge <input> <output> [<threads>]
+    gfa-merge <input> <output> [<threads> [<orphan>]]
 
     <input>      Input GFA file to read
     <output>     Output FastA file to be created
     <threads>    If passed, parallelize process with these many threads
+    <orphan>     If set to 'yes', it also outputs the sequences of all
+                 segments that are not included in any path
   HELP
   exit(1)
 end
 
 $stderr.puts "Loading GFA: #{input}"
+orphan ||= 'no'
 gfa = GFA.load_parallel(input, (threads || 1).to_i)
 
 $stderr.puts "Saving path sequences: #{output}"
 File.open(output, 'w') do |fasta|
+  in_path = {}
   gfa.paths.set.each do |path|
     fasta.puts '>%s' % path.path_name.value
     fasta.puts path.sequence(gfa)
+    path.segment_names_a.each { |i| in_path[i] = true }
+  end
+
+  if orphan == 'yes'
+    gfa.segments.set.each do |segment|
+      next if in_path[segment.name.value]
+      fasta.puts '>S:%s' % segment.name.value
+      fasta.puts segment.sequence.value
+    end
   end
 end
 

data/lib/gfa/version.rb

@@ -1,5 +1,5 @@
 class GFA
-  VERSION       = '0.6.3'
+  VERSION       = '0.7.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.6.3
+  version: 0.7.0
 platform: ruby
 authors:
 - Luis M. Rodriguez-R
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2023-10-06 00:00:00.000000000 Z
+date: 2023-10-10 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rgl