Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 700,996 | 0 | C | A | 54.5% | 20.4 / 35.3 | 22 | P533T (CCG→ACG) | nagE | PTS N‑acetyl glucosamine transporter subunit IIABC |
Reads supporting (aligned to +/- strand): ref base C (0/10); new base A (12/0); total (12/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.55e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCC > NZ_CP009273/700909‑701071 | ccGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAcca > 1:1360341/1‑89 (MQ=255) aTTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGAcc > 1:1534407/1‑89 (MQ=255) cTGGATCAGGTTCCTGACGAAGCATTCGACAGCAAGGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCa < 1:1394233/90‑1 (MQ=255) tGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAg < 1:1378560/90‑1 (MQ=255) ggATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGc < 2:1150105/90‑1 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGAAAAAAACGTCGTATTACCCGCCGCCGGGGCa > 1:220118/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACCGCAGAAAAAAACGTCGTATCCCCAGCCGGAGGGGCa > 1:635173/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACCACAAAAAAAAACGTCGTATCACCCGCCCCCGGGGCa > 1:1479975/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAGAAAAAAAAGTCGTATCACCAGCCGCCGGGGCa > 1:795195/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACAACAAAAAAAAAAGTCGGATCACCAGCCCCCGGGGCa > 1:1380080/1‑90 (MQ=255) cTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACaa < 2:102579/90‑1 (MQ=255) gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGGCAGATAAAAACGTCGTATCACCAGCCGCAGGGACAAACggg > 2:791810/1‑88 (MQ=255) gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAGAAAAAAAAGTCGGATCACCCACCGCCGGGGCAAACGTg > 1:417537/1‑90 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 1:615389/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 1:324826/90‑1 (MQ=255) cGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGGAAACCACAGATAAAAACGTCGTATCACCCGCCGCAGGGGCAATCGGGGAAAACtt > 2:680259/1‑90 (MQ=255) gCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTc < 2:752490/90‑1 (MQ=255) ccAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCa < 2:1360341/90‑1 (MQ=255) aaaGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACacc < 1:1505509/90‑1 (MQ=255) gggtgACGGTGTGGCGGTGAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGGCAATCGTGAAAAACTTCAACACCAACCAcgc > 2:893404/1‑90 (MQ=255) gtgACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGt < 2:1534407/90‑1 (MQ=255) gACGGTGTGGCGGGGGAAACCACAGAAAAAAAAGTCGTATCACCAGCCGCGGGGGGAAACGTGGAAAAATTTAACAACAACCACCCGGTc > 1:887696/1‑90 (MQ=255) gtgtGGCGGGGAAAACGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAATCGGGAAAAAATTTAAAACCAACCACCCGGTCTGcc > 1:1226527/1‑90 (MQ=255) | CCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCC > NZ_CP009273/700909‑701071 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |