| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 328,803 | 0 | T | G | 47.7% | 48.0 / 28.8 | 44 | T178P (ACC→CCC) | betI | choline‑inducible betIBA‑betT divergent operon transcriptional repressor |
| Reads supporting (aligned to +/- strand): ref base T (4/19); new base G (21/0); total (25/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.43e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCT > NC_000913/328656‑328950 | ggTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGACGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGttt < 1:822207/150‑1 (MQ=255) tctcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGTAGGGGAATAAAGCGGGGGtttt > 2:338566/1‑151 (MQ=255) cGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATcc < 1:889287/150‑1 (MQ=255) gCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCt < 1:20515/151‑1 (MQ=255) cTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGGAT‑TAACGCGGGGTTTTTTCAGCGGGTTGTCGccc > 2:773072/1‑149 (MQ=255) ggtggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAA‑TAGCGCGGGGTTTTTTCACCGGgtt > 2:738754/1‑142 (MQ=255) aGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCcgcgc < 2:781100/151‑1 (MQ=255) ggTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCcgcgcg < 1:652696/151‑1 (MQ=255) ggTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCCGGGGAT‑TAACGCGGGGTTTTTCCCCCGGGTTGCCGccc > 2:458524/1‑138 (MQ=255) ataACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCGCGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAg < 1:678394/151‑1 (MQ=255) taACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAG‑TAGCGCGGGGGTTTTTCAGCGGGTTGGCGCTCCGAGCCGCGcgccgc > 2:485669/1‑151 (MQ=255) ccaccaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCAt < 1:392202/151‑1 (MQ=255) ccaccaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCCGGGGAA‑TAGCGGGGGGTTTTTTCAACGGGTTGCCCCTCCGAGCCCCGCgccgcc > 2:805849/1‑149 (MQ=255) gTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCat < 2:121977/151‑1 (MQ=255) atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCatca < 1:360066/151‑1 (MQ=255) tataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGGAT‑TAGCGCGGGGTTTTTTCACCGGGTTGCCGCCCCGAGCCGCGCGCCGCCAAAATACa > 2:894779/1‑147 (MQ=255) aaGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGGAA‑TAACGCGGGGTTTTTTCAGCGGGTTGCCGCCCCCAGCCGCGCGCCGCCCTAAACCaacca > 2:746923/1‑143 (MQ=255) gctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAA‑TAGCGCGGGGGTTTTTCCGCGGGTTTTCGCCCAGAGCCGCGCGCAGCCCTAATCCatccaataaccc > 2:523775/1‑140 (MQ=255) ctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGGAT‑TAGCGCGGGGGTTTTTCAGCGGGTTGTCGCCCACAGCCGCGCGCCGGCACAATCCATCAATCAGCGCgg > 2:317035/1‑150 (MQ=255) gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGcc < 1:433581/150‑1 (MQ=255) gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCa < 2:655900/151‑1 (MQ=255) tCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCCGGGGAA‑TTGCGGGGGGTTTTTCCCACGGGTTGCCCCCCCCAGCCGCGCGCCGCCATAATACATCAAACAACGCgg > 2:539194/1‑146 (MQ=255) cATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAA‑TAACGCGGGGTTTTTTCAGCGGGTTGGCGCCCCCAGCCGCGCGCAGCCATAAACCATCAAACAACGcggcccggcc > 2:464457/1‑148 (MQ=255) aCATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCCGGGGAT‑TAGAGCGGGGTTTTTTCAACGGGTTGCCGCCCCGAGCCGCGCGCCGCCCTAAATCATCAATCAACGcggcccggcccgaaccc > 2:615832/1‑140 (MQ=255) cACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCt < 2:413692/150‑1 (MQ=255) gTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGc < 2:901408/151‑1 (MQ=255) cctcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTg < 1:880018/150‑1 (MQ=255) ctcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAT‑TAGCGCGGGGTTTTTTCAGCGGGTTTGCGCCCCGAGCCGCGCGCAGCCATAAACCATCAATCAACGCGGCCCGGCCGCAGCCCGCCTCCCCTGCCTGtt > 2:120266/1‑151 (MQ=255) ggTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAT‑TAGCGCGGGGTTTTTTCAGCGGGTTGTCGCTCAGAGCCGCGCGCCGCCAAAATCCCTCAATCAACGCGGCCAGGGCCGAGCCCGCCTCCTCTGCCTGtt > 2:275045/1‑148 (MQ=255) gTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTcgcg < 2:813517/151‑1 (MQ=255) ttAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGAAGGGAAA‑TAGCGCGGGGTTTTTTCACCGGGTTTTCGCTCCCAGCCCCGCGCCGCCCCAATCCCTAAATCACCGCGGCCACGCCCTAACCCCCCTCCCGTGccccgt > 1:504168/1‑142 (MQ=255) aaTCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAA‑TAGAGCGGGGTTTTTTCACCGGGTTGGCGCCCACAGCCGCGCGCCGCCCTAATCCATCAAACAACGCGGCCCGGCCGCAGCCCCCCTCCTGTGCCTGttttcgc > 2:910750/1‑144 (MQ=255) aTCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCCGGGGAT‑TAACGCGGGGTTTTTTCAGCGGGTTTTCGCCCCGAGCCGCGCGCCGCCCAAATCCATCAATCAGCGCGGCCAGGCCGCAGCCCCCTTCCTCTGCCCGTTTGCGCgg > 2:48679/1‑150 (MQ=255) tGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCCGGGGAT‑TAACGCGGGGTTTTTTCAACGGGTTTTCGCCCCGAGCCGCGCGCCGCCCCAATCCCTCAATCAACGCGGCCCCGCCGCAACCCCCCTCCCGTGCCCGTTTTCGCGGCCAtt > 1:450135/1‑150 (MQ=255) gggTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAT‑TAGCGCGGGGTTTTTTCCGCGGGTTGTCGCCCCGAGCCCCGCGCCGCCATAATCCATCAATCAACGCCGCCAGGCCGTAGCCCGCCTCCTCTGCCCGTTTGCGCGGCaaat > 1:746449/1‑147 (MQ=255) tGAGTGATAAAGTGGC‑GGGTCAGGGGAA‑TAGCGCGGGGTTTTTTCAACGGTTTTCCGCCCCGAGCCGCGCGCAGCCATAATCCATCAAACAACGCGGCCCGGCCGGAGCCCCCCTCCTGTGCCCGTTTGCGCGGCAATTCGCCACGgaa > 2:179079/1‑149 (MQ=255) tGAGTGATAAAGTGGC‑GGGTCAGGGGAA‑TAAAGCGGGGTTTTTTCAGCGGTTTGTCGCCCCGACGCGCGCGGCGGCATAAACCCCTAATCAACGCGGGCCGGCCGCAGGCCCCTCTCCGTGCCTGTTCGCGCGGCCAtt > 1:586362/1‑139 (MQ=255) aTAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTc < 1:824284/150‑1 (MQ=255) aaaGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCAc < 1:408446/150‑1 (MQ=255) tCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACa < 1:623003/150‑1 (MQ=255) aGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGc > 1:732577/1‑150 (MQ=255) gAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAgc > 1:102894/1‑151 (MQ=255) aaT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAgcg > 2:904094/1‑151 (MQ=255) aGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGcggcgg > 2:276243/1‑151 (MQ=255) cgGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGgctg < 2:114572/150‑1 (MQ=255) gggTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGgctgct < 1:396960/151‑1 (MQ=255) | GGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCT > NC_000913/328656‑328950 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |