Local sequence alignments statistics: deviations from Gumbel statistics in the rare-event tail

摘要

BackgroundThe optimal score for ungapped local alignments of infinitely long random sequences is known to follow a Gumbel extreme value distribution. Less is known about the important case, where gaps are allowed. For this case, the distribution is only known empirically in the high-probability region, which is biologically less relevant.