Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,141,009 | 0 | G | A | 22.2% | 178.3 / 19.4 | 72 | E28K (GAA→AAA) | yijC | DNA‑binding transcriptional repressor |
Reads supporting (aligned to +/- strand): ref base G (45/11); new base A (0/16); total (45/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.17e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTT > REL606/4140934‑4141084 | aTTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg < 2:2010100/76‑1 (MQ=255) aTTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg < 2:1986629/76‑1 (MQ=255) aTTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg > 1:1866060/1‑76 (MQ=255) ttCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa > 1:1159514/1‑76 (MQ=255) ttCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa < 1:528534/76‑1 (MQ=255) ttCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa > 2:919295/1‑76 (MQ=255) ttCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa < 1:1324987/76‑1 (MQ=255) tctcTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaa > 2:845621/1‑76 (MQ=255) ctctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaa > 1:1212938/1‑76 (MQ=255) tctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaa < 1:1941135/76‑1 (MQ=255) ctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaa > 1:633766/1‑76 (MQ=255) tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 2:1895011/76‑1 (MQ=255) ggTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc > 1:78227/1‑76 (MQ=255) tataGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAccc > 1:325473/1‑76 (MQ=255) tataGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAccc < 1:792902/76‑1 (MQ=255) ataGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg < 1:818741/76‑1 (MQ=255) aGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgcc < 1:619589/76‑1 (MQ=255) gTGCCAGCAGTATTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAACCccgccg < 2:1493267/76‑1 (MQ=255) tGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGt > 1:2201479/1‑76 (MQ=255) gcagTA‑TTTTGGAAAAGAATCAAGACATCGTGAGGGGCGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCt < 1:722812/76‑1 (MQ=255) cagTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtg > 1:1956190/1‑76 (MQ=255) cagTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtg > 1:462569/1‑76 (MQ=255) cagTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtg > 1:1082377/1‑76 (MQ=255) gTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggt > 1:1991166/1‑76 (MQ=255) gTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGAAAGAGCGCAAAAAAAAGAAAAAACCCGCCGTTCGCtggt < 2:63947/76‑1 (MQ=255) tACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtg > 2:2149043/1‑76 (MQ=255) tACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtg > 2:1093976/1‑76 (MQ=255) tACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtg < 2:1920744/76‑1 (MQ=255) aCTTTTGGCAAGTATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAAAAAAACCCCGCCGTTCGCtggtgg < 2:755433/76‑1 (MQ=255) aCTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtgg > 1:714220/1‑76 (MQ=255) ttttgGGAAAGGTTGC‑AGACATCGAGATGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGa < 1:2003507/75‑1 (MQ=255) ttttGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGaa > 1:1891192/1‑76 (MQ=255) tttgGGAAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGaa < 2:170188/76‑1 (MQ=255) tttGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAg > 1:1240128/1‑76 (MQ=255) tttGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAg > 2:188268/1‑76 (MQ=255) ggaaaGGATTC‑AGACTGCGTGAGGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCg < 1:1838141/73‑1 (MQ=255) ggaaaGGATTC‑AGACATGGTGAGGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCg < 1:1755578/73‑1 (MQ=255) ggaaaGGATCC‑AGACATCGGGAGGGGCAAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCg < 2:1193247/73‑1 (MQ=255) ggCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCg > 1:140056/1‑76 (MQ=255) gCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGc > 1:1976110/1‑76 (MQ=255) cAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCa > 1:38237/1‑76 (MQ=255) aaGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt > 2:1653337/1‑76 (MQ=255) aaGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt > 1:18656/1‑76 (MQ=255) aaGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt > 1:1672156/1‑76 (MQ=255) agaaTTC‑AGGCATGGGGAGGGGCAAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCAtt < 2:1576339/73‑1 (MQ=255) gATTC‑AGACATGGGGAGGGGAGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTa < 1:1374746/76‑1 (MQ=255) gATTC‑AGACATCGTGAGGGGGGTAAGGGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTa < 1:2066812/76‑1 (MQ=255) aTTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAg > 2:1199128/1‑76 (MQ=255) tC‑AGACATCGTGATGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc < 2:890641/76‑1 (MQ=255) c‑AGACATCGTGATGGGCGTAAGNGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa > 2:299748/1‑76 (MQ=255) c‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa > 1:2030247/1‑76 (MQ=255) c‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGACa > 2:2056971/1‑76 (MQ=255) aCATGGTGAGGGGCAAAAGACCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAtt < 2:953411/76‑1 (MQ=255) aCATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAtt > 2:1127286/1‑76 (MQ=255) cATCGTGATGGGCGTAAGAGCGCAACAAAACGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTa > 2:495659/1‑76 (MQ=255) aTCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTaa > 1:2094614/1‑76 (MQ=255) gggggCGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg < 2:861778/74‑1 (MQ=255) gggCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGAGGAAGCCGCATTTAGCCAATTAAGTGCTGaa > 2:2103952/1‑76 (MQ=255) cgaaaGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc < 2:886993/73‑1 (MQ=255) cgaaaGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc < 2:2082793/73‑1 (MQ=255) cGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc > 2:1614973/1‑76 (MQ=255) cGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc > 1:2201970/1‑76 (MQ=255) cGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc > 1:2052728/1‑76 (MQ=255) gTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCa > 1:1520123/1‑76 (MQ=255) tAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAg > 1:2004018/1‑76 (MQ=255) aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGc > 1:1804354/1‑76 (MQ=255) aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGc > 2:538395/1‑76 (MQ=255) agagCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCt < 1:924160/76‑1 (MQ=255) gagCGAAACAAAGAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt < 2:1747181/76‑1 (MQ=255) gagCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt < 2:2082462/76‑1 (MQ=255) gagCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt < 2:1176758/76‑1 (MQ=255) gagCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt < 2:1507648/76‑1 (MQ=255) agCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc > 2:829791/1‑76 (MQ=255) agCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc > 1:1579442/1‑76 (MQ=255) agCGAAAAAAAAAAAAAAACCNAGCCTTTCGCTGGTGGTAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc < 2:348076/76‑1 (MQ=255) gcgcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCg > 2:5419/1‑76 (MQ=255) gcgcAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCg < 2:2191290/76‑1 (MQ=255) gcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGcc > 2:195011/1‑76 (MQ=255) aaaaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAg < 2:1524645/73‑1 (MQ=255) aaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGcc < 2:48297/75‑1 (MQ=255) aaaaaaaaaaCCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt < 1:1870827/66‑1 (MQ=255) aaaaaaaaaaCCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTg < 2:481257/76‑1 (MQ=255) aaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt > 2:1870827/1‑66 (MQ=255) aaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTg > 2:1325390/1‑76 (MQ=255) aGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGt > 2:1672234/1‑76 (MQ=255) gAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGtt > 1:916764/1‑76 (MQ=255) | ATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTC‑AGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTT > REL606/4140934‑4141084 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |