| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 941,713 | 0 | T | G | 52.9% | 33.0 / 51.6 | 35 | V252G (GTG→GGG) | dmsA | dimethyl sulfoxide reductase, anaerobic, subunit A |
| Reads supporting (aligned to +/- strand): ref base T (1/15); new base G (18/0); total (20/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.62e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-06 | |||||||||||
GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG > NC_000913/941563‑941859 | gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGt < 1:3724359/151‑1 (MQ=255) cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTAc < 2:3335409/149‑1 (MQ=255) tccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTACt > 1:2444930/1‑148 (MQ=255) ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATc < 2:3052644/150‑1 (MQ=255) gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGGATGAGTGGC‑GGTGGGGTGACTTACTATCttt < 1:2971857/151‑2 (MQ=255) gaaaTGGAGGAAGGTTTGAACGGGAAAAACGGGGGCCGGTTAGTTGGCAACAGCGGGTCGGCGATCGAAAACAAAAAACTAGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGCGTCTTATTATCTTTaa > 1:1342941/2‑150 (MQ=255) aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATCTTTaa > 2:989140/1‑149 (MQ=255) gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATCTTGAACAgg > 1:3589682/1‑148 (MQ=255) cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACg < 1:3760850/151‑1 (MQ=255) aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTTCTATCTTGAAAAAGCCCGc > 1:3084023/1‑149 (MQ=255) aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATATTTAACAGGCACGcc > 1:3995623/1‑149 (MQ=255) ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑Tg < 2:3589878/150‑1 (MQ=255) cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTACTATCTTTAAAAAGGCCGCCAGGAAATTAA‑TGACCGc > 1:770123/1‑148 (MQ=255) cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGATTTCTATCTTGAAAAAGCACGCCAGAAAAATAA‑TGCCCCCATGATCCTCATCg > 1:3912658/1‑150 (MQ=255) gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCcgcg < 2:2838490/151‑1 (MQ=255) cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGTGGGGGGA‑TTACTTTTTTTAAAAGGCCAGCCAAAAAATTAATTGCACCCA‑GGTCATCATCGACTCCgcgcct > 1:3876643/1‑146 (MQ=255) cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATAcc < 1:1645891/151‑1 (MQ=255) cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACCTAAAATCCTTAAAAAGCACGGCAGAAAACTAA‑AGCCCGCCAGTTCCTCATCGA‑TCCGCCCCATAACGACAACGGTGcc > 1:4089713/1‑151 (MQ=255) gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGcca < 1:75599/151‑2 (MQ=255) gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGCGGGGGGGCTTACTAATTTTTACAGGCCCGCCAGGAATCTAA‑ATCCCGCAATATTCTCCTCGA‑TCCCCGCCATACCGACAACggctcc > 2:2335811/1‑146 (MQ=255) aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTTAAAAAGCACGCCAGAAATCTAA‑TGCCCCCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGGGcc > 2:2789363/1‑149 (MQ=255) aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGCCggg < 1:1540588/151‑1 (MQ=255) aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTTCCGGgc < 2:2355626/151‑1 (MQ=255) aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGCCGGGcgcg < 1:1621689/151‑1 (MQ=255) aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGCCGGGCGCGAAg < 2:3508089/151‑1 (MQ=255) aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATTTTGAAAAGGCACGCCAGAAATCTAA‑TGCCCGCCTGATCATCATCGA‑TCCGCGCTATAACGACACCGGGGCCGGGcgcg > 2:2791768/1‑148 (MQ=255) tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGGGCCGGGGGCGAAGATg > 2:3017362/1‑150 (MQ=255) aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTgg < 2:2477479/151‑1 (MQ=255) ataaCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAA‑TGCCCGCCTGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGCCGGGGGCGAAGATGAGTGGATCCCTATTCGTcc > 2:3760693/1‑151 (MQ=255) cTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCCTGAACAAGCACGCCCGAAAACTAA‑AGCCCGCATGAACATCATCGA‑TCCGCGCGATACCGGCACCCGTGCCGGGCGCGGAAGTGAGTGGGTCCCTATTTGTCCgggga > 1:1291278/1‑148 (MQ=255) ggCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc < 1:4030611/149‑1 (MQ=255) gAAACGCGAATGAGTGGC‑GGTGGGGGGACTTTCTATCTTTAACAGGCCAGCCAGAAATCTAA‑TGCCCGCCTGATCATCATCGA‑ACCGCGCTATACCGAAACCCGGGGCGGGGGCGGAGAAGAGTGGATCCCTATTCGTCCCGGAACAGATg > 1:2185899/1‑151 (MQ=255) aaaCGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGc < 2:2315487/151‑1 (MQ=255) cgcgAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGc > 2:2053079/1‑151 (MQ=255) gAGTGGC‑GGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGAACATCATCGA‑ACCGCGCTATACCGACACCCGGGCCGGGCGCGAAGATGAGTGGGTCCCCATTCGTCCGGGAACAGATGCCGCACTGGtt > 1:3514553/1‑151 (MQ=255) ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg > 1:905704/1‑151 (MQ=255) | GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAA‑TGCCCGCATGATCATCATCGA‑TCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG > NC_000913/941563‑941859 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |