[Bioperl-l] another searchIO bug? with blast report

Anders Stegmann anst at kvl.dk
Fri Feb 17 09:18:18 UTC 2006



>>>Anders Stegmann <anst at kvl.dk> 02/16/06 11:20 am >>>
Hi!

I am blasting a protein seq (query) against an identical seq with a
deletion of Aa nr 61 (subject).
Then I print out the type of nomatch Aa and its position.
The nomatch for the query seq is Aa G at position 61, which is correct.
The nomatch for the subject seq is V at position 60, which is definitely
not correct!?

Is this a bug?

testblast2.pl is the program to run

Q0045 is the query seq.

Q0045del61 is the subject seq (it has to be formated: formatdb -i
Q0045del61 -p T -o F).

Regards Anders.


-------------- next part --------------
A non-text attachment was scrubbed...
Name: Q0045
Type: application/octet-stream
Size: 873 bytes
Desc: not available
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20060217/74838520/attachment-0012.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Q0045del61
Type: application/octet-stream
Size: 872 bytes
Desc: not available
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20060217/74838520/attachment-0013.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: testblast2.pl
Type: application/octet-stream
Size: 6109 bytes
Desc: not available
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20060217/74838520/attachment-0014.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20060217/74838520/attachment-0004.html>


More information about the Bioperl-l mailing list