| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,118,209 | 0 | G | T | 39.1% | 33.1 / 39.1 | 46 | intergenic (‑444/‑311) | csgD/csgB | DNA‑binding transcriptional activator in two‑component regulatory system/curlin nucleator protein, minor subunit in curli complex |
| Reads supporting (aligned to +/- strand): ref base G (4/24); new base T (17/1); total (21/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.35e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TAAAACGTTCTTGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACG > REL606/1118134‑1118284 | tAAAACGTTCTTGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTtg < 2:61323/76‑1 (MQ=255) aaaaCGTTCTTGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTtgt < 2:3782549/76‑1 (MQ=255) aaaCGTTTTTTTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTtgtg < 2:619714/76‑1 (MQ=255) ttCTTGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCaa < 1:2369879/76‑1 (MQ=255) tCTTGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTTTGCAA‑g > 2:1032516/1‑76 (MQ=255) cTTGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑Ga < 2:2758845/76‑1 (MQ=255) ttGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GAc < 2:4191485/76‑1 (MQ=255) ttGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GAc < 1:1707671/76‑1 (MQ=255) tGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GACa < 1:2481171/76‑1 (MQ=255) ttttttCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGCGCGA‑GACata > 2:3343172/1‑76 (MQ=255) tttttCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGCGGCA‑GACatat > 2:3316659/1‑76 (MQ=255) tttttCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GACatat < 1:4785647/76‑1 (MQ=255) tttCTCCACCCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGAAA‑GACATATCg < 2:2985673/76‑1 (MQ=255) tttCTCCACCCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGAAA‑GACATATCg < 2:2182928/76‑1 (MQ=255) ttCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGGGCCA‑GAAAAATCgc > 2:39861/1‑76 (MQ=255) tctcCACCCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GACATATCgca < 1:1768862/76‑1 (MQ=255) ctcCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTTTGGCA‑GACATATTgcag > 1:3042573/1‑76 (MQ=255) tcCCCACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GACATATCgcagc < 1:751922/76‑1 (MQ=255) tcCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGCGCAA‑GACATATCgcaac > 2:2565807/1‑74 (MQ=255) ccACCCCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGAAA‑GACATATCgcagca < 2:789986/76‑1 (MQ=255) ccACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTTCGCAA‑GACATATAccagca > 2:2969150/1‑76 (MQ=255) ccACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGTGCAA‑GACATATCgccgca > 2:2170732/1‑76 (MQ=255) ccccTCCGTGGACAATTTTTTACTGCAAAAAGACGAGTTTT‑GTCACGGCTTGTGCGCAA‑GACATATCGCAGCAATc < 2:2861021/74‑1 (MQ=255) ccccTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGAAA‑GACATATCGCAGCAATc < 2:2974153/74‑1 (MQ=255) cacCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTTTGCAA‑GACCTATAGCCGCaaaa > 2:3707448/1‑74 (MQ=255) cacCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTTCGCGC‑AACATATCACAGTaaac > 1:3570820/1‑74 (MQ=255) cccTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GACATATCGCAGCAATCa < 1:309107/75‑1 (MQ=255) ccTGCGTGGACAATTTTTTACGGCAAAAAGACGAGGTTT‑GTCGCGGCTTGTGCTCAA‑GACATATCGCAGCAATCAg < 2:1723733/76‑1 (MQ=255) ccTCCGTGGACATTTTTTTACTGCAAAAAGACGAGGTTT‑GTCAGGGCTTGTGCAAAA‑GACATATCGCAGCAATCAg < 2:2109210/76‑1 (MQ=255) ccTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGAAA‑GACATATCGCAGCAATCAg < 1:1389454/76‑1 (MQ=255) cTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGGGCAA‑GGCCAATCGCAGCAACCAGc > 2:4742728/1‑76 (MQ=255) cTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGAAA‑GACATATCGCAGCAATCAGc < 1:2711065/76‑1 (MQ=255) tCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTTTGGCA‑AACATATCGCAGCAATCACCg > 1:444285/1‑76 (MQ=255) tCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGGGCAA‑GACAAATCGCAGCAAACAGCg > 1:5168478/1‑76 (MQ=255) tCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GACATATCGCAGCAATCAGCg < 2:1454548/76‑1 (MQ=255) cGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GACATATCGCAGCAATCAGCGAc < 1:1838007/76‑1 (MQ=255) tGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGGGGAA‑GACAAATCGCAGCAAACAGAGACgg > 2:3713365/1‑76 (MQ=255) tGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGGGCAA‑GACATATCGCAGCAAACAGCGAccg > 1:2557647/1‑74 (MQ=255) gggACAATTTTTTACTGCAAAAAGACGAGTTTT‑GTCACGGCTTTTGCGCAA‑GACATATCGCAGCAATCAGCGACgg < 2:72792/75‑1 (MQ=255) ggACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCCCGCCTTGTGCGAAA‑GACATATCGCAGCAATCAGCGACggg < 2:3546458/76‑1 (MQ=255) aaaaTTTTTTACTGCAAAAAGACGAGTTTT‑GTCGCGGCTTGTGCAAAA‑GACATATCGCAGCAATCAGCGACGGGCa < 2:4507838/74‑1 (MQ=255) aCAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGGGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCa < 2:3362557/76‑1 (MQ=255) tttttACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCCAAA‑GACATATCGCAGCAATCAGCGACGGGCaagaag < 1:2170732/76‑1 (MQ=255) tttttACTAAAAAAAGACGAGTTTT‑GTCACGGCTTGTGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCaagaag < 2:3129574/76‑1 (MQ=255) ttttACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTTCGGAAGGACAAAT‑GCAGGCATCAGCGACGGGCAagaaga > 2:4777159/1‑76 (MQ=255) tttACTGCAAAAAGACGAGGTTT‑GTCACGTCTTGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAgaagaa < 1:2472308/76‑1 (MQ=255) tttACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGGGCAA‑GAAATATCGCAGCAAACAGCGAAGGGCAAgaagaa > 2:4797370/1‑76 (MQ=255) tttACGGCAAAAAGACGAGGTTT‑GCCCCGGCTTGGGCAAAA‑GACATATCGCAGCAATCAGCGACGGGCAAgaagaa < 1:3352202/76‑1 (MQ=255) ttACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAAt < 2:510692/76‑1 (MQ=255) tACTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGCGCAA‑GACATATCGCAGCAAACAGCGACGGGGCAgaagaaag > 1:1072014/1‑74 (MQ=255) aCTGCAAAAAGACGAGGTTT‑GTCACGGCTTTTGTGCCA‑GACATATCGCAGCAAATAGCGACGGGCAAGCAGAATGa > 1:1211868/1‑76 (MQ=255) aCTGCAAAAAGACGAGGTTT‑GTCACGGCTTTGGTCGAA‑GACATATCGCAGCAATCAGAGCCGGGGAAAAAAAAAGa > 1:2098923/1‑76 (MQ=255) cTGAAAAAAGACGAGGTTT‑GTCACGTCTTGGGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGAc < 1:1402315/76‑1 (MQ=255) gaaaaaaGACGGGGTTT‑GTCACGGCTTGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTg < 2:219612/74‑1 (MQ=255) gaaaaaaGACGAGTTTT‑GTCACGGCTGGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTg < 2:4261651/74‑1 (MQ=255) gaaaaaaGACGAGGTTT‑GTCACGGCTTGTGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTg < 2:4775543/74‑1 (MQ=255) gCAAAAAGACGAGGTTT‑GTCACGGCTTTTGTGCCA‑GACATATCGCAGCAAACAGCGACGGGGAAGAAAAAAGAATg > 2:2885492/1‑76 (MQ=255) gCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGGAA‑AACAAATCGCAGCAATCAACGACGGGCAAGAAAAAAGAATg > 2:2140078/1‑76 (MQ=255) cAAAAAGACGAGTTTT‑GTCACGGCTTGTGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGt < 2:4123036/76‑1 (MQ=255) cAAAAAGACGAGGTTT‑GTCACGGCTTTTGGGCAA‑GACATATCGCCGGCATCAGCGGCGGGGAAGAAGAATGACTGt > 1:4389219/1‑76 (MQ=255) cAAAAAGACGAGGTTT‑GTCACGGCTTTTGCGCAA‑GAAATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGt > 1:2567103/1‑76 (MQ=255) aaaaaaCACGAGGTTT‑GCCGCGGCTTGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGt < 1:5153689/75‑1 (MQ=255) aaaaaGACGAGGTTT‑GTCACGGCTTTTTCGCAA‑GACATATCTCAGCAAACAACGACGGGGAAGAAAAAAGACTGTc > 2:345323/1‑76 (MQ=255) aaaaGACGAGGTTT‑TGCACGGCTTGTGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCt < 1:1305910/76‑1 (MQ=255) aaaaGACGAGGTTT‑GTCACGGCTTGTGCACAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCt < 2:2488664/76‑1 (MQ=255) aaaGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTg < 1:3929833/76‑1 (MQ=255) aaaGACGAGGTTT‑GTCACGGCTTGTGCCAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTg < 2:5079215/76‑1 (MQ=255) aaGACGAGGTTT‑GTCACGGCTTTTGCGCCA‑GACAAATCGCCGCAATCAGCGACCGGGAAGAAGAATGACTGTCttg > 2:656294/1‑74 (MQ=255) aCGAGGTTT‑GTCACGGCTTTTGGGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTgg > 2:2532736/1‑75 (MQ=255) ggggTTT‑GTCACGGCTGGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCtt < 2:599492/74‑1 (MQ=255) ggTTTGGTCACGGCTTGTGCACAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCttt < 2:3644240/76‑1 (MQ=255) ttt‑GTCACGGCTTGTGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTg < 1:1875627/76‑1 (MQ=255) t‑GTCACGGCTTTTGCGCAA‑GACAAATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTTGGGGtttttttt > 2:509763/1‑73 (MQ=255) ttCACGGTTGGGGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATa < 2:4300762/75‑1 (MQ=255) tcccGGCTGGGGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAg < 2:703725/73‑1 (MQ=255) ccACGGCTTGTGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAg < 1:4170068/75‑1 (MQ=255) cACGGCTTGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGGGGTTTTTGATAga > 1:1290070/1‑75 (MQ=255) cACGGCTTGTGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGc < 1:5117424/76‑1 (MQ=255) cgtcTTGTGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCgg < 2:737233/73‑1 (MQ=255) cGGCTTTTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTTATAGCgg > 2:2344201/1‑76 (MQ=255) cGGCTTGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCgg < 2:1218680/76‑1 (MQ=255) ggCTTGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGa > 2:4335554/1‑76 (MQ=255) gCTTGTGCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGaa < 2:4893601/76‑1 (MQ=255) gCTTGTGCAAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGaa < 1:2603786/76‑1 (MQ=255) ttGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGaaaa > 1:3559501/1‑76 (MQ=255) ttGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGaaaa > 1:3575192/1‑76 (MQ=255) ttGTGCAAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGaaaa < 2:2625323/76‑1 (MQ=255) tgtgCGAAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAAc < 1:1906186/76‑1 (MQ=255) gtgCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACg > 2:4654186/1‑76 (MQ=255) gtgCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACg > 2:4425033/1‑76 (MQ=255) gtgCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACg > 2:1116776/1‑76 (MQ=255) | TAAAACGTTCTTGTTTTTTCTCCACACCTCCGTGGACAATTTTTTACTGCAAAAAGACGAGGTTT‑GTCACGGCTTGTGCGCAA‑GACATATCGCAGCAATCAGCGACGGGCAAGAAGAATGACTGTCTGGTGCTTTTTGATAGCGGAAAACG > REL606/1118134‑1118284 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |