| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 328,791 | 0 | A | G | 51.4% | 17.9 / 44.6 | 35 | S182P (TCC→CCC) | betI | choline‑inducible betIBA‑betT divergent operon transcriptional repressor |
| Reads supporting (aligned to +/- strand): ref base A (3/14); new base G (18/0); total (21/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.93e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-04 | |||||||||||
GTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGG‑CAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCA > NC_000913/328657‑328937 | gTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTa < 2:465822/151‑1 (MQ=255) tcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGATTAGGGC‑gggggttttt > 2:248329/1‑143 (MQ=255) gAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAg < 1:543177/151‑1 (MQ=255) ggTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGAGGGGGTTAGAGC‑GGGGTTTTTTCAGCgg > 2:658913/1‑151 (MQ=255) gACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGAATAGCGC‑GGGGTTTTTTCAGCGGgtt > 2:33470/1‑149 (MQ=255) gACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGc < 2:341664/151‑1 (MQ=255) cTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAAGGC‑GGGGTTTTTTCAACGGGTTTcc > 2:316390/1‑147 (MQ=255) ggtggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCa < 2:233143/150‑1 (MQ=255) tggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTAGGGC‑GGGGTTTTTTCAGCGGTTTg > 2:130845/1‑141 (MQ=255) gCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAGTCGCGC‑GGGGTTTTTCGCGCCGGGTGCCGCCCag > 2:590653/1‑147 (MQ=255) caccaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGCGCAGGGGAGTTGCGC‑GGGGGTTTTTCAGCGGTTTgttgcctacgcccgc > 1:254501/1‑126 (MQ=255) aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTGAGGGGAATAGGGG‑GGGGGTTTTTCACCGGGTTGCCGCCCCGAGGCCCGCGCCGCCCTaa > 2:411688/1‑148 (MQ=255) aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCCGGGGCATAGCGG‑GGGGGTTATCCAGCCGTGTGCCCCGCCCCGCCGCGCGCAGCCCTAATcc > 1:129848/1‑151 (MQ=255) tGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCGGGGGGGTAACGCGGGGGTTTTTCACACGGTTT‑CCGCCCCGAGCCCCGCgccgcccaaaatc > 2:511118/1‑143 (MQ=255) gTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCa < 2:135475/150‑1 (MQ=255) ataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGAGGGGAATAAAGC‑GGGGGTTTTTCACCGGGTTTGCGCCCCCAGCCGCGCgccgccccaa > 2:651940/1‑138 (MQ=255) tCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAACGC‑GGGGTTTTTTCTGCGGTTTTGCCCTCCGAGCCGCGCGCCGCCCCAATCCATCAAACAACGCGGCCCGGCCgcacc > 1:65583/1‑146 (MQ=255) cGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGccc < 2:631820/151‑1 (MQ=255) ggACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTTATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGc < 1:340859/151‑1 (MQ=255) gACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCCGGGGGGTAACGCGGGGGTTTTTTAACCGGTTTTCC‑CTCCGAGCCGCGCGCCGCCCCAATCCATCAATAAACGCGGGCAGGCCggcaccccc > 2:699214/1‑142 (MQ=255) ctcctcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAGAGC‑GGGGTTTTTTCAACGGGTTGTCGCCCAGAGCCGCGCGCAGCCCCAAAACATCAAACAACGCCGGCCGGCCGCAGCCCCCCTCCCCTGcc > 2:592619/1‑150 (MQ=255) ggTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTGAGGGGATTTGAGC‑GGGGGTTTTTCAGCCGTTTGGCGCCCAGAGCCGCGCCCCGCCCTCAAACATCCCTACACGCGGCCCGGCCGGCGCCCGCGTCCCCTGCCTGTTct > 1:40403/1‑149 (MQ=255) ggTTAATCGGGGGGTAGAGGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGtt < 1:171875/148‑1 (MQ=255) tAATCGGGGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGTTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTcgcgc < 1:56346/150‑1 (MQ=255) tCGGAGGGTAGATGCTGAGTGATAAAGTGGCGGGGGATGGGAATAGGGG‑GGGGGTTTTCTAGCGGGTTGTCGCGCCGAGCCGCGCGCCGCCCCAATCCCTCCATCAACGGGGGCAGGGCCCAGCCCGCGTCCCGCGCCCGTTGGCCCgg > 1:48422/1‑149 (MQ=255) gTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGG‑CAAtt < 2:410322/151‑1 (MQ=255) tAGATGCTGAGTGATAAAGTGGCGGGGCCGGGGAATAACGT‑GGGGGTTTTTCCGCGGGTTGGCGCTCCGAGACCCCCGCCGCCACAAACCATCAACAAACGCGGCCCGGCCGGTACCCGCCTCCTGTGCCTGTTCTCGCGG‑cgactc > 1:457276/1‑142 (MQ=255) aGATGCTGAGTGATAAAGTGGCGGGTCAGGGGAATAACGC‑GGGGTTTTTTCCTCGGGTTTTCGCCCCGAACCGCGCGCCGCCCAAAACCCTCAATCAACGCGGCCCGGCCGTAGCCCCCCTTCTTTGCCTGTTTGCGCGGCCAATTcc > 2:342603/1‑147 (MQ=255) aTGCTGAGTGATAAAGTGGCGGGGCCGGGAAATAGCGC‑GGGGGTTTTTCAGCGGTTTTCCCCCCACAGCCGCGCGCCGCCCCAAACCATCCATCAACGCGGCCCGGCCGGAGCCCCCCTCCTCTGCCTGTTTGCGCCG‑CCATTCGcgac > 2:495473/1‑149 (MQ=255) tGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTGGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGG‑CAATTCGcgacga < 1:140240/150‑1 (MQ=255) gTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGG‑CAATTCGCGACGAAACTCGCTCACCAGa < 2:351936/151‑1 (MQ=255) ggCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGG‑CAATTCGCGACGAAACTCGCTCACCAGAtt < 1:447751/151‑1 (MQ=255) gggTCGGGGAATTAGCGC‑GGTTTTTATCCAGCGTTTTGCCCCTCAGAGCCGCGCCCAGCCATAATCCATCAAGCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGG‑CAATTCGCGACGAAACTCGCTCACCAGATTc < 1:27659/149‑1 (MQ=255) ggTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGG‑CAATTCGCGACGAAACTCGCTCACCAGATTCGAc > 1:545237/1‑151 (MQ=255) ggAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGG‑CAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCa > 2:510630/1‑149 (MQ=255) | GTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGG‑CAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCA > NC_000913/328657‑328937 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |