added fix to find_genes

[biopieces.git] / bp_bin / find_genes

diff --git a/bp_bin/find_genes b/bp_bin/find_genes
index f8fc0e8110b862758776ba689529d58f8ceff810..16cba9ba59573c37690b6a5297ce9e160388c67e 100755 (executable)
--- a/bp_bin/find_genes
+++ b/bp_bin/find_genes
@@ -55,6 +55,10 @@ prodigal = Prodigal.new(infile, outfile, options)
 prodigal.run
 
 prodigal.each do |record|
+  record[:DEFINITION].match /seqhdr="([^"]+)/ do |m| 
+    record[:S_ID] = $1
+  end
+
   bp.puts record
 end