print_counts( per_species_counter, log, Constants::LINE_DELIMITER )
per_species_counter = 0
end
- puts " " + species_counter.to_str + ": " + current_species + " [" + my_readlink + "]"
- log << species_counter.to_str + ": " + current_species + " [" + my_readlink + "]" + Constants::LINE_DELIMITER
+ puts " " + species_counter.to_s + ":" + current_species + " [" + my_readlink + "]"
+ log << species_counter.to_s + ": " + current_species + " [" + my_readlink + "]" + Constants::LINE_DELIMITER
species_counter += 1
end
puts " " + seq_name
Util.fatal_error( PRG_NAME, "error: seq name \"" + seq_name + "\" not unique" )
end
- # if current_msa.find_by_name_start( seq_name, true ).size > 0
- # begin
- # seq = current_msa.get_by_name_start( seq_name, true ).copy
- # rescue ArgumentError => e
- # Util.fatal_error( PRG_NAME, "error: " + e.to_s )
- # end
- # else
- # # Not found, try finding by partial match.
- # begin
- # seq = current_msa.get_by_name( seq_name, true, true )
- # rescue ArgumentError => e
- # Util.fatal_error( PRG_NAME, "error: " + e.to_s )
- # end
- # end
-
normalized_id = per_species_counter.to_s( 16 ).upcase +
"_" + current_species
rescue Exception => e
Util.fatal_error( PRG_NAME, "error: " + e.to_s )
end
- if @seqs <= 10000000
+ if @seqs <= 600000
@file_to_msa[ input ] = msa
@seqs += msa.get_number_of_seqs
- puts "total seqs in memory: " + @seqs.to_str
+ puts " total seqs in memory: " + @seqs.to_s
end
msa
end