[EMBOSS] [BiO BB] Building an alignment from BLAST hsp

Catherine Letondal letondal at pasteur.fr
Tue Mar 28 07:25:07 UTC 2006


On Mar 27, 2006, at 8:03 PM, Ryan Golhar wrote:

> Hi Peter,
>
>> You are quite right that EMBOSS may align the sequences completely
>> differently - unless the HSPs are very significant and cover most
>> of the sequence this will be true of any attempt to simply realign.
>> There has to be some way to pass on the HSPs as fixed positions,
>> as in the BioPerl solution.
>
> I looked at a bioperl method, but can't seem to find something that 
> will
> accomplish this.
>
>> However, it could make a nice EMBOSS application - the only question
>> would be how you would like to specify the HSPs. Perhaps we could read
>
>> BLAST output (in some specified format), or perhaps some other way to
>> give the input alignments.
>
> Yes, I agree.  I suppose the best way would be to specify the two
> sequences and the blast output.  The application could then construct 
> an
> alignment based on a particular HSP (probably the first one, or 
> whatever
> the user specifies).
>

Have you tried this:
http://bioweb.pasteur.fr/seqanal/interfaces/seqsblast.html

It is based on bioperl. check "Get HSP" option (you can even extend it).

Best,

--
Catherine Letondal -- Institut Pasteur -- Computing Center




More information about the EMBOSS mailing list