Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,882,526 | 0 | T | G | 26.1% | 134.6 / 10.0 | 69 | I21S (ATT→AGT) | manZ | mannose‑specific enzyme IID component of PTS |
Reads supporting (aligned to +/- strand): ref base T (36/15); new base G (0/18); total (36/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.02e-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:431330/76‑1 (MQ=255) aCTGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAtt < 2:1026901/76‑1 (MQ=255) aCTGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAtt > 2:887935/1‑76 (MQ=255) ggACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGt > 2:982113/1‑76 (MQ=255) ggACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGt < 1:844714/76‑1 (MQ=255) ggACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGt > 1:1842773/1‑76 (MQ=255) ggACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGt < 1:844687/76‑1 (MQ=255) gACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTg < 1:652933/76‑1 (MQ=255) aCTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTgg > 2:16038/1‑76 (MQ=255) aCTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTgg > 2:1249446/1‑76 (MQ=255) cTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGc > 1:1159499/1‑76 (MQ=255) tAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCg > 2:1849247/1‑76 (MQ=255) aaCAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGt > 1:1018417/1‑76 (MQ=255) cAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCt > 2:1225040/1‑76 (MQ=255) cAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCt < 2:640549/76‑1 (MQ=255) aGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCtt > 1:972738/1‑76 (MQ=255) gTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTcc > 1:454856/1‑76 (MQ=255) tGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCt > 1:1209880/1‑76 (MQ=255) gAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTg > 1:845825/1‑76 (MQ=255) gAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTg > 1:1698250/1‑76 (MQ=255) gCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCg < 2:416271/76‑1 (MQ=255) gCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCg < 2:1897194/76‑1 (MQ=255) gAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGtt > 1:231037/1‑76 (MQ=255) aaaTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTc > 1:687739/1‑76 (MQ=255) aaaTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTc > 1:1634096/1‑76 (MQ=255) aTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTa > 1:1226424/1‑76 (MQ=255) aTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTa < 1:109367/76‑1 (MQ=255) aTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTa < 1:1288862/76‑1 (MQ=255) tGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTaa < 2:1733937/76‑1 (MQ=255) ggTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAAc < 1:1661421/76‑1 (MQ=255) gTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAAcc > 1:1802107/1‑76 (MQ=255) tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc > 1:494038/1‑76 (MQ=255) tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc > 2:79111/1‑76 (MQ=255) tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc > 2:919971/1‑76 (MQ=255) tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc > 1:1982616/1‑76 (MQ=255) gATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCCGCGTTCTAACctct > 2:306254/1‑76 (MQ=255) aTACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCtt < 2:1344398/76‑1 (MQ=255) tACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTc > 2:1058289/1‑76 (MQ=255) aCAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTcc < 1:919971/76‑1 (MQ=255) cAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCa > 2:410303/1‑76 (MQ=255) cAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCa < 1:571291/76‑1 (MQ=255) aaCTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAg > 2:1932669/1‑76 (MQ=255) cTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAggg < 2:692493/76‑1 (MQ=255) cAAACTACCCCCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGtt < 2:579014/76‑1 (MQ=255) cAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGtt > 1:1226692/1‑76 (MQ=255) cAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGtt < 2:517130/76‑1 (MQ=255) aaaCTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTc > 2:1348238/1‑76 (MQ=255) cTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATg < 1:1671465/76‑1 (MQ=255) cTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCCCTTCCACAGTTCATg > 1:1366120/1‑76 (MQ=255) tACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATgg > 1:504081/1‑76 (MQ=255) accaccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGa > 1:845725/1‑76 (MQ=255) accaccGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGa < 1:651575/76‑1 (MQ=255) caccGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAAc < 1:1953359/76‑1 (MQ=255) accGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACt < 2:1223515/76‑1 (MQ=255) ccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACtt < 2:1603705/76‑1 (MQ=255) cGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTc > 1:309317/1‑76 (MQ=255) cGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTc < 2:1620280/76‑1 (MQ=255) gagaAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCg < 2:59895/76‑1 (MQ=255) agaAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGa > 2:1930766/1‑76 (MQ=255) aaaaaaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAAc < 2:1149699/76‑1 (MQ=255) aaaaaaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAAc < 2:1761548/76‑1 (MQ=255) aaaaaaCTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAAc < 1:1356404/76‑1 (MQ=255) aaaaCTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGt < 2:1774056/76‑1 (MQ=255) aaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGATCATGGAACTTCGAACGTAt > 1:1024656/1‑76 (MQ=255) actcactcAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATg > 2:561673/1‑76 (MQ=255) actcactcAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATg < 1:104011/76‑1 (MQ=255) ctcactcaAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGc < 1:1275679/76‑1 (MQ=255) ctcactcaAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGc < 1:1379399/76‑1 (MQ=255) tcactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCa > 1:1241401/1‑76 (MQ=255) tcactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCa > 1:1752949/1‑76 (MQ=255) aaGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTg > 1:1563569/1‑76 (MQ=255) aGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTgg < 2:1358766/76‑1 (MQ=255) aGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTgg < 1:1052505/76‑1 (MQ=255) aGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTgg < 2:112055/76‑1 (MQ=255) tGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGt > 1:163501/1‑76 (MQ=255) atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt < 2:859984/76‑1 (MQ=255) atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt > 2:1905695/1‑76 (MQ=255) atagTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt < 2:360565/76‑1 (MQ=255) tatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTc > 1:552928/1‑76 (MQ=255) tatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTc < 2:1408659/76‑1 (MQ=255) atTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCCGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCt > 2:1295391/1‑76 (MQ=255) atTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCt < 1:1722501/76‑1 (MQ=255) atTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCt > 1:1793106/1‑76 (MQ=255) ttCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCTg > 1:564632/1‑76 (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 < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |