[nextos]

BLAST estimates are derived from extreme value theory and large deviations, which is a very elegant area of probability and statistics.

That's the key part, I think, being able to estimate how unique each alignment is without having to simulate the null distribution, as it was done before with FASTA.

The index also helps, but the speedup comes mostly from the other part.