[Bioperl-l] Retrive the sequences from Blast report

Jan Aerts jan.aerts at wur.nl
Thu Jul 29 02:19:15 EDT 2004


Of course, there's the wonderfull HOWTO of SearchIO at
http://bioperl.org/HOWTOs/html/SearchIO.html

jan.

On Thu, 2004-07-29 at 02:46, Sean Davis wrote:
> Jian,
> 
> If you have something like (from Bio::Search::Hit::HitI documentation):
> 
>     my $searchio = new Bio::SearchIO(-format => 'blast', -file =>
> 'result.bls');
>     my $result = $searchio->next_result;
>     my $hit    = $result->next_hit;
> 
>     $hit_name = $hit->name();
> 
>     $desc = $hit->description();
> 
> Depending on where you are getting your blast results, you can probably grab
> the name, description, locus, or accession and use Bio::DB::Genbank to look
> up the sequence.  See Bio::DB::GenBank for details.
> 
> Sean
> ----- Original Message -----
> From: "Sun, Jian" <jsun at utdallas.edu>
> To: <bioperl-l at bioperl.org>
> Sent: Wednesday, July 28, 2004 5:16 PM
> Subject: [Bioperl-l] Retrive the sequences from Blast report
> 
> 
> > Dear Bioperl:
> >     Does anybody could provide me some advice on how to get or retrive the
> sequence from a Blast search report, which will be saved in a sequence file
> and then used for Clustalw alignment?
> >
> > Any help will be appreciated
> > Jian
> >
> >
> > _______________________________________________
> > Bioperl-l mailing list
> > Bioperl-l at portal.open-bio.org
> > http://portal.open-bio.org/mailman/listinfo/bioperl-l
> >
> 
> 
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l


More information about the Bioperl-l mailing list