[Bioperl-l] strand and frames in SearchIO
Andy Nunberg
anunberg@oriongenomics.com
Mon, 13 May 2002 14:21:19 -0500
I am trying to retrieve strand and frame info for each HSP. I have created
the blast object from a
TBLASTX report using Bio::SearchIO with -format 'psiblast'.
when iterating over the HSPs i call the following method to get strand
information
my ($q_strand,$s_strand)=$hsp->strand();
I can see that sometimes I get subject strand info ($s_strand) and
sometimes I do not.
I also cannot seem to get frame information if I call the frame method on
the hsp object
ie ($q_frame,$s_frame)=$hsp->frame; Although I am not sure which of the
frame methods I am calling when I do this. (I am not clear on exactly
which frame method is inherited)
what is the appropriate method for getting strand and frame info for hsp's
using SearchIO?
Thanks
Andy
*******************************************************************
Andy Nunberg, Ph.D
Computational Biologist
Orion Genomics, LLC
(314) 615-6989
http://www.oriongenomics.com