X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bp_bin%2Fassemble_seq_velvet;h=e2753ce91ad65023d47627bd1c372553075fb14e;hb=2f0fd91b461033529a4a72e161bd133252a22eb6;hp=abe6bd59f76ea342de5130d6ef783c2d9c53306a;hpb=b1fadc0ce9224b6a13f8deae332bb902a54c65af;p=biopieces.git diff --git a/bp_bin/assemble_seq_velvet b/bp_bin/assemble_seq_velvet index abe6bd5..e2753ce 100755 --- a/bp_bin/assemble_seq_velvet +++ b/bp_bin/assemble_seq_velvet @@ -101,7 +101,7 @@ class Velvet count = 0 n50 = 0 - Fasta.open(file, mode="r") do |fasta_io| + Fasta.open(file, "r") do |fasta_io| fasta_io.each do |entry| total += entry.length lengths << entry.length @@ -144,7 +144,7 @@ Dir.mkdir(options[:directory]) unless Dir.exists?(options[:directory]) file_fasta = File.join(options[:directory], "sequence_in.fna") Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output| - Fasta.open(file_fasta, mode="w") do |fasta_io| + Fasta.open(file_fasta, "w") do |fasta_io| input.each_record do |record| if record[:SEQ_NAME] and record[:SEQ] seq = Seq.new_bp(record) @@ -161,7 +161,7 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output| kmer = file_contigs.match(/_\d+/) - Fasta.open(file_contigs, mode="r") do |fasta_io| + Fasta.open(file_contigs, "r") do |fasta_io| fasta_io.each do |entry| entry.seq_name << "_kmer#{kmer}" output.puts entry.to_bp