[Bioperl-l] PullHSP bug

Sendu Bala bix at sendu.me.uk
Wed Dec 3 00:25:21 UTC 2008


Chris Fields wrote:
> Sendu,
> 
> I found a silent bug while fixing some issues with SimpleAlign which are 
> causing blast_pull tests to fail (HSP-related).  It appears that the 
> query sequence ID is not being caught

I don't see this as a bug of PullHSPI. The blastn report being parsed in 
that failing test has no query sequence id, so obviously PullHSPI can't 
give an ID to the query sequence it creates.

Why does SimpleAlign now fall over when an added sequence has no id? I'd 
suggest the old behaviour be retained...



More information about the Bioperl-l mailing list