Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 154,497 | 0 | T | G | 35.4% | 79.2 / 16.8 | 48 | A191A (GCA→GCC) | yadL | predicted fimbrial‑like adhesin protein |
Reads supporting (aligned to +/- strand): ref base T (6/25); new base G (17/0); total (23/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.38e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA > REL606/154424‑154572 | cGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGc < 2:372700/76‑1 (MQ=255) gTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTgaa > 1:705164/1‑74 (MQ=255) ttGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAg > 1:1367590/1‑76 (MQ=255) tGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCGGa > 1:473460/1‑76 (MQ=255) aGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaa < 2:593681/76‑1 (MQ=255) taTGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTa > 1:1183636/2‑76 (MQ=255) aTGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAc > 1:1403245/1‑76 (MQ=255) aTGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAAttttc > 1:1744669/1‑74 (MQ=255) atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGc < 1:1549027/76‑1 (MQ=255) atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGCGGAGAAAAATTAccccc > 2:705912/1‑74 (MQ=255) tatataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCa < 2:361858/76‑1 (MQ=255) tatataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCa < 2:373381/76‑1 (MQ=255) tatataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCa > 2:1027254/1‑76 (MQ=255) tataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGt < 2:312984/76‑1 (MQ=255) ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑g < 1:980265/76‑1 (MQ=255) aGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCCGGAAAATTACCCGCAGT‑GAc > 1:512049/1‑76 (MQ=255) ttCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACAt < 2:961301/76‑1 (MQ=255) tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTTCCCGCAGA‑GACCtt > 2:962638/1‑76 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTgngt < 2:1254471/76‑4 (MQ=255) gATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCCCCGT‑GACATTGCTTa > 1:497023/1‑76 (MQ=255) aTCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGT‑GACATTGCTTAc > 2:1426227/1‑76 (MQ=255) cTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACgg < 2:1367590/76‑1 (MQ=255) cTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGG‑GACATTGCTTACgg > 1:882337/1‑76 (MQ=255) cTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAAATTACCCCCAGT‑GACATTGCTTACgg > 1:943115/1‑76 (MQ=255) tGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTTCGGt > 2:147312/1‑76 (MQ=255) tGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGt > 1:154579/1‑76 (MQ=255) tCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTACCCGCAGT‑GACATTGCTTTCGGttt > 2:277308/1‑76 (MQ=255) cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGG‑GACATTGCTTTCGGtttt > 1:26980/1‑76 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTa < 1:1052876/76‑1 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAAATTTCCCCCCGTGGCCATTGTTTCGGGtttt > 2:821762/1‑76 (MQ=255) gTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGAGAAAATTTACCCCCAGT‑GACATTGCTTACGGtttttt > 1:1571291/1‑74 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGT‑GACATTGCTTTCGGTtttttttt > 2:1060724/1‑76 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGGtttttttt > 2:1180921/1‑76 (MQ=255) gttagttaTAGGTTACTGTAAATGTGGCAGGGGCGGAAAATTTACCCCCCGG‑GGCATTGCTTTCGGGTtttttttt > 1:895841/1‑76 (MQ=255) ttagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTAtttttt < 2:567952/76‑1 (MQ=255) ttagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTAtttttt < 2:157635/76‑1 (MQ=255) ttagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAGT‑GACATTGCTTACGGTtttttttttt > 1:1354022/1‑76 (MQ=255) tagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTAttttttt < 2:523996/76‑1 (MQ=255) ttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTcga < 2:1436263/76‑1 (MQ=255) ttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTcga < 2:962949/76‑1 (MQ=255) ttaTAGGTTACTGTAAATGTGGCAGGGGCCGAAAATTTAACCCCGGG‑GAGATTTTTTACGGGTttttttttttttt > 1:919686/1‑73 (MQ=255) ttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGT‑GACCTTTCTTACGGtttttttttttttgc > 2:1672709/1‑73 (MQ=255) ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTcgacg < 1:523916/76‑1 (MQ=255) taGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGGTTTTTTTTTTTtgacga > 1:1262864/1‑76 (MQ=255) aGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGT‑GACATTGCTTCCGGGTTTATTTTTTTCgacgac > 2:1483172/1‑76 (MQ=255) aGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCGACC‑GACCTTGCTTACGGTTTTTTTTTTTTCgacgac > 2:1508858/1‑76 (MQ=255) ggTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGAcc < 2:799927/76‑1 (MQ=255) ggTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGAcc < 1:1570200/76‑1 (MQ=255) tACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGt < 1:1138599/76‑1 (MQ=255) aCTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTc < 1:519490/76‑1 (MQ=255) aCTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTc > 2:671193/1‑76 (MQ=255) cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg < 2:1207788/76‑1 (MQ=255) gTAAATGTGGCAGGTGCAGAAAATTTTCCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGCCCAGTCGAg < 2:824492/76‑1 (MQ=255) gTAAATGTGGCAGGGGCAGAAAATTTACCCCCCGT‑GACATTGCTTACGGTTTTTTTTTTTTTTACGACCCGTCGAg > 1:934873/1‑76 (MQ=255) aaaTGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCCGTCGAGCa > 2:724477/1‑76 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt > 2:1432093/1‑76 (MQ=255) cAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCgg < 2:1028169/76‑1 (MQ=255) aGGTGCAGAAAATTTTCCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGa < 1:1319772/76‑1 (MQ=255) ggTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGat < 2:1070560/76‑1 (MQ=255) tGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGataa > 2:358258/1‑76 (MQ=255) tGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGataa > 1:859326/1‑76 (MQ=255) | CGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGT‑GACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA > REL606/154424‑154572 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |