committing major ruby overhaul

[biopieces.git] / bp_bin / join_seq

diff --git a/bp_bin/join_seq b/bp_bin/join_seq
index de80ad2a53246c67385bc4c3a0d28a9f721e87f4..f071f41d0dcc6a5165267a85f426a719f2a21fa8 100755 (executable)
--- a/bp_bin/join_seq
+++ b/bp_bin/join_seq
@@ -34,28 +34,29 @@ require 'maasha/seq'
 
 casts = []
 
-bp = Biopieces.new
-
-options = bp.parse(ARGV, casts)
+options = Biopieces.options_parse(ARGV, casts)
 
 seq = Seq.new(nil, "")
 
-bp.each_record do |record|
-  if record.has_key? :SEQ
-    unless seq.seq_name 
-      seq.seq_name = record[:SEQ_NAME]
+Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
+  input.each_record do |record|
+    if record.has_key? :SEQ
+      unless seq.seq_name 
+        seq.seq_name = record[:SEQ_NAME]
+      end
+      seq.seq << record[:SEQ]
     end
-    seq.seq << record[:SEQ]
+
+    output.puts record
   end
 
-  bp.puts record
+  if seq.seq != ""
+    new_record = seq.to_bp
+    new_record[:REC_TYPE] = "JOIN"
+    output.puts new_record
+  end
 end
 
-if seq.seq != ""
-  new_record = seq.to_bp
-  new_record[:REC_TYPE] = "JOIN"
-  bp.puts new_record
-end
 
 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<