Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 700,996 | 0 | C | A | 50.0% | 23.7 / 16.5 | 26 | P533T (CCG→ACG) | nagE | PTS N‑acetyl glucosamine transporter subunit IIABC |
Reads supporting (aligned to +/- strand): ref base C (2/11); new base A (13/0); total (15/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.02e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAA > NZ_CP009273/700911‑701085 | gATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGAc < 1:344154/90‑1 (MQ=255) aTTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACa > 2:727567/1‑90 (MQ=255) tGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAt < 2:109094/90‑1 (MQ=255) gATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAtc > 1:337539/1‑90 (MQ=255) gATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGCGGAAACGCCAGATaaaagc > 2:467410/1‑88 (MQ=255) cGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATcgtcgt < 2:504774/90‑1 (MQ=255) ggATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAAAGTCGTATCACCAGc > 2:300813/1‑90 (MQ=255) cAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCa < 1:285670/90‑1 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGGAAACCACAGATAAAAACGTCGTATCACCCCCCCCCGGGACa > 2:68324/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCa > 1:14086/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACCACAGATAAAAAAGTCGTATTACCCGCCGCAGGGGCa > 2:144069/1‑90 (MQ=255) tGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAAAAAAAAAAGTCGGATCACCAGCCGCCGGGGCaaa > 2:761546/1‑89 (MQ=255) cGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGt < 1:477604/90‑1 (MQ=255) gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCCACAGATAAAAACGGCGTATCACCAGCCGCCGGGACAAACggg > 1:367008/1‑88 (MQ=255) gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACCACAGATAAAAACGGCGGATCACCAGCCGCAGGGGCAATAGTg > 2:494240/1‑90 (MQ=255) aaGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGGAAACCACCGATAAAAAAGTCGGATCACCCGCCGCGGGGGCAAACGTGa > 2:578398/1‑90 (MQ=255) aGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaa < 2:732451/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 1:462743/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 2:523635/90‑1 (MQ=255) cGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCtt < 1:82350/90‑1 (MQ=255) ccAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCa < 1:195075/90‑1 (MQ=255) aGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAac < 1:727567/90‑1 (MQ=255) gCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAATCCTCCTATCACCACCCCCAGGGACGAACGGGAAAATCTTCaaaa > 2:872405/1‑88 (MQ=255) aCGGTGTGGCGGTGAAAACGACAAATAAAAAAGTCGTATCCCCAACCCCCGGGGCAAACGGGAAAAAATTTAACACCCACCCcgcggttt > 1:807829/1‑86 (MQ=255) gtgtGGCGGTGAAAACGACAGATAAAATCGTCGTATCACCAGCCCCAGGGGCAAACGTGAAAATCTTCAACACCAACCACGCGTTCTGcc > 1:808209/1‑90 (MQ=255) gtgtGGCGGTGAAAACGACAGATAAAAACGGCGTATCACCCGCCGCAGGGGCAAACGGGAAAAACTTCAACACCAACCACCCGTTCTGcc > 2:174723/1‑90 (MQ=255) ccGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGaaaaa < 1:283448/90‑1 (MQ=255) | GATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAA > NZ_CP009273/700911‑701085 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |