| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,564,941 | 0 | C | G | 27.3% | 117.3 / 10.9 | 44 | R29P (CGC→CCC) | ECB_01464 | putative fimbrial usher protein |
| Reads supporting (aligned to +/- strand): ref base C (29/3); new base G (0/12); total (29/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.16e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAGCCCTGCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGTGACAGATGTTG > REL606/1564892‑1564990 | aaGCCCTGCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCcgcg > 1:7201571/1‑51 (MQ=255)aaGCCCTGCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCcgcg > 1:7232679/1‑51 (MQ=255)aaGCCCTGCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCcgcg < 1:135567/51‑1 (MQ=255)aaGCCCTGCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCcgcg > 1:6684854/1‑51 (MQ=255) aGCCCTGCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGa < 1:4698366/51‑1 (MQ=255) aGCCCTGCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGa < 1:4991592/51‑1 (MQ=255) gCCCTGCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAg > 1:6477389/1‑51 (MQ=255) cTGCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGCCa > 1:1038088/1‑51 (MQ=255) cTGCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGCCa > 1:4693455/1‑51 (MQ=255) gCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGCCAtt > 1:7082215/1‑51 (MQ=255) aTTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCt > 1:5495305/1‑51 (MQ=255) aTTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCt > 1:1486642/1‑51 (MQ=255) ttCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCaa > 1:2310995/1‑51 (MQ=255) ttCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCaa > 1:2405258/1‑51 (MQ=255) tCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGTCATTCTGCCAAc < 1:3210150/51‑1 (MQ=255) tCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAAc > 1:2477973/1‑51 (MQ=255) tCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAAc > 1:2930113/1‑51 (MQ=255) cttcttCCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACAt > 1:5274974/1‑51 (MQ=255) cttcttCCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACAt > 1:3350348/1‑51 (MQ=255) cttcCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATaaa > 1:4408/1‑51 (MQ=255) cttcCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATaaa > 1:2091429/1‑51 (MQ=255) ttcCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATAAAt > 1:2553807/1‑51 (MQ=255) ttcCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATAAAt > 1:5322019/1‑51 (MQ=255) cccACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATAAAttt > 1:6263824/1‑51 (MQ=255) cccACAGATCAGGAAAAACATAGCCGGGAGCCATTCTGCCAACATaaa < 1:3581494/48‑1 (MQ=255) cccACAGATCAGGAAAAACATAGCCGGGAGCCATTCTGCCAACATAAAtt < 1:3544217/50‑1 (MQ=255) cccACAGATAAGGAAAAACATAGCCGGGAGCCATTCTGCCAACATAAAttt < 1:293484/51‑1 (MQ=255) ccACAGATCAGGAGAAACATAGCCGGGAGCCATTCTGCCAACATAAATTTg < 1:6606158/51‑1 (MQ=255) cacaGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATAAATTTGc > 1:3566127/1‑51 (MQ=255) acaGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATAAATTTGCg > 1:3713080/1‑51 (MQ=255) caGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATAAATTTGCgg > 1:2946319/1‑51 (MQ=255) aGATCAGGAGAAACATAGCCGGGAGCCATTCTGCCAACATAAATTTGCg < 1:5184310/49‑1 (MQ=255) gATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATAAATTTGCGGaa > 1:4442920/1‑51 (MQ=255) aTCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATAAATTTGCGGAAc > 1:626005/1‑51 (MQ=255) aTCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATAAATTTGCGGAAc > 1:127789/1‑51 (MQ=255) ggAGAAACATAGCCGCGAGCCATTCTGCCAACATAAATTTGCGGAACGgag > 1:3420550/1‑51 (MQ=255) aaaaaCATAGCCGGGAGCCATTCTGCCAACATAAATTTGCGGAACGGagag < 1:3479950/49‑1 (MQ=255) aaaaaCATACCCGGGAGCCATTCTGCCAACATAAAtt < 1:4731321/35‑1 (MQ=255) gaAACATAGCCGCGAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGt > 1:2723489/1‑51 (MQ=255) aaaaCATAGCCGGGAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGt < 1:2177930/50‑1 (MQ=255) aaaCATAGCCGGGAGCCATTCTGCCAACATAAATTTGCGGAAc < 1:7172088/43‑1 (MQ=255) aaaCATAGCCGGGAGCCATTCTGCCAACATAAATTTGCGGAACg < 1:7297165/44‑1 (MQ=255) aaaCATAGCCGGGAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGt < 1:4033466/50‑1 (MQ=255) aaaCATAGCCGCGAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGTg > 1:2941655/1‑51 (MQ=255) tAGCCGCGAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGTGACAGa > 1:7132574/1‑51 (MQ=255) gCCGCGAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGTGACAGATg > 1:1378739/1‑51 (MQ=255) ccGGGAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGTGACAGATGt < 1:5667976/51‑1 (MQ=255) ccGCGAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGTGACAGATGt > 1:6258681/1‑51 (MQ=255) cgcgAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGTGACAGATGtt > 1:2004001/1‑51 (MQ=255) gcgAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGTGACAGATGTTg > 1:2990094/1‑51 (MQ=255) gcgAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGTGACAGATGTTg > 1:2721104/1‑51 (MQ=255) | AAGCCCTGCATTTATTCCTTCTTCCCACAGATCAGGAGAAACATAGCCGCGAGCCATTCTGCCAACATAAATTTGCGGAACGGAGAGTGACAGATGTTG > REL606/1564892‑1564990 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |