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

Annotation of /FigKernelScripts/test_figfam_peg.pl

Parent Directory Parent Directory | Revision Log Revision Log


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

1 : arodri7 1.6 use Data::Dumper;
2 : overbeek 1.1 use FIG;
3 :     my $fig = new FIG;
4 :    
5 :     $usage = "usage: test_figfam_peg PEG [FigFamsData]";
6 :    
7 :     (
8 :     ($peg = shift @ARGV)
9 :     )
10 :     || die $usage;
11 :    
12 : olson 1.5 $ff_data = (@ARGV > 0) ? $ARGV[0] : &FIG::get_figfams_data();
13 : overbeek 1.1 use FigFams;
14 :     $figfams = new FigFams($fig,$ff_data);
15 :     $seq = $fig->get_translation($peg);
16 :    
17 :     $ENV{'DEBUG'} = 1;
18 : arodri7 1.6 my ($in_ff) = $figfams->families_containing_peg($peg);
19 :     print "In FIGfams? ";
20 :     if ($in_ff){
21 :     my $famO = new FigFam($fig,$in_ff,$ff_data);
22 :     print "yes, $in_ff, " . $famO->family_function . "\n";
23 :     }
24 :     else{
25 :     print "no\n";
26 :     }
27 : overbeek 1.1
28 : arodri7 1.6 print "Sequence Annotation: " . $fig->function_of($peg) . "\n\n";
29 : overbeek 1.4 ($got,$sims) = $figfams->place_in_family($seq,1);
30 : overbeek 1.1 if (defined($got))
31 :     {
32 :     print &Dumper($sims,$got->family_id,$got->family_function);
33 :     }
34 :     else
35 :     {
36 : overbeek 1.4 print "Sorry, no call, trying loose check\n";
37 : overbeek 1.2 ($got,$sims) = $figfams->place_in_family($seq,1,1);
38 :     if (defined($got))
39 :     {
40 :     print &Dumper($sims,$got->family_id,$got->family_function);
41 :     }
42 :     else
43 :     {
44 : overbeek 1.4 print "Sorry, no call\n";
45 : overbeek 1.2 }
46 : overbeek 1.1 }
47 :    
48 :    

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3