[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.7, Mon Oct 8 16:36:33 2007 UTC revision 1.8, Thu Oct 11 14:38:33 2007 UTC
# Line 46  Line 46 
46    
47  foreach my $db (@databases){  foreach my $db (@databases){
48    
49      if($db){ $tempf = "/vol/metagenome-48-hour/Data/tblastx-template-nr";}      if($db eq "nr_2007_0208"){ $tempf = "/vol/metagenome-48-hour/Data/tblastx-template-nr";}
50      else{$tempf = "/vol/metagenome-48-hour/Data/tblastx-template-others";}      else{$tempf = "/vol/metagenome-48-hour/Data/tblastn-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;
# Line 63  Line 63 
63          map {print OUT ">$_\n",$fa->{$_},"\n"} @subset;          map {print OUT ">$_\n",$fa->{$_},"\n"} @subset;
64          close OUT;          close OUT;
65    
66          my $job = `/decypher/cli/bin/dc_template -query sequences_to_blast -template $tempf -targ $db`;          my $job = `/decypher/cli/bin/dc_template -query sequences_to_blast -template $tempf -targ $db -priority 9`;
67          $job =~ s/^OK1\s+//;          $job =~ s/^OK1\s+//;
68          chomp($job);          chomp($job);
69          push(@jobs,$job);          push(@jobs,$job);

Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3