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