[Bioperl-l] another searchIO bug?
Anders Stegmann
anst at kvl.dk
Thu Feb 16 10:20:06 UTC 2006
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/20060216/5062b2cb/attachment-0012.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/20060216/5062b2cb/attachment-0013.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/20060216/5062b2cb/attachment-0014.obj>
More information about the Bioperl-l
mailing list