Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 941,713 | 0 | T | G | 42.2% | 66.3 / 12.5 | 45 | V252G (GTG→GGG) | dmsA | dimethyl sulfoxide reductase, anaerobic, subunit A |
Reads supporting (aligned to +/- strand): ref base T (3/23); new base G (19/0); total (22/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.32e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT > NC_000913/941566‑941862 | actactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCAGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTATGGTAATAACCATGGCGAAACGCGAATCAGTGGCGTTGGGGTg < 1:66381/149‑1 (MQ=255) actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt < 2:489670/149‑1 (MQ=255) tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtttc > 2:429621/1‑148 (MQ=255) gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTGTATc > 2:665724/1‑148 (MQ=255) cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTATTATCTTTaa > 1:474720/1‑150 (MQ=255) aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑g < 1:576237/151‑1 (MQ=255) cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTATCTTTaaaa > 1:817049/1‑143 (MQ=255) cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTATTATCTTTAAAA‑gg > 1:849/1‑147 (MQ=255) aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCa < 1:179279/151‑1 (MQ=255) aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCa < 1:778157/151‑1 (MQ=255) aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCCATCTTTAACA‑AGCACGcc > 2:690582/1‑149 (MQ=255) ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGTCTTAATATCTTTaaaaaggacgac > 1:159346/1‑138 (MQ=255) ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTTCTATCTTTAACA‑AGGACGccc > 2:449147/1‑148 (MQ=255) tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGaa < 1:229612/136‑1 (MQ=255) ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAA‑AGCACGCCAGaaa > 2:462643/1‑150 (MQ=255) tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAAt < 2:268231/150‑1 (MQ=255) tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAAt < 2:24854/150‑1 (MQ=255) gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAATTACTATCTTTAAAA‑AGCACGCCAGaaa > 2:179909/1‑148 (MQ=255) ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACA‑GGCACGCCCGaaaaataa > 1:347898/1‑143 (MQ=255) taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAA‑AGCACGCCCGAAAACTaaaggc > 2:380899/1‑147 (MQ=255) gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGatcatcat < 1:275916/151‑1 (MQ=255) gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGatcatc < 2:49747/149‑1 (MQ=255) gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACA‑GGCACGCCAGAAATCTAAAGGCCGCATGgtcattat > 1:194002/1‑149 (MQ=255) aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAt < 2:254699/151‑1 (MQ=255) aTGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTAATAACTTGAACA‑GGGACGCCCGAAAACGTATGGCCGCATGatgattat > 1:239082/1‑145 (MQ=255) tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTAACTTGAAAA‑GGGACGCCAGAAATCTAATGcccccgcgagca > 2:311451/1‑131 (MQ=255) aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtat < 1:82946/150‑1 (MQ=255) cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc < 1:282703/151‑1 (MQ=255) atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTTAAAA‑GGCACGCCCGAAAACTAATTCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCgg > 1:321555/1‑151 (MQ=255) atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTg < 2:459225/151‑1 (MQ=255) cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAA‑GGCACGCCAGAAATATAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGcc > 1:453712/1‑151 (MQ=255) aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGaa < 2:399555/150‑1 (MQ=255) gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAA‑CGCCCGCCCGAAATTTAAAGCCCGCAAGATCATCATAGATTCGCCCTATACCCACACCGGGGCCGGGCGGGGAGATg > 1:356311/1‑148 (MQ=255) aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTATATCTTTAACA‑AGCACGCCCGAAAACTAAAGCCCGCATGATCATCATTGATCCGCGCTATACCCAAAACGGTGCCGGGCGGGAAGATGAGTgg > 1:657832/1‑151 (MQ=255) cTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAA‑GGCACCCCAGAAAAATAAAGCCCCCATGATCATCAACGATCCCCGCTATACCGACACCCGTGCCGGGGGCGAAGATGAGTGGATcc > 2:588108/1‑151 (MQ=255) ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc < 2:168707/151‑1 (MQ=255) ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCAACGATCCGAGCTATGCCGACTCCGGTGCCGGGCGCGAAGATGAGTGGATCCATa < 1:449919/148‑1 (MQ=255) cccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGa < 1:64774/151‑1 (MQ=255) ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg < 2:462458/151‑1 (MQ=255) gCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGa < 2:347935/151‑1 (MQ=255) cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt < 1:824136/151‑1 (MQ=255) gcgAATGAGTGGCGGTGGGGGGACCTATTATCTTTAAAACGGCCCG‑CAAAAAACTAATGCCCGCAAGAACCTCATCGTTCCCCGCTCTACCCCCACCCGGGGCCGGCGCGAAAATGAGCGGATCCCCCTTCCTCCCGGGAAAGaag > 1:797251/1‑144 (MQ=255) gAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGtt < 1:519978/151‑1 (MQ=255) ggCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc < 1:840534/150‑1 (MQ=255) gTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCg > 1:662381/1‑150 (MQ=255) tGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt > 1:656113/1‑150 (MQ=255) | ACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT > NC_000913/941566‑941862 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |