added missing files

[biopieces.git] / bp_bin / assemble_seq_velvet

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