| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,369,508 | 0 | C | G | 43.6% | 63.5 / 10.1 | 39 | A39G (GCG→GGG) | pspD | peripheral inner membrane phage‑shock protein |
| Reads supporting (aligned to +/- strand): ref base C (4/18); new base G (17/0); total (21/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.81e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGAT > NC_000913/1369368‑1369658 | gTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGGTaaaa > 2:598630/1‑148 (MQ=255) tCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGGGTGGCGGGCGGGGCGATATAAACAGt > 1:503405/1‑150 (MQ=255) ccGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGGGGGGGGTaaaaaaaagt > 1:1573/1‑144 (MQ=255) cGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGGGTGGCGGGGTGGGGGGCTAAATCAGtt > 1:483647/1‑149 (MQ=255) gTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGATaaaaaaaggt > 2:757043/1‑142 (MQ=255) tCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGGTAAAATAAggt > 2:65579/1‑145 (MQ=255) cTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGGTAAAAACAGTTTCTcccc > 2:386754/1‑150 (MQ=255) gAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGAcc < 2:568370/150‑1 (MQ=255) ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGTCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCt < 1:689230/151‑1 (MQ=255) aTACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTgct < 2:735383/151‑1 (MQ=255) cGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTgctggctg < 1:50143/151‑1 (MQ=255) gCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCCGGGGGGGAAAAAAAGGTGCTCGCCCCCCCCCGGAAATGTTgctgg > 2:596756/1‑147 (MQ=255) cTGGCTACAGGCCGGGCAAAATATAAAGCCTTGTTTCAAATTAGCAGCCAAGCTGATACTTCTTACCACACTGCACTATGGCCCGGCGGGTGTGGCTTGCTGTGCTATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTgctggctt < 1:50154/149‑2 (MQ=255) ggCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCtgtg < 2:796825/149‑1 (MQ=255) aGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGGGGGAAAAAAAAGTGGTCGCCGACCGCTGGAAATGTTGGTGGGTGTGGGG‑GT‑ggg > 2:226075/1‑150 (MQ=255) aaaGGTAAAGGCTAGTTTCAAATTGTCAGGCAAGCTGGTACTTCGTACCTCACTTGGCCATGGCCCGGCGGGTGTGGCGCGCTGGGCGACAAAATCCGTTCCTCGCCCGCCGCTGAACATGTTGCTGGCTGTGGCG‑CT‑TGACCCGCTGTTa < 1:270490/151‑1 (MQ=255) ggTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGGTAAAATTAATTGGTCCCCGACCGCCGGAAAAGGTGCTGGCTGTGGGG‑CC‑GGGAACGCTGttta > 2:201553/1‑147 (MQ=255) tAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGGGGTAAAAACAGTTTGTCGCCCACCCCCGGAAAAGTTGGTGGCTGTGGGG‑CC‑GGGAACGCCGttta > 2:647576/1‑145 (MQ=255) aaGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGGGGAAAAAACAAGTTCTCCCCCGCCCCCGAAAAAGTTGTTGGCTGTGGGGCCG‑GGAACCGC‑GTTAAAGCggg > 2:637994/1‑151 (MQ=255) cTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGATAAAAAAAGTTGTCCCCCGCCCGCTGGAAATGTTTCTGGGTGTGGGG‑CC‑GGGACCGCTGGTaaaatggg > 2:496866/1‑141 (MQ=255) ggTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCta < 2:8030/151‑1 (MQ=255) aaGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGGGAAAAAACGTTTGTCGCCCATCCCCGAAAAAGTTTTTGGCTGTTGGG‑CT‑CGAATCCCTTTTAAGTTCGGGTGCTAACAAAAAGGGACAGC‑Gtt > 2:580090/1‑151 (MQ=255) aGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCCGGGGGGTAAAAATAGGTGGTCTCCGGCCGCCGGAAAAGTTTCTGGCTGTGGGG‑CTGGGAACCCCGGTTAA‑TAGGGGTGCTTATAAAATGGCACAGC‑ggt > 2:440231/1‑148 (MQ=255) gCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCCGGGGGGTAAAAACAGGTGCTCCCCCACCCCTGAAAATGTTGCTGGGTGTGGGG‑CT‑GGAAACGCTCTTTAGTCGGGGGGCTaaaaaaaagggccagccgttt > 2:496617/1‑136 (MQ=255) aCTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGGGGTAAAAACAGTTTCTCGCCCACCCCCGGAAAAATTTCTGGGTGTGGGG‑CT‑GGGACCGCTGTTTAGACGGGGTGCTAAAAAAATGGGCCAGC‑GGTATaaaa > 2:88140/1‑149 (MQ=255) aCTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGGTAAAAACAGTTGCCCGCCGACCGCTGAAAAAGTTGGTGGCTGTGGGG‑CT‑GGAACCGCTGTTTAGTCGGGGCGCTAATAAAATGGGCCCAC‑GTTTTAAAAGgggggggggggtt > 2:665573/1‑146 (MQ=255) cTATGGCCCGGCGGGTGTGGCGGGGTGGGGGGTAAAAACAGGTGCTCCCCCACCGCTGGAAAAGTTGGTGGCTGTGGGG‑CC‑GGGAACCCCGTTTAGTCGGGGTGCTAAAAAAAATGGACAGCGGTTTTAAAAGGGGGGGGGGGtttg > 2:610008/1‑146 (MQ=255) ggCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAAt < 2:31591/150‑1 (MQ=255) cccGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGACACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGt < 2:174272/151‑1 (MQ=255) ggTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAAtct < 2:656722/151‑1 (MQ=255) gtgtGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAAtctc > 2:257673/1‑151 (MQ=255) gtgtGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAAtctc < 2:93221/151‑1 (MQ=255) gtgtGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCAGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAAtctc < 1:519418/151‑1 (MQ=255) ggCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAgg < 2:666673/150‑1 (MQ=255) gCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCg < 2:702435/151‑1 (MQ=255) cGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCg < 1:369882/150‑1 (MQ=255) cGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGt > 2:787961/1‑151 (MQ=255) gggCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATaa < 2:145267/149‑1 (MQ=255) ggCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAatgg < 1:64952/151‑1 (MQ=255) cGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAatggat > 1:278351/1‑151 (MQ=255) cGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAatggat < 2:601994/151‑1 (MQ=255) | GTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGAT > NC_000913/1369368‑1369658 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |