cgi-bin/tools/sixframe_translate.pl

   1 use CatalystX::GlobalContext qw( $c );
   2 # instead of linking to this, consider just embedding the translations
   3 # directly into the relevant detail page
   4
   5 use strict;
   6 use warnings;
   7
   8 use CGI;
   9 use Bio::PrimarySeq;
  10
  11 my $cgi = CGI->new;
  12 my $seq = get_legacy_est( $cgi ) || get_legacy_unigene( $cgi ) || get_direct_seq( $cgi )
  13     or $c->throw( message => 'No sequence found.', is_error => 0 );
  14
  15 $c->forward_to_mason_view( '/tools/sixframe_translate_standalone.mas',
  16                            seq => $seq,
  17                            blast_url => '/tools/blast'
  18                          );
  19
  20
  21 # ========= helper functions =======
  22
  23 sub get_legacy_est {
  24     my ($cgi) = @_;
  25
  26     my $est_id = $cgi->param('est_id')
  27         or return;
  28
  29     my $est = $c->dbic_schema('SGN::Schema')
  30                 ->resultset('Est')
  31                 ->find( $est_id )
  32         or return;
  33
  34     return Bio::PrimarySeq->new( -id   => "SGN-E".($est_id+0),
  35                                  -seq  => $est->hqi_seq,
  36                                );
  37 }
  38
  39 sub get_legacy_unigene {
  40     my ($cgi) = @_;
  41
  42     my $unigene_id = $cgi->param('unigene_id')
  43         or return;
  44
  45     my $unigene = $c->dbic_schema('SGN::Schema')
  46                     ->resultset('Unigene')
  47                     ->find( $unigene_id )
  48         or return;
  49
  50     return Bio::PrimarySeq->new( -id   => 'SGN-U'.( $unigene_id+0 ),
  51                                  -seq  => $unigene->seq,
  52                                );
  53 }
  54
  55 sub get_direct_seq {
  56     my ($cgi) = @_;
  57
  58     my $seq = $cgi->param('seq')
  59         or return;
  60
  61     my ( $seqid, $sequence ) = split '#', $seq;
  62     return Bio::PrimarySeq->new( -id => $seqid, -seq => $seq );
  63
  64 }
  65