inprogress
authorcmzmasek@gmail.com <cmzmasek@gmail.com@ca865154-3058-d1c3-3e42-d8f55a55bdbd>
Wed, 16 Jan 2013 03:08:37 +0000 (03:08 +0000)
committercmzmasek@gmail.com <cmzmasek@gmail.com@ca865154-3058-d1c3-3e42-d8f55a55bdbd>
Wed, 16 Jan 2013 03:08:37 +0000 (03:08 +0000)
forester/ruby/evoruby/exe/select_same_gn.rb

index c27526a..1ecd373 100755 (executable)
@@ -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