[Bioperl-l] SearchIO, StandAloneBlastPlus

michael watson (IAH-C) michael.watson at bbsrc.ac.uk
Wed Mar 17 20:47:57 UTC 2010


I think that relies on the blast database being built with the "-o T" option, which is not the default for formatdb....
________________________________________
From: bioperl-l-bounces at lists.open-bio.org [bioperl-l-bounces at lists.open-bio.org] On Behalf Of Smithies, Russell [Russell.Smithies at agresearch.co.nz]
Sent: 17 March 2010 19:48
To: 'Janine Arloth'; 'bioperl-l at lists.open-bio.org'
Subject: Re: [Bioperl-l] SearchIO, StandAloneBlastPlus

If you're running blast locally, use fastacmd to extract the sequences from the blast database.
Eg fastacmd -d nr -S AC147927

Russell Smithies

Bioinformatics Applications Developer
T +64 3 489 9085
E  russell.smithies at agresearch.co.nz

Invermay  Research Centre
Puddle Alley,
Mosgiel,
New Zealand
T  +64 3 489 3809
F  +64 3 489 9174
www.agresearch.co.nz



> -----Original Message-----
> From: bioperl-l-bounces at lists.open-bio.org [mailto:bioperl-l-
> bounces at lists.open-bio.org] On Behalf Of Janine Arloth
> Sent: Monday, 15 March 2010 9:16 p.m.
> To: bioperl-l at lists.open-bio.org
> Subject: [Bioperl-l] SearchIO, StandAloneBlastPlus
>
> Hello,
>
> exists a possibility to get/extract the whole hit sequences? (Not only the
> hit string from the alignment with $hsp->$hit_string;)
>
> Best regards
>
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
=======================================================================
Attention: The information contained in this message and/or attachments
from AgResearch Limited is intended only for the persons or entities
to which it is addressed and may contain confidential and/or privileged
material. Any review, retransmission, dissemination or other use of, or
taking of any action in reliance upon, this information by persons or
entities other than the intended recipients is prohibited by AgResearch
Limited. If you have received this message in error, please notify the
sender immediately.
=======================================================================

_______________________________________________
Bioperl-l mailing list
Bioperl-l at lists.open-bio.org
http://lists.open-bio.org/mailman/listinfo/bioperl-l




More information about the Bioperl-l mailing list