From f32f933e55cc44d3b49be95aa636a6706ac1a784 Mon Sep 17 00:00:00 2001
From: martinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Date: Wed, 9 Jan 2013 14:24:28 +0000
Subject: [PATCH] adding kmer to contigs

git-svn-id: http://biopieces.googlecode.com/svn/trunk@2069 74ccb610-7750-0410-82ae-013aeee3265d
---
 bp_bin/assemble_seq_velvet | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/bp_bin/assemble_seq_velvet b/bp_bin/assemble_seq_velvet
index d6106fc..abe6bd5 100755
--- a/bp_bin/assemble_seq_velvet
+++ b/bp_bin/assemble_seq_velvet
@@ -159,8 +159,11 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
 		velvet.run_velvetg(options[:ins_length], options[:type].match("Paired"))
 		file_contigs = velvet.pick_best_assembly
 
+    kmer = file_contigs.match(/_\d+/)
+
 		Fasta.open(file_contigs, mode="r") do |fasta_io|
 		  fasta_io.each do |entry|
+        entry.seq_name << "_kmer#{kmer}"
 			  output.puts entry.to_bp
 	    end
 		end
-- 
2.39.2