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
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)
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