| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 701,009 | 0 | T | A | 46.2% | 18.7 / 13.0 | 26 | I537N (ATC→AAC) | nagE | PTS N‑acetyl glucosamine transporter subunit IIABC |
| Reads supporting (aligned to +/- strand): ref base T (2/12); new base A (12/0); total (14/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.53e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGATC > NZ_CP009273/700923‑701097 | tGTCGTGGCACTGGATCAGGTTCCTGATGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtcgt < 1:123628/90‑1 (MQ=255) tGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTAt < 1:95580/90‑1 (MQ=255) ggCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtc > 2:1710614/1‑82 (MQ=255) ggCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtc < 1:1710613/82‑1 (MQ=255) ggCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGAAAAAAAAGCCGTATc > 1:1084200/1‑90 (MQ=255) cACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGGGACGGGGTGGCGGGGAAAACGACAGGAAAAATGCTCGTagtgc > 1:804509/1‑86 (MQ=255) cAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCa < 1:1090364/90‑1 (MQ=255) aGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGAAAAAAAAGTCGTATCACCCGCCGGAg > 1:615268/1‑90 (MQ=255) ttCCTGACGAAGAATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGAAAAAAACGTCCTATCACCAACCGCAgggg > 2:859965/1‑89 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAACCCACCGAGAAAAAAGTTGTATAACCCGCCGCGGGGGCa > 1:83248/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCa > 2:154492/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAGATAAAAAAGTCGGATCACCCGCCGCCggggga > 2:1418447/1‑87 (MQ=255) cTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGCGGAACCCACAAGAAAAAACATCGTATCACCAGCCGCCgggggca > 2:414191/1‑86 (MQ=255) tGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAAt < 1:1623960/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 1:334161/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 2:1367847/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAGATAAAAACGTCGGATCACCAGCCGCAGGGACAATCGTGgaaa > 1:1885419/1‑90 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACCACAGATAAAAACGGCGTATCACCAGCCGCCGGGGCAATCGTGgaaa > 1:1732813/1‑90 (MQ=255) cATTCGCCAGCAAAGCGGTGGGGGACGGGGGGCCGGGGAAACCAAAAAAAAAAAACCTGGTTTAACCACCCCCAGGGGAAAACGGGaaaa > 1:886358/1‑90 (MQ=255) ccAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCttaa < 1:1101358/90‑3 (MQ=255) aGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaa < 2:1860665/90‑1 (MQ=255) aCGGTGTGGCGGGGAAAACCACAGATAAAAAAGTCGTATCACCAGCCCCCGGGGCAAACGTGAAAAACTTCAACAACAACCCCGCGGTCt > 1:1802649/1‑90 (MQ=255) gtgtGGCGGTGAAAACGACAAATAAAATCGTCGTATCACCAGCCGCAGGGGCAAACGTGGAAAACTTCAACACCAAACACGCCGTCTGcc > 1:1155696/1‑90 (MQ=255) ggCGGTGAAACCGACAGATATAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGa < 2:1732814/90‑1 (MQ=255) ggCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGa < 2:886359/90‑1 (MQ=255) ggCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGa < 2:1885420/90‑1 (MQ=255) gCGGTGAAAACGACAGAAAAAAACGTCGTATCACCAGCCGCCGGGGCAAACGGGAAAAACTTCAAAACCAACCACCCGGTCTGCCTGGaa > 1:463235/1‑90 (MQ=255) aaaaTCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGgag < 2:1155697/90‑1 (MQ=255) aTCGTCGTATCCCCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGAtc < 1:414191/90‑1 (MQ=255) | TGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGATC > NZ_CP009273/700923‑701097 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |