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

Annotation of /FigKernelScripts/make_close.pl

Parent Directory Parent Directory | Revision Log Revision Log


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

1 : overbeek 1.1
2 : overbeek 1.2 # usage: make_close [G1 G2 G3 ...]
3 :    
4 : overbeek 1.1 use FIG;
5 :     my $fig = new FIG;
6 :     $dist = 5000;
7 :    
8 : overbeek 1.2 my @genomes = (@ARGV > 0) ? @ARGV : $fig->genomes;
9 :    
10 :     foreach $genome (@genomes)
11 : overbeek 1.1 {
12 :     print STDERR "processing $genome\n";
13 :     if ($fig->is_prokaryotic($genome))
14 :     {
15 :     foreach $peg ($fig->all_features($genome,"peg"))
16 :     {
17 : olson 1.3 my($peg_gnum, $peg_pnum) = $fig->map_peg_to_ids($peg);
18 : overbeek 1.1 my $loc = $fig->feature_location($peg);
19 :     if ($loc)
20 :     {
21 :     my($contig,$beg,$end) = &FIG::boundaries_of($loc);
22 :     if ($contig && $beg && $end)
23 :     {
24 :     my $min = &FIG::min($beg,$end) - $dist;
25 :     my $max = &FIG::max($beg,$end) + $dist;
26 :     my $feat;
27 :     ($feat,undef,undef) = $fig->genes_in_region($genome,$contig,$min,$max);
28 :     foreach $fid (@$feat)
29 :     {
30 : olson 1.3
31 : overbeek 1.1 if (($fid ne $peg) && (&FIG::ftype($fid) eq "peg"))
32 :     {
33 : olson 1.3 my($fgnum, $fpnum) = $fig->map_peg_to_ids($fid);
34 :     print "$peg_gnum\t$peg_pnum\t$fgnum\t$fpnum\n";
35 : overbeek 1.1 }
36 :     }
37 :     }
38 :     }
39 :     }
40 :     }
41 :     }
42 :    

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3