Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,815,692 | 0 | A | G | 63.8% | 37.3 / 67.5 | 58 | E76G (GAA→GGA) | rpoZ | DNA‑directed RNA polymerase subunit omega |
Reads supporting (aligned to +/- strand): ref base A (10/11); new base G (19/18); total (29/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-02 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGTAGGCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAAT > NZ_CP009273/3815558‑3815818 | ggTAGGCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGcc > 1:160389/1‑139 (MQ=255) ggTAGGCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGcc > 2:627888/1‑139 (MQ=255) ggTAGGCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGcc > 2:549531/1‑139 (MQ=255) ggTAGGCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGcc > 2:473431/1‑139 (MQ=255) gTAGGCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCg < 1:603646/139‑1 (MQ=255) ggCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCTAGGAACAGCAAGAGCAGGAAGCCGCTg < 1:304789/139‑1 (MQ=255) ggCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTg < 1:69413/139‑1 (MQ=255) ggCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCt < 1:602635/138‑1 (MQ=255) gCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTg > 2:235811/1‑138 (MQ=255) gCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTg > 2:778878/1‑138 (MQ=255) gCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTg > 1:759986/1‑138 (MQ=255) aGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAGCCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACa < 1:278290/139‑1 (MQ=255) gATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAg < 2:332093/139‑1 (MQ=255) gATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAg < 2:759986/139‑1 (MQ=255) ccGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCg > 1:156991/1‑139 (MQ=255) gCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGtt < 2:117225/139‑1 (MQ=255) tACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGc > 1:425011/1‑139 (MQ=255) ccGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTa < 1:778878/139‑1 (MQ=255) aaCGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAg > 2:240577/1‑139 (MQ=255) aTAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGgtcg < 2:599053/139‑1 (MQ=255) tAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGgtcgt > 1:41235/1‑139 (MQ=255) aTCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACaaa < 2:351155/139‑1 (MQ=255) aTCGCGCTGCGCGAAATCAAAGAAGGTCTGAGCAACAACCAGATCCTCGACTTTCGCAAACGCCAGGAACACCAACAGCAGGGAGCCGCTGAATTACAAGCCTTTACCGCTATTGCTGAAGGTCGTCGTTAATCACaaa < 1:549531/139‑1 (MQ=255) gcgcTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTAc < 2:746249/104‑1 (MQ=255) gcgcTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTAc > 1:746249/1‑104 (MQ=255) cTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGt < 1:289859/139‑1 (MQ=255) gcgcGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTa > 2:792549/1‑104 (MQ=255) gcgcGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTa < 1:792549/104‑1 (MQ=255) cgcgAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTc < 2:822356/137‑1 (MQ=255) cgcgAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTc < 2:463923/137‑1 (MQ=255) cgcgAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTc < 1:350461/137‑1 (MQ=255) cgcgAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCg < 2:425011/138‑1 (MQ=255) cgcgAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGc > 1:378263/1‑139 (MQ=255) gcgAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGcc > 1:705048/1‑139 (MQ=255) cgAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGccc < 2:762122/139‑1 (MQ=255) aagaagGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCt > 1:567361/1‑139 (MQ=255) ggTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTg < 1:235811/139‑1 (MQ=255) caacCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTg > 2:688229/1‑128 (MQ=255) caacCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTg < 1:688229/128‑1 (MQ=255) caacCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAAt > 2:266456/1‑139 (MQ=255) aacCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAAt > 1:266987/1‑138 (MQ=255) ccTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGAtt < 2:597836/139‑1 (MQ=255) cGACGTTCGCGAACGCCAGGAACAGAAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCa > 1:93570/1‑138 (MQ=255) gACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCaaa > 1:527931/1‑139 (MQ=255) aCGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGa > 1:634773/1‑70 (MQ=255) aCGTTCGCGAACGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGa < 2:634773/70‑1 (MQ=255) gcgAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAacctacc < 2:705048/138‑1 (MQ=255) gcgAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAacctacc < 1:399990/138‑1 (MQ=255) cGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCgg < 2:527931/139‑1 (MQ=255) cGCCAGGAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCg > 1:648912/1‑138 (MQ=255) ggAACAGCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGAc < 1:240577/139‑1 (MQ=255) gCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAAt > 2:522324/1‑138 (MQ=255) gCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAAt > 1:319485/1‑138 (MQ=255) gCAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAAt > 1:21401/1‑138 (MQ=255) gCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAAt > 2:169300/1‑138 (MQ=255) gCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAAt > 1:392554/1‑138 (MQ=255) cAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAAc < 2:566182/106‑1 (MQ=255) cAAGAGCAGGGAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAAc > 1:566182/1‑106 (MQ=255) | GGTAGGCGGAAAGGATCCGCTGGTACCGGAAGAAAACGATAAAACCACTGTAATCGCGCTGCGCGAAATCGAAGAAGGTCTGATCAACAACCAGATCCTCGACGTTCGCGAACGCCAGGAACAGCAAGAGCAGGAAGCCGCTGAATTACAAGCCGTTACCGCTATTGCTGAAGGTCGTCGTTAATCACAAAGCGGGTCGCCCTTGTATCTGTTTGAAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAAT > NZ_CP009273/3815558‑3815818 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |