my $blast = Bio::EnsEMBL::Analysis::RunnableDB::BlastmiRNA->new
(
-analysis => $analysis,
-db => $db,
-input_id => 'contig::AL1347153.1.3517:1:3571:1'
);
$blast->fetch_input;
$blast->run;
my @output =@{$blast->output};
sub fetch_input
{ my ($self) = @_;
my $dna_db = $self->get_dbadaptor($DNA_DBNAME) ;
$self->db->dnadb($dna_db);
my $slice = $self->fetch_sequence($self->input_id, $self->db,'');
$self->query($slice);
my %blast = %{$self->BLAST_PARAMS};
my $parser = $self->make_parser;
my $filter;
if($self->BLAST_FILTER){
$filter = $self->make_filter;
}
my $seq = $self->query;
my $runnable = Bio::EnsEMBL::Analysis::Runnable::BlastmiRNA->new
(
-query => $seq,
-program => $self->analysis->program_file,
-parser => $parser,
-filter => $filter,
-database => $self->analysis->db_file,
-analysis => $self->analysis,
-params =>\% blast,
);
$self->runnable($runnable);
return 1; } |