[BiO BB] Building an alignment from BLAST hsp
letondal at pasteur.fr
Tue Mar 28 02:25:07 EST 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
> 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
> alignment based on a particular HSP (probably the first one, or
> the user specifies).
Have you tried this:
It is based on bioperl. check "Get HSP" option (you can even extend it).
Catherine Letondal -- Institut Pasteur -- Computing Center
More information about the BBB