do you want to search for (nearly) exact matches in a chromosom? Then you really should use a suffix array implementation. It is O(m + log(n)) where m is query size and n is sequence size. So probably MUCH faster than your solution. /msg me if you are interested and I send you my mature vmatch module.