1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
#!/usr/bin/env perl
use strict;
use warnings;
use lib ($ENV{EUK_MODULES});
use Fasta_reader;
use CdbTools;
my $usage = "usage: dbFrom dbTo\n\n";
my $dbFrom = $ARGV[0] or die $usage;
my $dbTo = $ARGV[1] or die $usage;
my $fasta_reader = new Fasta_reader($dbFrom);
while (my $seq_obj = $fasta_reader->next()) {
my $acc = $seq_obj->get_accession();
my $sequence = uc $seq_obj->get_sequence();
$sequence =~ s/\./-/g;
my $other_seq = uc &cdbyank_linear($acc, $dbTo);
$other_seq =~ s/\./-/g;
if ($sequence eq $other_seq) {
print "$acc\tSAME\n";
}
else {
print "$acc\tDIFFERENT\n";
}
}
exit(0);
|
