[Bioperl-l] HSP alignments using AlignIO

Charles Hauser chauser@duke.edu
09 Sep 2002 10:00:47 -0400


I am trying to capture blast alignments to be loaded into a database.  
Using the code below I am getting 'generic' ID's, e.g.query_[BlastHSP] & hit_[BlastHSP].  How do I
have it return the true blast query and blast hit? 

I would also like to save these as <pre>  </pre> html - ideas?
 

my $algn = Bio::AlignIO->new(-file => ">out.aln",
			     -format => 'clustalw',
			     -idlength => 25);


$algn->write_aln($hsp->get_aln());




CLUSTAL W(1.81) multiple sequence alignment


query_[BlastHSP] 7/1-24 SSPRAAPPLGRGWPP-----PGASGCRAPSP-PPAHRAPGTTP---------*AAPAA
hit_[BlastHSP] 7/1-24   STPPDASPAGASSLPEVPAAPPAEDAPAPSPEPDSPPRPETAPAPEATALLGDAAAIA
			*:*  *.* * .  *     * *... **** * :   * *:*          **. *

regards,

Charles