Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,882,526 | 0 | T | G | 35.0% | 96.7 / 13.1 | 60 | I21S (ATT→AGT) | manZ | mannose‑specific enzyme IID component of PTS |
Reads supporting (aligned to +/- strand): ref base T (27/12); new base G (0/21); total (27/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.81e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACTGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCTG > REL606/1882451‑1882601 | aaCTGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGatat < 1:841320/76‑1 (MQ=255) aCTGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAtt < 1:488345/76‑1 (MQ=255) tGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGCGATATTCg > 1:152315/1‑76 (MQ=255) aCTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTgg > 1:1380839/1‑76 (MQ=255) cTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGc < 2:398995/76‑1 (MQ=255) cTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGc < 1:372378/76‑1 (MQ=255) tAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCg > 1:1110035/1‑76 (MQ=255) aCAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTc < 1:663369/76‑1 (MQ=255) aGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCtt < 2:1023390/76‑1 (MQ=255) ggTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTc < 2:1296397/76‑1 (MQ=255) gTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTcc > 1:1525778/1‑76 (MQ=255) gTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTcc > 1:1578957/1‑76 (MQ=255) tGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCt < 2:82802/76‑1 (MQ=255) tGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCt < 2:7858/76‑1 (MQ=255) gAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCGTg < 2:1198893/76‑1 (MQ=255) gAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGtt < 1:1328059/76‑1 (MQ=255) gAAATGGCTGATACAACTCAAACTACCACCGACAAAAAACTCACTCAAACCGATATTCGCGGCGCCTCCCTGCGtc > 2:72631/1‑75 (MQ=255) aaaTGGTTGATACAACTCAAACTACCACCGAGAAACAACTCACTCAAAGTGATATTCGTGGCGTATTCCTGCGTTc > 2:1651748/1‑76 (MQ=255) aaTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCt > 2:1089146/1‑76 (MQ=255) aaTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCt < 2:1085532/76‑1 (MQ=255) aaTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCt < 2:479349/76‑1 (MQ=255) tGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTaa > 2:450547/1‑76 (MQ=255) tGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTaa > 2:458572/1‑76 (MQ=255) tGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTaa > 1:255431/1‑76 (MQ=255) ggTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAAc > 1:221811/1‑76 (MQ=255) ggTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAAc < 2:290360/76‑1 (MQ=255) ttGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACct > 1:577358/1‑76 (MQ=255) aTACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCtt > 1:1042061/1‑76 (MQ=255) tACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTc > 2:636078/1‑76 (MQ=255) tACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTc < 2:1659784/76‑1 (MQ=255) aaCTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAg < 2:20514/76‑1 (MQ=255) aaCTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAg < 2:1215978/76‑1 (MQ=255) cAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGtt < 2:759670/76‑1 (MQ=255) cAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGtt < 2:517897/76‑1 (MQ=255) aaaCTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTc > 1:685678/1‑76 (MQ=255) accaccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGa > 2:606199/1‑76 (MQ=255) accaccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGa < 2:1610560/76‑1 (MQ=255) accaccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGa > 1:523637/1‑76 (MQ=255) accaccGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGa < 2:1261952/76‑1 (MQ=255) ccaccGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGaa < 1:407451/76‑1 (MQ=255) ccaccGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGaa < 1:1306273/76‑1 (MQ=255) ccaccGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGaa < 2:1256201/76‑1 (MQ=255) accGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACt > 1:886966/1‑76 (MQ=255) gaAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGaa < 2:556808/76‑1 (MQ=255) aaaaaaCTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAAc < 2:1060305/76‑1 (MQ=255) aaaaaCTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACg < 1:1411316/76‑1 (MQ=255) aaaaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGt < 2:629193/76‑1 (MQ=255) aaaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTa > 2:1015541/1‑76 (MQ=255) aaaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTa > 1:18652/1‑76 (MQ=255) aaaCTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTa < 2:1357642/76‑1 (MQ=255) aaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTAt > 2:1630006/1‑76 (MQ=255) aaCTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTAt < 1:1126562/76‑1 (MQ=255) ctcactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGc < 2:1051745/76‑1 (MQ=255) cactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAg > 2:838310/1‑76 (MQ=255) cactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAg > 2:122641/1‑76 (MQ=255) actcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAgg > 1:961514/1‑76 (MQ=255) ctcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGc > 2:76401/1‑76 (MQ=255) ctcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGc < 1:613717/76‑1 (MQ=255) tcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCa < 1:225498/76‑1 (MQ=255) tcaAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCa < 1:764003/76‑1 (MQ=255) tcaAAGTGATAGTCGTGGCGTCTTCCTGCGTCCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCa < 2:717839/76‑1 (MQ=255) aaaGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACt > 2:1416584/1‑76 (MQ=255) aaaGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACt > 2:24978/1‑76 (MQ=255) aaGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTg > 1:809859/1‑76 (MQ=255) aaGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTg < 1:918951/76‑1 (MQ=255) atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt > 1:41875/1‑76 (MQ=255) atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt > 1:788255/1‑76 (MQ=255) atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt < 2:1655894/76‑1 (MQ=255) atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt < 2:1472561/76‑1 (MQ=255) atagTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt < 1:857213/76‑1 (MQ=255) atagTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt < 2:979026/76‑1 (MQ=255) tatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTc > 1:632212/1‑76 (MQ=255) ttCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCTg < 1:1321529/76‑1 (MQ=255) | AACTGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCTG > REL606/1882451‑1882601 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |