Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,882,526 | 0 | T | G | 22.8% | 175.7 / 10.9 | 79 | I21S (ATT→AGT) | manZ | mannose‑specific enzyme IID component of PTS |
Reads supporting (aligned to +/- strand): ref base T (44/17); new base G (0/18); total (44/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.65e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCT > REL606/1882455‑1882600 | ggACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGt > 2:105776/1‑76 (MQ=255) gACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTg < 1:1819492/76‑1 (MQ=255) gACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTg < 1:1500463/76‑1 (MQ=255) aCTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTgg < 1:896149/76‑1 (MQ=255) aCTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTgg > 1:1668745/1‑76 (MQ=255) cTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGc < 1:1207933/76‑1 (MQ=255) tAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCg > 1:1625997/1‑76 (MQ=255) aaCAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATGCGTGGCGt > 1:1320020/1‑76 (MQ=255) cAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCt > 1:1015409/1‑76 (MQ=255) aGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCtt > 1:400955/1‑76 (MQ=255) aGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCtt < 2:1343877/76‑1 (MQ=255) ggTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTc < 2:520335/76‑1 (MQ=255) ggTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTc < 2:432578/76‑1 (MQ=255) gTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTcc > 1:986400/1‑76 (MQ=255) tGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCt > 1:476333/1‑76 (MQ=255) gAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTg < 2:262029/76‑1 (MQ=255) gCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCg > 2:700332/1‑76 (MQ=255) cGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGt > 1:1248171/1‑76 (MQ=255) cGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGt > 2:1504443/1‑76 (MQ=255) cGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGt > 2:1696595/1‑76 (MQ=255) gAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGtt > 1:776961/1‑76 (MQ=255) gAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGtt > 1:748520/1‑76 (MQ=255) aaTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCt < 2:1019868/76‑1 (MQ=255) aTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTa < 2:601763/76‑1 (MQ=255) aTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTa < 2:2059063/76‑1 (MQ=255) tGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTaa > 1:1829302/1‑76 (MQ=255) tGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTaa < 1:452905/76‑1 (MQ=255) tGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTaa < 1:1576788/76‑1 (MQ=255) ttGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACct > 2:515334/1‑76 (MQ=255) ttGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACct < 1:626121/76‑1 (MQ=255) tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAAcccc > 1:1278530/1‑74 (MQ=255) tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc > 2:1085553/1‑76 (MQ=255) tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc > 1:910600/1‑76 (MQ=255) tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc > 2:638588/1‑76 (MQ=255) tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc > 2:1044198/1‑76 (MQ=255) gATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctct > 1:1223217/1‑76 (MQ=255) gATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctct > 1:1554432/1‑76 (MQ=255) aTACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCtt > 1:414495/1‑76 (MQ=255) aTACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCtt > 2:1227914/1‑76 (MQ=255) aTACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCtt < 2:84224/76‑1 (MQ=255) tACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTc < 2:909599/76‑1 (MQ=255) aCAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTcc > 2:1862781/1‑76 (MQ=255) aCAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTcc > 2:70610/1‑76 (MQ=255) aCAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTcc < 2:2009466/76‑1 (MQ=255) cAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCa < 2:27449/76‑1 (MQ=255) cAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCa < 2:1389481/76‑1 (MQ=255) cAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCa < 2:388093/76‑1 (MQ=255) aaCTCAAACTGCCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAg < 2:1887147/76‑1 (MQ=255) aCTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAgg > 2:968032/1‑76 (MQ=255) cTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAggg > 1:1626544/1‑76 (MQ=255) cTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAggg < 2:2016697/76‑1 (MQ=255) cAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGtt > 2:64147/1‑76 (MQ=255) aaaCTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTc > 1:1206061/1‑76 (MQ=255) aaaCTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTc > 2:2026451/1‑76 (MQ=255) aaaCTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTc < 2:1914727/76‑1 (MQ=255) aaaCTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTc < 1:457554/76‑1 (MQ=255) aCTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCAt < 1:1040578/76‑1 (MQ=255) cTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATg < 2:1233846/76‑1 (MQ=255) tACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATgg > 2:49313/1‑76 (MQ=255) ccaccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGaa > 2:1079636/1‑76 (MQ=255) caccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATgggcc > 1:1306093/1‑73 (MQ=255) caccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAAc > 2:1721996/1‑76 (MQ=255) ccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACtt < 2:670647/76‑1 (MQ=255) ccGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACtt < 2:1110708/76‑1 (MQ=255) gaAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGaa < 1:1833649/76‑1 (MQ=255) aaaaaaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAAc > 2:849737/1‑76 (MQ=255) aaaaCTCACTCAAAGTGAGATTCGGGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGt < 2:307128/76‑1 (MQ=255) aaCTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTAt < 2:234630/76‑1 (MQ=255) tcactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCa > 1:980341/1‑76 (MQ=255) tcactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCa < 2:887697/76‑1 (MQ=255) tcactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCa < 1:1584378/76‑1 (MQ=255) cactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAg > 1:804372/1‑76 (MQ=255) cactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATCGAACTTCGAACGTATGCAg > 2:1352867/1‑76 (MQ=255) cactcaAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAg < 2:730700/76‑1 (MQ=255) caAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCAc < 2:2086847/76‑1 (MQ=255) aaGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTg > 1:360412/1‑76 (MQ=255) gTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTggg > 2:916460/1‑76 (MQ=255) tGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGt > 1:1737746/1‑76 (MQ=255) gATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGtt > 1:1942931/1‑76 (MQ=255) atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt > 2:1101135/1‑76 (MQ=255) atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt > 2:128452/1‑76 (MQ=255) atagTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt < 2:2212649/76‑1 (MQ=255) atagTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt < 2:2212668/76‑1 (MQ=255) tatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTc < 2:2161165/76‑1 (MQ=255) tatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTc > 2:1637294/1‑76 (MQ=255) atTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCt < 2:1657364/76‑1 (MQ=255) atTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCt > 1:1150778/1‑76 (MQ=255) | GGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCT > REL606/1882455‑1882600 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |