adding kmer to contigs

git-svn-id: http://biopieces.googlecode.com/svn/trunk@2069 74ccb610-7750-0410-82ae-013aeee3265d

diff --git a/bp_bin/assemble_seq_velvet b/bp_bin/assemble_seq_velvet
index d6106fc0f3fbe7c996282437d5c4c6d6d15ea628..abe6bd59f76ea342de5130d6ef783c2d9c53306a 100755 (executable)
--- 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