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

Ryan Golhar golharam at umdnj.edu
Mon Mar 27 18:03:39 UTC 2006


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).

Ryan




More information about the EMBOSS mailing list