[Bio] / FortyEightMeta / mg_compute_sims_on_timelogic.pl Repository:
ViewVC logotype

Diff of /FortyEightMeta/mg_compute_sims_on_timelogic.pl

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

revision 1.3, Thu Sep 27 19:48:15 2007 UTC revision 1.5, Wed Oct 3 15:42:05 2007 UTC
# Line 5  Line 5 
5  use FIG_Config;  use FIG_Config;
6  use File::Basename;  use File::Basename;
7  use GenomeMeta;  use GenomeMeta;
 #use lib "/home/mkubal/public_html/FIGdisk/dist/releases/dev/FortyEight";  
 #use lib "/home/mkubal/public_html/FIGdisk/dist/releases/dev/PPO";  
8    
9  use DBMaster;  use DBMaster;
10  use Job48;  use Job48;
11  #use SGE;  
12    my $STAGE = "sims";
13    
14  my $fig = new FIG;  my $fig = new FIG;
15    
# Line 41  Line 40 
40  my $fasta = $meta->get_metadata("preprocess.fasta_file");  my $fasta = $meta->get_metadata("preprocess.fasta_file");
41    
42    
43  my $tempf = "/vol/metagenome-48-hour/Data/tblastx-template";  my $tempf;
 #my $tempf = "/home/mkubal/TimeLogic_Templates/tblastx-template";  
44    
45  my @databases = ("nr_2007_0208","16s","lsu","gg");  my @databases = ("nr_2007_0208","16s","lsu","gg");
46    
47  foreach my $db (@databases){  foreach my $db (@databases){
48    
49        if($db){ $tempf = "/vol/metagenome-48-hour/Data/tblastx-template-nr";}
50        else{$tempf = "/vol/metagenome-48-hour/Data/tblastx-template-others";}
51    
52      my $fa=raelib->read_fasta("$fasta"); # this is a generic method of reading a fasta file into a ref to a hash      my $fa=raelib->read_fasta("$fasta"); # this is a generic method of reading a fasta file into a ref to a hash
53      my @keys=keys %$fa;      my @keys=keys %$fa;
54    
# Line 153  Line 154 
154          }          }
155          close(IN);          close(IN);
156      }      }
157    
158        $job->meta->set_metadata("status.$STAGE", "complete");
159        $job->meta->set_metadata("$STAGE.running", "no");
160    
161  }  }

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

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3