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

Diff of /FigKernelScripts/check_for_new_genomes.pl

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1, Wed Nov 26 16:21:22 2008 UTC revision 1.2, Tue May 19 01:37:39 2009 UTC
# Line 41  Line 41 
41      &test(@$x,$file_of,$replace,$fig);      &test(@$x,$file_of,$replace,$fig);
42  }  }
43    
44  my %oldH = map { $_ =~ /^(\S+)/; $1 } `cat $old`;  my %oldH = map { $_ =~ /^(\S+)/; $1 => 1 } `cat $old`;
45  my %newH = map { $_ =~ /^(\S+)/; $1 } `cat $new`;  my %newH = map { $_ =~ /^(\S+)/; $1 => 1 } `cat $new`;
46    
47  my($from,$to,$n);  my($from,$to,$n);
48  foreach $_ (keys(%$replace))  foreach $_ (keys(%$replace))
# Line 287  Line 287 
287          }          }
288      }      }
289      my $frac = sprintf("%0.3f",$tot/$all);      my $frac = sprintf("%0.3f",$tot/$all);
290      print STDERR "coverage: $frac\n";  #   print STDERR "coverage: $frac\n";
291      return $frac;      return $frac;
292  }  }
293    
# Line 323  Line 323 
323  #    my $gs1 = $fig->genus_species($genome1);  #    my $gs1 = $fig->genus_species($genome1);
324  #    my $gs2 = $fig->genus_species($genome2);  #    my $gs2 = $fig->genus_species($genome2);
325  #    print STDERR "$genome1 [$gs1] => $genome2 [$gs2]\n";  #    print STDERR "$genome1 [$gs1] => $genome2 [$gs2]\n";
326      $replace->{$genome1} = {$genome2};      $replace->{$genome1} = $genome2;
327  }  }

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3