my $blast = Bio::EnsEMBL::Analysis::RunnableDB::BlastGenscanPep->
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 $slice = $self->fetch_sequence($self->input_id, $self->db);
$self->query($slice);
my %blast = %{$self->BLAST_PARAMS};
my $logic_name = $BLAST_AB_INITIO_LOGICNAME;
$logic_name = 'Genscan' if(!$logic_name);
my $pta = $self->db->get_PredictionTranscriptAdaptor;
my $pts = $pta->fetch_all_by_Slice($self->query, $logic_name);
my $parser = $self->make_parser;
my $filter;
if($self->BLAST_FILTER){
$filter = $self->make_filter;
}
foreach my $t(@$pts){
my $runnable = Bio::EnsEMBL::Analysis::Runnable::BlastTranscriptPep->
new(
-transcript => $t,
-query => $self->query,
-program => $self->analysis->program_file,
-parser => $parser,
-filter => $filter,
-database => $self->analysis->db_file,
-analysis => $self->analysis,
%blast,
);
$self->runnable($runnable);
} } |