Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 154,497 | 0 | T | G | 32.2% | 121.5 / 14.5 | 59 | A191A (GCA→GCC) | yadL | predicted fimbrial‑like adhesin protein |
Reads supporting (aligned to +/- strand): ref base T (3/37); new base G (19/0); total (22/37) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.10e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATA > REL606/154425‑154571 | gTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCa < 2:122099/76‑1 (MQ=255) tGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGa > 2:133875/1‑76 (MQ=255) aGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGaaa > 1:920198/1‑76 (MQ=255) aGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGaaa > 1:181708/1‑76 (MQ=255) gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaaa < 2:1414722/76‑1 (MQ=255) tGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTTTAAATGTGGCAGGGGCAGAAAAt > 1:1074829/1‑76 (MQ=255) gggATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGGAAAtt > 2:1315844/1‑76 (MQ=255) gggATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAAtt > 2:1315837/1‑76 (MQ=255) gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTa < 2:298727/76‑1 (MQ=255) gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTa < 2:571552/76‑1 (MQ=255) aTGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTAc > 2:1377949/1‑76 (MQ=255) tGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAcc < 2:682277/76‑1 (MQ=255) tGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAcc < 1:422630/76‑1 (MQ=255) tGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAcc < 2:682294/76‑1 (MQ=255) atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTAACCGc > 2:1697984/1‑76 (MQ=255) tatataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCa < 1:545370/76‑1 (MQ=255) atataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAg < 2:227279/76‑1 (MQ=255) atataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGAGGCAGAAAATTTACCCGCAg > 1:1697232/1‑76 (MQ=255) tataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAgg > 2:567130/1‑75 (MQ=255) tataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTAACCGCAGt > 1:295824/1‑76 (MQ=255) taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAGTGa > 2:367713/1‑76 (MQ=255) tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt < 1:1201804/76‑1 (MQ=255) cTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTGc > 2:430483/1‑76 (MQ=255) tGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCt < 2:1302355/76‑1 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCtt < 2:1271220/76‑1 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCtt < 1:927606/76‑1 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTACCCGCAGTGACATTGCtt > 2:123084/1‑76 (MQ=255) gATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTa < 2:866211/76‑1 (MQ=255) gATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTa < 1:925599/76‑1 (MQ=255) gATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTa > 1:1184318/1‑76 (MQ=255) cTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGTGACATTGCTTACgg > 2:21791/1‑76 (MQ=255) tGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTACCCGCAGGGACATTGCTTTCGGt > 1:507411/1‑76 (MQ=255) tGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTGCTTACGGt > 1:1313536/1‑76 (MQ=255) gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtt < 1:984825/76‑1 (MQ=255) gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAGTGACATTGCTTTCgggt > 1:800074/1‑74 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTa < 1:1727079/76‑1 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTa < 2:379557/76‑1 (MQ=255) gTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAt > 1:1424973/1‑76 (MQ=255) gTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAGTGACATTGCTTACGGtttttt > 1:333163/1‑74 (MQ=255) cTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGtttttttt > 2:1635755/1‑76 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt < 2:588612/76‑1 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt < 1:1454903/76‑1 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt < 2:834123/76‑1 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt > 2:384013/1‑76 (MQ=255) gttagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttttt < 2:922374/76‑1 (MQ=255) gttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcg < 1:806024/76‑1 (MQ=255) gttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtttttttttttttg > 2:358772/1‑74 (MQ=255) ttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTTTTTTTTTtga > 1:1478544/1‑76 (MQ=255) tataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgac < 2:1588002/76‑1 (MQ=255) tataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgac < 1:1342062/76‑1 (MQ=255) ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt < 2:1486487/51‑1 (MQ=255) ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacg < 1:233753/76‑1 (MQ=255) ataGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACAtt > 1:1486487/1‑51 (MQ=255) taGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCCGTGACATTTCTTTCGGGTTTTTTTTTTTttacga > 2:789474/1‑76 (MQ=255) ggTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGAcc < 1:516844/76‑1 (MQ=255) ggTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGAcc < 2:769172/76‑1 (MQ=255) tGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGa < 2:967147/76‑1 (MQ=255) aaaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCa < 1:41163/76‑1 (MQ=255) aaaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCa < 2:761547/76‑1 (MQ=255) tgtgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTc < 1:1345528/76‑1 (MQ=255) tgtgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTc < 2:1384178/76‑1 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt < 2:306499/76‑1 (MQ=255) ggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATc < 2:624598/76‑1 (MQ=255) aGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGa > 2:848890/1‑76 (MQ=255) ggTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGat < 1:424678/76‑1 (MQ=255) gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata > 1:299460/1‑76 (MQ=255) gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata > 1:1599486/1‑76 (MQ=255) | GTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATA > REL606/154425‑154571 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |