Table 1 Fit parameters of the modified Gumbel distribution Eq. (18) using the BLOSUM62 scoring matrix and affine gap costs with α = 10, β = 1 . 104 describes the estimated value of λ2 using the scaling relation Eq. (19). Fit parameters for other scoring systems are provided as supplementary material to this artilce [see additional file 1].
From: Local sequence alignments statistics: deviations from Gumbel statistics in the rare-event tail
L, M | λ | 104 λ2 | K | S0 |
| 104 |
|---|---|---|---|---|---|---|
40 | 0.3272 ± 0.108% | 8.6347 ± 0.412% | 0.1028 ± 0.65% | 15.597 ± 0.0676% | 79.05 | 8.1560 ± 12.485% |
60 | 0.3034 ± 0.086% | 6.2007 ± 0.285% | 0.0751 ± 0.60% | 18.455 ± 0.0645% | 49.40 | 6.1711 ± 12.907% |
80 | 0.2892 ± 0.070% | 4.8781 ± 0.222% | 0.0612 ± 0.53% | 20.644 ± 0.0540% | 21.67 | 5.0458 ± 13.280% |
100 | 0.2747 ± 0.072% | 4.3187 ± 0.330% | 0.0472 ± 0.58% | 22.413 ± 0.0611% | 39.42 | 4.3056 ± 13.627% |
150 | 0.2541 ± 0.083% | 3.2974 ± 0.529% | 0.0303 ± 0.61% | 25.682 ± 0.0422% | 39.46 | 3.2047 ± 14.437% |
200 | 0.2432 ± 0.063% | 2.6343 ± 0.344% | 0.0241 ± 0.52% | 28.257 ± 0.0412% | 10.47 | 2.5806 ± 15.214% |
250 | 0.2359 ± 0.071% | 2.1999 ± 0.454% | 0.0198 ± 0.60% | 30.196 ± 0.0459% | 9.40 | 2.1701 ± 15.984% |
300 | 0.2303 ± 0.061% | 1.9101 ± 0.348% | 0.0174 ± 0.54% | 31.934 ± 0.0408% | 2.00 | 1.8758 ± 16.758% |
350 | 0.2261 ± 0.046% | 1.6404 ± 0.239% | 0.0153 ± 0.41% | 33.334 ± 0.0300% | 1.27 | 1.6525 ± 17.544% |
400 | 0.2224 ± 0.052% | 1.4806 ± 0.266% | 0.0136 ± 0.49% | 34.556 ± 0.0369% | 1.36 | 1.4762 ± 18.347% |
600 | 0.2140 ± 0.062% | 1.0206 ± 0.384% | 0.0106 ± 0.64% | 38.561 ± 0.0472% | 2.15 | 1.0250 ± 21.787% |
800 | 0.2090 ± 0.063% | 0.7660 ± 0.419% | 0.0088 ± 0.67% | 41.320 ± 0.0457% | 1.82 | 0.7691 ± 25.697% |