added S_BEG and S_END keys to digest_seq

[biopieces.git] / bp_bin / digest_seq

diff --git a/bp_bin/digest_seq b/bp_bin/digest_seq
index bf5c730b1398470ab5d7acaefb66919a8ba7a035..9fe81e9b83ad60162abc7711b95a3465182a865e 100755 (executable)
--- a/bp_bin/digest_seq
+++ b/bp_bin/digest_seq
@@ -48,6 +48,14 @@ bp.each_record do |record|
 
     digest.each do |subseq|
       new_record            = subseq.to_bp
+
+      if new_record[:SEQ_NAME] =~ /\[(\d+)-(\d+)\]$/
+        s_beg = $1
+        s_end = $2
+      end
+
+      new_record[:S_BEG]    = s_beg 
+      new_record[:S_END]    = s_end
       new_record[:REC_TYPE] = "DIGEST"
       bp.puts new_record
     end