[Bioperl-l] parse Blast Output and Composition Based Statistics
parameter
Hubert Prielinger
hubert.prielinger at gmx.at
Thu Jan 12 18:57:16 EST 2006
Hello,
I want to know, if there is a possibility to get from a Blast Outputfile
the whole Sequence of a protein not only the best local alignment...
for example:
>ref|XP_480077.1| hypothetical protein [Oryza sativa (japonica
cultivar-group)]
dbj|BAD33542.1| hypothetical protein [Oryza sativa (japonica
cultivar-group)]
Length=95
Score = 24.1 bits (47), Expect = 493
Identities = 6/7 (85%), Positives = 7/7 (100%), Gaps = 0/7 (0%)
Query 2 KKRRRWW 8
K+RRRWW
Sbjct 87 KRRRRWW 93
and now, if I parse the file, I want to get the whole Sequence of this
hypothetical protein....is that possible with hsp for example, or any
other way....
my second question is:
I do my blast search with bioperl and the remoteblast module.....each
parameter is working very well, except the composition based statistics
parameter....
it looks like that:
my $factory =
$Bio::Tools::Run::RemoteBlast::HEADER{'COMPOSITION_BASED_STATISTICS'} =
'yes';
it should work like that, but it doesn't....
Thanks for your help in advance......
regards
Hubert
More information about the Bioperl-l
mailing list