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

Annotation of /FigKernelScripts/relevant_pchsG.pl

Parent Directory Parent Directory | Revision Log Revision Log


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

1 : overbeek 1.1
2 :     use FIG;
3 :     my $fig = new FIG;
4 :    
5 :     use DBrtns;
6 :     my $dbh = $fig->db_handle;
7 :    
8 :     $usage = "usage: relevant_pchs Genome";
9 :    
10 :     (
11 :     ($genome = shift @ARGV)
12 :     )
13 :     || die $usage;
14 :    
15 :     if (($res = $dbh->SQL("SELECT r1.peg1,r1.peg2,r2.peg2,r3.peg2, r2.iden, r3.iden, r2.paraN, r3.paraN
16 :     FROM close_pegs r1, condensed_sims r2, condensed_sims r3, close_pegs r4
17 :     WHERE r1.peg1 LIKE 'fig|$genome%' AND
18 :     r1.peg1 = r2.peg1 AND
19 :     r1.peg2 = r3.peg1 AND
20 :     r4.peg1 = r2.peg2 AND
21 :     r4.peg2 = r3.peg2")) &&
22 :     (@$res > 0))
23 :     {
24 :     foreach $x (sort { ($a->[0] cmp $b->[0]) or ($a->[1] cmp $b->[1]) or ($b->[4] <=> $a->[4]) } @$res)
25 :     {
26 :     print join("\t",@$x),"\n";
27 :     }
28 :     }

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3