require 'biopieces'
require 'fasta'
-require 'pp'
class Velvet
def initialize(directory, sequence_file, verbose)
def fasta_n50(file)
total = 0
lengths = []
+ count = 0
+ n50 = 0
Fasta.open(file, mode="r") do |fasta_io|
fasta_io.each do |entry|
end
end
- n50(total, lengths)
- end
-
- def n50(total, lengths)
- count = 0
- n50 = 0
-
lengths.sort.reverse.each do |length|
count += length