[Bioperl-l] Strand direction in Bl2seq

Julio Fernandez Banet julipao at terra.es
Wed May 28 01:46:20 EDT 2003


Hello.
First of all say thanks for your help and keep the good work as I found this
mailing list very helpful while trying to write my first lines of code.

I'm writing some code using bl2seq.
My problem is when I try to get the direction of the alignment using
$hit_direction=Hsp->strand('hit') it seems it doesn't work and I only get
hit as a result and not the direction in the format 1, -1, 0 as I good like
to get.
I tried writing it in diferent ways:
$hit_direction=Hsp->hit->strand();

I don't know if there any way to use it with bl2seq instead of blast, but if
this is not the way to get the direction of the alignment any help on how to
know if the strands are Plus/Plus or Plus/Minus would be really helpful.

Thanks a lot for your help



More information about the Bioperl-l mailing list