[Bioperl-l] about bioperl issue
Hsueh
lovebaby39 at gmail.com
Thu Dec 10 09:22:14 EST 2009
Dear
The following is code.
--------------------------------------------------------------------------------
my at params_rb = ( 'program' => 'blastn',
'database' => 'DB\\RB_GUS\\RB_GUS');
my $factory_rb = Bio::Tools::Run::StandAloneBlast->new(@params_rb);
my $input_rb = Bio::Seq->new(-id =>"test_query",
-seq => $testline2);
my $blast_report_rb = $factory_rb->blastall($input_rb);
while (my $result_rb = $blast_report_rb-> next_result ) {
while (my $hit_rb = $result_rb->next_hit()){
while (my $hsp_rb = $hit_rb->next_hsp()){
print $hit_rb->name,"\nevalue = " , $hsp_rb->evalue , "\t score = " , $hsp_rb->score , "\n" ;
#print " ",$hit->name,"\n";
}
}
}
--------------------------------------------------------------------------------
I know how to get "name", "evalue" and "score", but I don't know how to get the word which is in red color. (or please see attachment.)
------------------------------------------------------------------------------------------------------------------
Query: 147 ctctttactcttaggtttacccgccggggtatcgtggcaaacaaggatagtttaaacaga 206
|||||| |||||||||||||||||| |||| || |||||| |||||||||||| ||
Sbjct: 114 ctcttttctcttaggtttacccgccaatatatcctgtcaaacactgatagtttaaactga 173
------------------------------------------------------------------------------------------------------------------
I will appreciate if you could tell me how to do it.
Thank you.
Reginald Hsueh
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: R20080801-1.seq.txt
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20091210/0431bad7/attachment.txt>
More information about the Bioperl-l
mailing list