Bio::EnsEMBL::Analysis::Runnable::BlastTranscriptDNA

  my $blast = Bio::EnsEMBL::Analysis::Runnable::BlastTranscriptDNA->
  new(
      -transcript => $transcript,
      -program => 'wublastn',
      -database => 'embl_vertrna',
      -options => 'hitdist=40 -cpus=1',
      -parser => $bplitewrapper,
      -filter => $featurefilter,
     );
  $blast->run;
  my @output =@{$blast->output};

  my ( $self, $features )  = @_;
  
  # for each feature
  my @features = sort{ $a->start <=> $b->start} @$features;
  for my $feature ( @features ) {
    my %exon_hash = ();
    # for each ungapped piece in it
    my @ungapped = $feature->ungapped_features;
    for my $ugFeature ( @ungapped ) {
      my @split = $self->transcript->pep2genomic($ugFeature->start(),
                                              $ugFeature->end());
      
      my $cdna_total = 1;
      foreach my $gcoord ( @split ) {
        if($gcoord->isa('Bio::EnsEMBL::Mapper::Gap')) {
          $cdna_total += $gcoord->end - $gcoord->start + 1;
          next;
        }
        
        my $cdna_start = $cdna_total;
        my $gstart  = $gcoord->start;
        my $gend    = $gcoord->end;
        my $gstrand = $gcoord->strand;
        my $cdna_end = $gend - $gstart + $cdna_start;
        $cdna_total += $gend - $gstart + 1;
        
        #determine which exon this genomic coordinate overlaps
        my $exon;
        foreach my $e (@{$self->transcript->get_all_Exons}) {
          if($gstart >= $e->start && $gend <= $e->end) {
            $exon = $e;
            last;
          }
        }


        # first, eat away non complete codons from start
        while(( $cdna_start - 1 ) % 3 != 0 ) {
          $cdna_start++;
          if( $gstrand == 1 ) {
            $gstart++;
          } else {
            $gend--;
          }
        }
        
        # and from end
        while( $cdna_end  % 3 != 0 ) {
          $cdna_end--;
          if( $gstrand == 1 ) {
            $gend--;
          } else {
            $gstart++;
          }
        }
        
        if( $cdna_end <= $cdna_start ) {
          next;
        }
        #recalculate the hit coordinates.  They may be split up into
        # seperate features by introns
        my ($hstrand, $hstart, $hend);
        $hstrand = $feature->hstrand();
        if($hstrand == 1) {
          $hstart = $cdna_start - 1 + $ugFeature->hstart();
          $hend = $cdna_end - 1 + $ugFeature->hstart();
        } else {
          $hend = $ugFeature->hend() - $cdna_start + 1;
          $hstart = $ugFeature->hend() - $cdna_end + 1;
        }
        my $fp = $self->feature_factory->create_feature_pair
          ($gstart, $gend, $gstrand, $feature->score, $hstart,
           $hend, $hstrand, $feature->hseqname, $feature->percent_id, 
           $feature->p_value);
        #store generated feature pairs, hashed on exons
        push( @{$exon_hash{$exon}}, $fp );
      }
    }
    # Take the pieces for each exon and make gapped feature
    foreach my $ex ( keys %exon_hash ) {
      my $dna_align_feature = Bio::EnsEMBL::DnaDnaAlignFeature->new
        (-features => $exon_hash{$ex});
      $self->output([$dna_align_feature]);
    }
  }

  my ($class,@args) = @_;
  my $self = $class->SUPER::new(@args);
  my ($pt) = rearrange(['TRANSCRIPT'], @args);
  $self->transcript($pt);
  return $self;

  my ($self, $arr_ref, $reset) = @_;

  if(!$self->{'output'}){
    $self->{'output'} = [];
  }
  if($arr_ref){
    throw("Must pass Runnable:output an arrayref not a ".$arr_ref)
      unless(ref($arr_ref) eq 'ARRAY');
    if ($reset) {
      $self->{'output'} = $arr_ref;
    } else {
      push(@{$self->{'output'}}, @$arr_ref);
    }
  }
  return $self->{'output'};

  my ($self, $dir) = @_;

  $self->workdir($dir) if($dir);

  my $pep = $self->transcript->translate;
  $pep->id($self->transcript->dbID);

  if($pep->length <= 3){
    #transcripts this length cause problems for blast
    return;
  }

  my $query = $self->query;
  $self->query($pep);  
  $self->SUPER::run($dir);
  $self->query($query);

  my $out = $self->output;
  $self->output([], 1);
  $self->align_hits_to_query($out);

  my ($self, $pt) = @_;
  if($pt){
    throw("BlastGenscanDNA:transcript must be a ".
          "Bio::EnsEMBL::Transcript ") 
      unless($pt->isa("Bio::EnsEMBL::Transcript"));
    $self->{'prediction_transcript'} = $pt;
  }
  return $self->{'prediction_transcript'};