[Bio] / FigKernelScripts / similar_to.pl Repository:
ViewVC logotype

Annotation of /FigKernelScripts/similar_to.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (view) (download) (as text)

1 : overbeek 1.1 use FIG;
2 :     my $fig = new FIG;
3 :    
4 :     use Sim;
5 :    
6 :     my $usage = "usage: bbhs [CutOff] < PEG > PEG1-Sc-PEG";
7 :    
8 :     my($cutoff,$pair,$peg);
9 :     $cutoff = shift @ARGV;
10 :     $cutoff = defined($cutoff) ? $cutoff : 1.0e-10;
11 :    
12 :     my $id;
13 :     @pegs = map { chop; $_ =~ /(\S+)$/;
14 :     $id = $1;
15 :     if ($id =~ /^fig\|/)
16 :     {
17 :     $id;
18 :     }
19 :     else
20 :     {
21 :     if ($peg = $fig->by_alias($id))
22 :     {
23 :     $peg;
24 :     }
25 :     else
26 :     {
27 :     ();
28 :     }
29 :     }
30 :     } <STDIN>;
31 :    
32 :     %pegs = map { $_ => 1 } @pegs;
33 :    
34 :     foreach $peg (keys(%pegs))
35 :     {
36 :     foreach $sim ($fig->sims($peg,1000,$cutoff,"fig"))
37 :     {
38 :     push(@hits,[$peg,$sim->psc,$sim->id2]);
39 :     }
40 :     }
41 :    
42 :     foreach $x (sort { $a->[1] <=> $b->[1] } @hits)
43 :     {
44 :     print join("\t",@$x),"\n";
45 :     }

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3