| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,699,267 | 0 | G | T | 44.0% | 39.7 / 59.9 | 50 | A48A (GCG→GCT) | yrfC | predicted fimbrial assembly protein |
| Reads supporting (aligned to +/- strand): ref base G (5/23); new base T (22/0); total (27/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.87e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCCTCTGCTGCTGGCCGTCGGGATAACGCTAATACTGCGTCTGACAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGC > minE/2699206‑2699331 | cgcCTCTGCTGCTGGCCGTCGGGATAACGCTAATACTGCGTCTGACAGGCAGCGCCGAAGCGCGCATAGa < 1:4204287/70‑1 (MQ=255) tAACGCTAATACTGCGTCTGTCAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGaacaa > 1:706352/1‑71 (MQ=255) tAACGCTAATACTGCGTCTGACAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGaacaa > 1:1054311/1‑71 (MQ=255) tAACGCTAATACTGCGTCTGACAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGaacaa > 1:4307782/1‑71 (MQ=255) tAACGCTAATACTGCGTCTGACAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGaacaa > 1:2971661/1‑71 (MQ=255) tAACGCTAATACTACGTCTGACAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGaacaa > 1:3113567/1‑71 (MQ=255) cTGCGTCTGACAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGGGGaaca < 1:562803/59‑1 (MQ=255) cTGCGTCTGACAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGaaca < 1:2235664/59‑1 (MQ=255) aCAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACa < 1:47102/71‑1 (MQ=255) aCAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACa < 1:4210761/71‑1 (MQ=255) aCAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACa < 1:3807399/71‑1 (MQ=255) aCAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACa < 1:2816409/71‑1 (MQ=255) aCAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACa < 1:2375882/71‑1 (MQ=255) aCAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACa < 1:2141949/71‑1 (MQ=255) aCAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACa < 1:1388096/71‑1 (MQ=255) aCAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACa < 1:1394148/71‑1 (MQ=255) aCAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACa < 1:951839/71‑1 (MQ=255) aCAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACa < 1:1925408/71‑1 (MQ=255) cAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACa < 1:4342660/70‑1 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCa > 1:634646/1‑39 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCTGAACAACAACTCGCCCGCAGCTTa > 1:264474/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACACCTCGCCCGCAGCTTa > 1:1737108/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:1763265/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:1497664/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:4006836/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:3553823/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:349819/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:3454043/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:1762023/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:2121833/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:704751/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:2663485/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:448943/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:2399651/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:1938992/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:1984751/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:2000781/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:2218762/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:924764/1‑67 (MQ=255) aGGCAGCGCCGAAGCTCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:2155764/1‑67 (MQ=255) aGGCAGCGCCGAAGCTAGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTa > 1:3019918/1‑67 (MQ=255) gcCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGaa < 1:2452746/71‑1 (MQ=255) cGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGc < 1:1127707/71‑1 (MQ=255) cGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGc < 1:3158768/71‑1 (MQ=255) cGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGc < 1:3115273/71‑1 (MQ=255) cGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGc < 1:2292139/71‑1 (MQ=255) cGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGc < 1:2195231/71‑1 (MQ=255) cGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGc < 1:1512492/71‑1 (MQ=255) cGAAGCGCGCATAGACGCAGTTTTCCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGc < 1:632900/71‑1 (MQ=255) gAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGc < 1:2787741/70‑1 (MQ=255) | CGCCTCTGCTGCTGGCCGTCGGGATAACGCTAATACTGCGTCTGACAGGCAGCGCCGAAGCGCGCATAGACGCAGTTTTGCTTCAGGCGGAACAACAACTCGCCCGCAGCTTACAGATAACGAAGC > minE/2699206‑2699331 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |