[Bio] / SeedViewer / MetagenomeAnalysis.pm Repository:
ViewVC logotype

Diff of /SeedViewer/MetagenomeAnalysis.pm

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

revision 1.4, Tue Apr 22 19:56:18 2008 UTC revision 1.5, Tue Apr 29 14:08:42 2008 UTC
# Line 497  Line 497 
497    
498  =pod  =pod
499    
500    =item * B<get_recruitment_plot_data> (I<genome>)
501    
502    Given a genome id (83333.1), this method returns all sequence ids,
503    the alignment length, the match id and the taxonomy string for all
504    sequences which match the criteria and have their tax_str start equal
505    the genome tax string I<filter>.
506    
507    =cut
508    
509    sub get_recruitment_plot_data {
510      my ($self, $genome) = @_;
511    
512      my $table = $self->dbtable;
513      my $dbid  = $self->get_dataset_id("SEED");
514      my $where = $self->get_where_clause();
515      $where = ($where) ? "and $where" : '';
516    
517      my ($tax_id) = $self->dbh->selectrow_array("select tax_str from rdp_to_tax where seq_num='". $genome . "'");
518    
519      my $sth = $self->dbh->prepare("select s.id1, s.id2, s.b2, s.e2, s.logpsc from ( select id1, min(rank_psc) as rank from $table where dbid=$dbid $where group by id1) as b inner join $table as s on b.id1=s.id1 and b.rank=s.rank_psc where dbid=$dbid and s.tax_str='$tax_id'");
520    
521      $sth->execute;
522      my $result = $sth->fetchall_arrayref();
523    
524      return $result;
525    
526    }
527    
528    
529    
530    
531    =pod
532    
533  =item * B<get_hits_for_sequence> (I<seq_id>, I<dataset>, I<limit>)  =item * B<get_hits_for_sequence> (I<seq_id>, I<dataset>, I<limit>)
534    
535  Given a sequence id I<seq_id> (id1) and a dataset name (db_id), this method returns  Given a sequence id I<seq_id> (id1) and a dataset name (db_id), this method returns

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

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3