[Bio] / Sprout / AliasCrunch.pl Repository:
ViewVC logotype

Diff of /Sprout/AliasCrunch.pl

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

revision 1.2, Wed Aug 19 17:07:54 2009 UTC revision 1.3, Thu Mar 25 16:19:32 2010 UTC
# Line 296  Line 296 
296      # Get the parameters.      # Get the parameters.
297      my ($mergeH, $stats) = @_;      my ($mergeH, $stats) = @_;
298      # Loop through the organism directories.      # Loop through the organism directories.
299      for my $orgDir (grep { $_ =~ /^\d+\.\d+$/ } OpenDir($FIG_Config::organisms)) {      for my $orgDir (sort grep { $_ =~ /^\d+\.\d+$/ } OpenDir($FIG_Config::organisms)) {
300          Trace("Processing $orgDir.") if T(3);          Trace("Processing $orgDir.") if T(3);
301          $stats->Add(orgDirGenomes => 1);          $stats->Add(orgDirGenomes => 1);
302          # We need to process all of this organism's TBL files.          # We need to process all of this organism's TBL files.
# Line 323  Line 323 
323                              $stats->Add(orgDirNormal => 1);                              $stats->Add(orgDirNormal => 1);
324                              # Yes. Write it normally.                              # Yes. Write it normally.
325                              WriteToMerge($mergeH, $alias, B => $aliasType, $fid);                              WriteToMerge($mergeH, $alias, B => $aliasType, $fid);
326                            } elsif ($alias =~ /^LocusTag:(.+)/) {
327                                # No, but this is a specially-marked locus tag.
328                                $normalized = $1;
329                                $stats->Add(orgDirLocus => 1);
330                                WriteToMerge($mergeH, $normalized, B => 'LocusTag', $fid);
331                          } elsif ($normalized = AliasAnalysis::IsNatural(LocusTag => $alias)) {                          } elsif ($normalized = AliasAnalysis::IsNatural(LocusTag => $alias)) {
332                              # No, but this is a natural locus tag.                              # No, but this is a natural locus tag.
333                              $stats->Add(orgDirLocus => 1);                              $stats->Add(orgDirLocus => 1);
# Line 491  Line 496 
496                      if ($type eq 'RefSeq' && $alias =~ /^[A-Z][CMT]/) {                      if ($type eq 'RefSeq' && $alias =~ /^[A-Z][CMT]/) {
497                          $stats->Add(correspondingContig => 1);                          $stats->Add(correspondingContig => 1);
498                      } else {                      } else {
499                          # We're okay, so normalize the alias.                          # Check for a locus tag disguised as a CMR ID.
500                          my $normalized = AliasAnalysis::Normalize($type => $alias);                          my $realType = $type;
501                            if ($type eq 'CMR' && $alias =~ /^[A-Z]{2,3}_\d+$/) {
502                                $realType = 'LocusTag';
503                                $stats->Add(correspondingLocus => 1);
504                            } else {
505                          $stats->Add(correspondingNormal => 1);                          $stats->Add(correspondingNormal => 1);
506                            }
507                            # Normalize the alias.
508                            my $normalized = AliasAnalysis::Normalize($type => $alias);
509                          # Write it out once for each FIG ID.                          # Write it out once for each FIG ID.
510                          for my $fid (@fids) {                          for my $fid (@fids) {
511                              WriteToMerge($mergeH, $normalized, A => $type, $fid);                              WriteToMerge($mergeH, $normalized, A => $type, $fid);

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

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3