[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