[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