Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 328,791 | 0 | A | G | 48.9% | 44.6 / 23.5 | 47 | S182P (TCC→CCC) | betI | choline‑inducible betIBA‑betT divergent operon transcriptional repressor |
Reads supporting (aligned to +/- strand): ref base A (4/20); new base G (23/0); total (27/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.27e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCG > NC_000913/328647‑328941 | cGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAgc < 1:373191/151‑1 (MQ=255) gTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGcgcg < 2:694186/151‑1 (MQ=255) aTGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGGCAGGGGAT‑TAGCGCgggg > 1:676382/1‑149 (MQ=255) ggTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGACGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGttt < 1:822207/150‑1 (MQ=255) tctcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGTAGGGGAATAAAGCGGGGGtttt > 2:338566/1‑151 (MQ=255) tctcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGGAT‑TAGCGCgggggtttt > 1:250942/1‑145 (MQ=255) cGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATcc < 1:889287/150‑1 (MQ=255) gAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGGAA‑TAGCGggggggtttttc > 2:792812/1‑141 (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) | CGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCAGGGAAT‑TAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCG > NC_000913/328647‑328941 |
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: ‑ |