From 2795611c1d3ebe07e0d4b70cecc3f4fdafafdfe5 Mon Sep 17 00:00:00 2001 From: "cmzmasek@gmail.com" Date: Wed, 16 Jan 2013 03:08:37 +0000 Subject: [PATCH] inprogress --- forester/ruby/evoruby/exe/select_same_gn.rb | 8 ++++++++ 1 file changed, 8 insertions(+) diff --git a/forester/ruby/evoruby/exe/select_same_gn.rb b/forester/ruby/evoruby/exe/select_same_gn.rb index c27526a..1ecd373 100755 --- a/forester/ruby/evoruby/exe/select_same_gn.rb +++ b/forester/ruby/evoruby/exe/select_same_gn.rb @@ -24,6 +24,7 @@ module Evoruby exit end + all_names = Set.new gn_to_seqs = Hash.new unique_genes_msa = Msa.new longest_non_unique_genes_msa = Msa.new @@ -36,6 +37,13 @@ module Evoruby for i in 0 ... msa.get_number_of_seqs() seq = msa.get_sequence( i ) name = seq.get_name + if all_names.has?( name ) + puts "error: " + name + " is not unique" + else + all_names.put( name ) + end + + if fragment_re.match( name ) puts "ignored because fragment: " + name frag_counter += 1 -- 1.7.10.2