| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 700,996 | 0 | C | A | 43.2% | 61.2 / 12.9 | 45 | P533T (CCG→ACG) | nagE | PTS N‑acetyl glucosamine transporter subunit IIABC |
| Reads supporting (aligned to +/- strand): ref base C (5/20); new base A (19/0); total (25/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.34e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAA > NZ_CP009273/700907‑701084 | cGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAc > 1:1248843/1‑90 (MQ=255) aTTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAACCGACa > 2:2269476/1‑90 (MQ=255) gATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtc < 2:1248843/90‑1 (MQ=255) tGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGCCAGaaaaaaaagg > 1:26652/1‑86 (MQ=255) tGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGaaaaaaaagt > 1:1081706/1‑86 (MQ=255) cGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGATAAAAAcgtcgt > 1:1174538/1‑90 (MQ=255) cGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGACAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCAACAGATAAAAAagtcgt > 2:2200358/1‑90 (MQ=255) cGTGGCACTGGATCAGGTTCCTAACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAGATAAAAaagtcgg > 1:2027/1‑89 (MQ=255) ggCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGATAAAATAGTCGTATc > 1:2469328/1‑90 (MQ=255) cTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAGATAAAAAAGTCGTTTCAccc > 1:2232020/1‑89 (MQ=255) tGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAg < 1:1617473/90‑1 (MQ=255) tGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAg < 1:764463/90‑1 (MQ=255) tGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAg < 2:1399512/90‑1 (MQ=255) aTCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCg < 2:1490002/90‑1 (MQ=255) cAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCa < 2:1679506/90‑1 (MQ=255) cAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCa < 1:676223/90‑1 (MQ=255) gTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAggg < 1:2269476/90‑1 (MQ=255) ttCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGAAAAAAAAGTCGTATTACCCGCCGCAgggg > 1:2303057/1‑89 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACa < 1:382855/90‑1 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGGCa > 2:1546786/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAAAGAAAAAAGTCCGTTCACCACCCCCAGGGACa > 1:1807126/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACCAAAGATAAAAACGTCTTATCACCAGCCGCGGGGGCa > 2:614044/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGAAAAAAAAGTCGTTTCACCAGCCCCAGGGGCa > 1:1695978/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAAAAAAAATCCTCGTATATCCCGCCGCCGGGGCa > 1:1548864/1‑90 (MQ=255) gaagaaTTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGGCAGATAAAAAAGACGTGTCACCCGCCGCCGGGGGAATCGTg > 1:326284/1‑90 (MQ=255) aaGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAATGTCGTATCACCAGCCCCAGGGGCAATCGTGa > 2:2218324/1‑90 (MQ=255) aaGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACCACAGAAAAAAACGTCGTTTCACCCGCCGCCGGGGCGATCCTgg > 1:2162374/1‑89 (MQ=255) aGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaa < 2:1081706/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 1:378464/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAATCGCCAGATAAAAACGTCGGATTACCATCCGCAGGGGCCAGCGTGaaaa > 2:1782536/1‑90 (MQ=255) cATTCGCCAGCAAAGCGGTGGGCGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 1:2370600/90‑1 (MQ=255) aGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCCCCAGCCGCAGGGACAATCGTGAAAATCTTCAac < 2:1807126/90‑1 (MQ=255) aGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAac < 2:1548864/90‑1 (MQ=255) gCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGAAAAAAAAGTCGTATCCCCCACCGCCGGGAGAAACGTGAAAATCTTCAACaccaac > 1:1877167/1‑90 (MQ=255) cGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaacc < 1:2218324/90‑1 (MQ=255) gggtgACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCCCAGGGGCAATCGTGAAAAACTTCAAAACCAACCAcgc > 2:351216/1‑90 (MQ=255) gggtgACGGTGTGGCGGTGAAAACGACAGATAAAAACGTCGTATCACCAGCCGCAGGGACAATCGTGAAAAACTTCAACACCAACCAcgc > 2:44017/1‑90 (MQ=255) gtgACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGt < 2:1225288/90‑1 (MQ=255) aCGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCt < 2:1905254/90‑1 (MQ=255) aCGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAAAACCAACCACGCGTTCt < 2:326284/90‑1 (MQ=255) cGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGaaa < 2:1695978/90‑1 (MQ=255) gTGAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGGCAATCGGGAAAAACTTCAACACCAACCACCCGGTCTGCCTGGaaaac > 1:2043780/1‑88 (MQ=255) gTGAAAACGACAGATAAAAACGTCGTATCACCAGCCGAAGGGGGAATCGTGGAAAACTTCAACACCAACCACGCGTTCTGCCTGGgaaac > 1:1326378/1‑88 (MQ=255) gAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGGCAATCGTGAAAAAATTTAAAAAAAAACACCCGTTCTGCCTGGAAAACGa > 1:1864505/1‑90 (MQ=255) gAAAACGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAAACGTGAAAAAATTTAAAACCAACCACGCGGTCTGCCTGGAAAACGa > 1:565767/1‑90 (MQ=255) aCCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGaaaa < 2:1128685/90‑1 (MQ=255) | CGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAA > NZ_CP009273/700907‑701084 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |