Those pesky bugs

I have found a bug in the code featured in the previous post where some of the reads hitting several anchors were mapped incorrectly. So to set the record straight, here are the updated mapping rates for the simple code that does not yet handle mismatches or indels:

Dataset	% reads w/ a unique correct mapping	% reads w/ at last one correct mapping
1Mln reads, 100bp	91.07	100.00
300K reads, 300bp	93.32	100.00
10K reads, 1000bp	94.84	100.00

I wish that all bugs were like this and automagically improved performance :) Well, onto the mismatches then...