Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 154,497 | 0 | T | G | 23.7% | 181.2 / 22.5 | 76 | A191A (GCA→GCC) | yadL | predicted fimbrial‑like adhesin protein |
Reads supporting (aligned to +/- strand): ref base T (8/50); new base G (18/0); total (26/50) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.25e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA > REL606/154423‑154572 | aCGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTg < 1:1050068/76‑1 (MQ=255) gTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCa > 1:1286583/1‑76 (MQ=255) tGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGa < 1:1054755/76‑1 (MQ=255) tGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGgaga > 2:678814/1‑76 (MQ=255) gAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaa < 2:1919414/76‑1 (MQ=255) gAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaa < 1:1632626/76‑1 (MQ=255) gAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGaa > 2:1968916/1‑76 (MQ=255) aGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGaaa > 2:1507033/1‑76 (MQ=255) gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaaa < 2:984871/76‑1 (MQ=255) gggATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGaaaaat > 1:1188180/1‑74 (MQ=255) gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTa < 1:1512953/76‑1 (MQ=255) aTGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTAc > 1:192133/1‑76 (MQ=255) atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGc < 1:488737/76‑1 (MQ=255) atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGc < 1:774274/76‑1 (MQ=255) ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCACAAAATTTACCCGCAGTg < 2:50683/76‑1 (MQ=255) taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGa < 2:743814/76‑1 (MQ=255) taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGa < 2:58871/76‑1 (MQ=255) taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGa < 2:518706/76‑1 (MQ=255) gTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGTGACa > 1:248843/1‑76 (MQ=255) ttCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGGCAt > 1:699184/1‑76 (MQ=255) tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt < 2:1104464/76‑1 (MQ=255) tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt < 2:479389/76‑1 (MQ=255) ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTg < 1:904213/76‑1 (MQ=255) ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTg > 2:939309/1‑76 (MQ=255) tGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTACCCGCAGTGACATTGCt > 2:1959683/1‑76 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCtt < 2:685175/76‑1 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCCGAAAAATTACCCCCCGTGACACTGCtt > 1:488093/1‑76 (MQ=255) gATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTa < 2:1342813/76‑1 (MQ=255) gATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTATCCGCAGGGACATTGtttt > 1:1072599/1‑75 (MQ=255) aTCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTAc < 1:1074681/76‑1 (MQ=255) aTCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTTCCCCCAGTGGGATTTCtttc > 2:1515098/1‑74 (MQ=255) aTCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTAc > 2:631562/1‑76 (MQ=255) aTCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTAc > 2:1712539/1‑76 (MQ=255) tCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACg < 2:232277/76‑1 (MQ=255) tCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACg < 2:232303/76‑1 (MQ=255) gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCCGAAAATTTACCCGCAGTGACATTGCTTACGGtt > 1:250302/1‑76 (MQ=255) cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTGCTTTCGGtttt > 2:912227/1‑76 (MQ=255) gTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAt < 1:1979150/76‑1 (MQ=255) tttGTTAGTTATAGG‑TACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttt < 1:480499/74‑1 (MQ=255) tCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtt < 1:1960383/76‑1 (MQ=255) cTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtttttttt > 2:1650325/1‑76 (MQ=255) cTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGTGACATTGCTTACGGtttttttt > 1:1339299/1‑76 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt < 2:772257/76‑1 (MQ=255) gttagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAGTGACATTGCTTACGGtttttttttt > 1:1356521/1‑76 (MQ=255) ttagttaTAGGTTACTGTAAATGTGGCAGGGGGAGGAAATTTTCCCGCAGGGGCCTTGCTTTTGGGtttttttttt > 1:1782518/1‑76 (MQ=255) tagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttttttt < 2:804711/76‑1 (MQ=255) tagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttttttt < 2:1712723/76‑1 (MQ=255) tagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGGTtttttttttt > 1:71196/1‑76 (MQ=255) tagttaTAGGTTACTGTAAATGTGGCAGGGGGAGAAAAATTACCCCCGATGATAATTCTTACTGGTtttttttttt > 1:780279/1‑76 (MQ=255) ttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTTTTTTTTTtga > 1:1222561/1‑76 (MQ=255) ttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcga < 2:1346630/76‑1 (MQ=255) ttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcga < 1:279710/76‑1 (MQ=255) tataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgac < 2:1250890/76‑1 (MQ=255) ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacg < 2:1911046/76‑1 (MQ=255) ataGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGTGACATTGCTTTCGGGTTTTTTTTTTTtgacg > 1:1644487/1‑76 (MQ=255) taGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGGTTTTTTTTTTTcgacga > 2:158950/1‑76 (MQ=255) gTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCa < 2:934008/76‑1 (MQ=255) gTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCa < 2:386417/76‑1 (MQ=255) ttACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAg < 1:1089621/76‑1 (MQ=255) tACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTGCTTACGGGTTTTTTTTTTTTTACGACCCGt > 1:798315/1‑76 (MQ=255) aCTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTc < 1:805238/76‑1 (MQ=255) cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg < 2:1415422/76‑1 (MQ=255) cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg < 2:548081/76‑1 (MQ=255) tGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGa < 2:1786095/76‑1 (MQ=255) tGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGa < 2:625059/76‑1 (MQ=255) tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc < 1:1481401/76‑1 (MQ=255) aaaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCa < 1:1462915/76‑1 (MQ=255) aTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACt > 1:822755/1‑76 (MQ=255) tgtgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTc < 2:289409/76‑1 (MQ=255) tgGGAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt < 2:1451926/76‑1 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt < 1:1891773/76‑1 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt > 2:1226229/1‑76 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt < 1:1220695/76‑1 (MQ=255) ggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATc > 1:487657/1‑76 (MQ=255) gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg < 1:1898797/76‑1 (MQ=255) gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg < 2:804581/76‑1 (MQ=255) gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg < 2:909735/76‑1 (MQ=255) gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg < 1:1304516/76‑1 (MQ=255) cAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCgg < 1:1708719/76‑1 (MQ=255) aGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGa > 1:1690973/1‑76 (MQ=255) aGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGa < 2:1280876/76‑1 (MQ=255) ggTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGat < 2:1644487/76‑1 (MQ=255) gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata > 2:573606/1‑76 (MQ=255) tGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGataa < 2:1113550/76‑1 (MQ=255) | ACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA > REL606/154423‑154572 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |