[Bioperl-l] another searchIO bug?

Anders Stegmann anst at kvl.dk
Thu Feb 16 05:20:06 EST 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.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-0001.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-0002.obj 


More information about the Bioperl-l mailing list