[Bioperl-l] Generic Blast parsing issue

CALLEY_JOHN_N@Lilly.com CALLEY_JOHN_N@Lilly.com
Fri, 10 May 2002 10:56:50 -0500


I'm not sure if this is a bug or a missing feature. I'm a long time user 
of the old Blast.pm module. I recently showed a colleague how to set 
things up and suggested that he use the new generic SearchIO feature 
instead of doing things the old way. Unfortunately he needs the offset of 
the HSP in the query and hit sequences. In the Blast.pm module this is 
accessible via the $hsp->range or $hsp->start methods. In the GenericHSP 
object there are no direct equivalents, but there is the get_aln call. 
Since SimpleAlign uses LocatableSeqs internally I hoped that the start and 
end maintained there, might have what I needed. Unfortunately, they are 
always set to 1 and the length of the sequence. Is this a bug? Or a 
misunderstanding on my part? Do I need to give up on using the generic 
modules?

Thanks,
  John Calley