1 #!/usr/local/bin/ruby -w
3 require 'lib/evo/sequence/sequence'
4 require 'lib/evo/msa/msa'
5 require 'lib/evo/msa/msa_factory'
6 require 'lib/evo/io/writer/fasta_writer'
7 require 'lib/evo/io/parser/fasta_parser'
20 msa1 = f.create_msa_from_file( input1, FastaParser.new() )
21 msa2 = f.create_msa_from_file( input2, FastaParser.new() )
23 puts "error: " + e.to_s
31 for i in 0 ... msa1.get_number_of_seqs()
32 seq = msa1.get_sequence( i )
35 in_both.add_sequence(seq)
37 only_in_1.add_sequence(seq)
41 for j in 0 ... msa2.get_number_of_seqs()
42 seq = msa2.get_sequence( j )
44 unless msa1.has?( name )
45 only_in_2.add_sequence(seq)
49 puts "only in 1: " + only_in_1.get_number_of_seqs.to_s
50 puts "only in 2: " + only_in_2.get_number_of_seqs.to_s
51 puts "in both : " + in_both.get_number_of_seqs.to_s
54 w.write(only_in_1, "only_in_1.fasta")
55 w.write(only_in_2, "only_in_2.fasta")
56 w.write(in_both, "in_both.fasta")