| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 154,497 | 0 | T | G | 21.9% | 194.7 / 10.3 | 73 | A191A (GCA→GCC) | yadL | predicted fimbrial‑like adhesin protein |
| Reads supporting (aligned to +/- strand): ref base T (6/51); new base G (16/0); total (22/51) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.42e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATA > REL606/154427‑154571 | tGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGa > 1:1085729/1‑76 (MQ=255) gAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaa < 2:1199321/76‑1 (MQ=255) gAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGaa > 1:264270/1‑76 (MQ=255) gAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGaa > 1:386107/1‑76 (MQ=255) aGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaa < 1:821505/76‑1 (MQ=255) aGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaa < 1:1195054/76‑1 (MQ=255) gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaaa < 1:973384/76‑1 (MQ=255) gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaaa < 1:1081544/76‑1 (MQ=255) gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaaa < 1:442344/76‑1 (MQ=255) gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaaa < 1:1756128/76‑1 (MQ=255) gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaaa < 1:1278584/76‑1 (MQ=255) gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGaaaa > 1:835889/1‑76 (MQ=255) tGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAt > 2:1151284/1‑76 (MQ=255) ggATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAAttt < 1:1797900/76‑1 (MQ=255) ggATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGGAAAttt > 1:1777307/1‑76 (MQ=255) gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTa < 1:1455013/76‑1 (MQ=255) aTGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAc < 1:1208551/76‑1 (MQ=255) aTGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAc < 1:1208558/76‑1 (MQ=255) cATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCg < 1:277567/76‑1 (MQ=255) atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTAccccc > 2:1395516/1‑74 (MQ=255) tatataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCCCa > 1:829353/1‑76 (MQ=255) atataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAg < 2:893715/76‑1 (MQ=255) ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTg < 1:353248/76‑1 (MQ=255) taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTAACCCCAGGGa > 2:404860/1‑76 (MQ=255) aGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGAc < 1:1272265/76‑1 (MQ=255) aGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGGAAATTTACCCCCAGGGAc > 1:266533/1‑76 (MQ=255) tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACNtt < 2:1337539/76‑1 (MQ=255) tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt < 1:671096/76‑1 (MQ=255) ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTTACCCCCGGGACATTg > 2:702726/1‑76 (MQ=255) ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGTGGCATTg > 1:1669636/1‑76 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTTCtt > 2:636297/1‑76 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCtt < 2:1114974/76‑1 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTAACCCCAGTGACCTTGCtt > 2:636316/1‑76 (MQ=255) aTCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTTCtttc > 2:1247737/1‑74 (MQ=255) tGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGt < 1:324583/76‑1 (MQ=255) tGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGGGACATTGCTTACGGt > 2:1510094/1‑76 (MQ=255) gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtt > 1:1789884/1‑76 (MQ=255) gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtt < 2:1065455/76‑1 (MQ=255) tCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGttt < 2:862013/76‑1 (MQ=255) cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGTAGAAAATTTATCCGCAGTGACATTTTTTATGGgttt > 1:541453/1‑76 (MQ=255) cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtttt < 2:279892/76‑1 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCTGTGACATTGCTTACGGTTTTa < 2:1654408/76‑1 (MQ=255) tCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGGGATTGCTTTCGGTtttttt > 2:1428511/1‑76 (MQ=255) cTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTACCCGCAGTGACATTTCTTTCGGGTtttttt > 1:1571012/1‑76 (MQ=255) cTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAACTTTACCCGCAGTGACATTGCTTACGGtttttttt > 1:1650364/1‑76 (MQ=255) cTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGGTtttttt > 1:763197/1‑76 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt < 1:1590508/76‑1 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt < 1:1630619/76‑1 (MQ=255) ttagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttttt > 2:1846609/1‑76 (MQ=255) ttagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGTGACATTGCTTACGGTtttttttttt > 1:317185/1‑76 (MQ=255) tagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttttttt < 2:1658461/76‑1 (MQ=255) agttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTc < 1:1383665/76‑1 (MQ=255) gttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcg < 2:910873/76‑1 (MQ=255) ttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcga < 1:1360105/76‑1 (MQ=255) ttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTTCTTTCGGGTtttttttttttta > 2:1173051/1‑73 (MQ=255) tataGGTTACTGTAAATGTGGCAGGGGCTGAAAAATTACCCGCAGTGGCACTGCTTACGGGTTTATTTTTTTcgaa > 1:1725305/1‑75 (MQ=255) gTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCa < 2:410459/76‑1 (MQ=255) gTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCa < 2:526092/76‑1 (MQ=255) gTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCa < 2:1774538/76‑1 (MQ=255) ttACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGGTTTTTTTTTTTTGACGACCAg > 2:1160181/1‑76 (MQ=255) tACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGGTTTATTTTTTTCGACGACCAgg > 1:703892/1‑75 (MQ=255) cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg < 2:1540167/76‑1 (MQ=255) cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg < 2:158413/76‑1 (MQ=255) cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg < 1:1567010/76‑1 (MQ=255) gTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAg < 1:1093556/76‑1 (MQ=255) tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc < 2:450054/76‑1 (MQ=255) tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc < 1:11414/76‑1 (MQ=255) tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc < 2:1527834/76‑1 (MQ=255) tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc < 2:1385357/76‑1 (MQ=255) tAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTTGACGACCAGTCGAGc > 2:1062633/1‑76 (MQ=255) aTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCAca < 2:484621/76‑2 (MQ=255) tgtgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTc < 2:112251/76‑1 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt < 2:1565655/76‑1 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt > 1:1471898/1‑76 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt > 1:281363/1‑76 (MQ=255) ggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATc < 1:1563790/76‑1 (MQ=255) ggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATc < 1:1289913/76‑1 (MQ=255) ggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATc < 1:726386/76‑1 (MQ=255) ggCAGGTGCAGAAAAGGTACCCGAAGTGACATTGCGGACGGTTTTAGTTTGTTCGACGACCAGGCGAGCACTCATc < 2:1521657/76‑1 (MQ=255) cAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCgg < 1:525364/76‑1 (MQ=255) ggTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGat < 2:924671/76‑1 (MQ=255) gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata < 1:795192/76‑1 (MQ=255) | TGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATA > REL606/154427‑154571 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |