Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,141,009 | 0 | G | A | 27.1% | 117.9 / 16.8 | 48 | E28K (GAA→AAA) | yijC | DNA‑binding transcriptional repressor |
Reads supporting (aligned to +/- strand): ref base G (25/10); new base A (0/13); total (25/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.93e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAG > REL606/4140934‑4141082 | aTTCTCTGGTATAGTGGCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg < 2:1270714/76‑1 (MQ=255) aTTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg > 2:430642/1‑76 (MQ=255) ttCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa > 1:1069154/1‑76 (MQ=255) ttCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa > 1:85519/1‑76 (MQ=255) tctcTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaa > 1:116809/1‑76 (MQ=255) ctctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaa < 2:633217/76‑1 (MQ=255) tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 1:529603/76‑1 (MQ=255) tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 1:1215167/76‑1 (MQ=255) tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 2:2110028/76‑1 (MQ=255) tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 2:1722281/76‑1 (MQ=255) ggTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc < 1:2015835/76‑1 (MQ=255) gTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAcc < 1:1893623/76‑1 (MQ=255) tataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAccc < 1:481635/76‑1 (MQ=255) ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg < 1:1138705/76‑1 (MQ=255) taGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgc < 2:2150370/76‑1 (MQ=255) tatttttGGAAAGGATCCAAACATGGTGAGGGGCATAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCtggtg < 2:2174254/73‑1 (MQ=255) aCTTTTGGAAAGGATTCAGACATCGTGATGGGCGAAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCtggtgg < 2:1335328/76‑1 (MQ=255) tttttGGCAAGGTTTCAGACATGGTGAGGGGCGAAAGAGCCCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGa < 2:1650668/75‑1 (MQ=255) cTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGa > 1:1067146/1‑76 (MQ=255) ttttGGCAAGGATACAGACATGGTGAGGGGCGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGaa < 2:1218502/76‑1 (MQ=255) tttgGGAAAGGATTCAGACATGGGGAGGGGCAAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGaa < 2:1550921/76‑1 (MQ=255) tttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAg > 1:701738/1‑76 (MQ=255) tttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAg > 1:1195112/1‑76 (MQ=255) ttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc > 1:245809/1‑76 (MQ=255) ggCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCg > 1:2026647/1‑76 (MQ=255) cAAGGATTCAGACATCGTGATGGGCATAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCa > 1:1589417/1‑76 (MQ=255) aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt > 1:1164925/1‑76 (MQ=255) aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt > 1:1953308/1‑76 (MQ=255) aaGGATACAGACATCGTGAGGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCAt < 2:972608/76‑1 (MQ=255) atacAGACATCGTGGGGGGCAAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAg < 2:570740/73‑1 (MQ=255) aTTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAg > 2:1524665/1‑76 (MQ=255) aTTCAAACATCGGGAGGGGCAAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAg < 1:2166764/76‑1 (MQ=255) cAGACATGGTAGGGGGCGAAAGCGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa < 1:795397/76‑1 (MQ=255) cAGACATCGTGAGGGGCGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa < 1:278333/76‑1 (MQ=255) aaCATCGTAAGGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAt < 2:1340863/75‑1 (MQ=255) aCATGGTGAGGGGGGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAtt < 1:1572042/76‑1 (MQ=255) tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGc > 2:1224464/1‑76 (MQ=255) ggCGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAAc < 2:1173160/76‑1 (MQ=255) gTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCa > 1:67741/1‑76 (MQ=255) tAAGAGCGCAACAAAAAGAGAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAACAGCTGAACGCAg > 2:1104342/1‑76 (MQ=255) tAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAg > 2:1647672/1‑76 (MQ=255) tAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAg > 1:1014458/1‑76 (MQ=255) aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGACCGCAGc > 2:1687706/1‑76 (MQ=255) gagCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt < 1:1163817/76‑1 (MQ=255) gagCACAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt < 2:157181/76‑1 (MQ=255) gcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGcc > 1:760273/1‑76 (MQ=255) caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCa > 1:406802/1‑76 (MQ=255) caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTACGCCa > 1:1328680/1‑76 (MQ=255) caaaaaAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCa < 1:2191524/76‑1 (MQ=255) aacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAg > 2:1364143/1‑76 (MQ=255) aacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAg > 2:1145482/1‑76 (MQ=255) aaaaaaaGAAAACCCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGc < 1:272397/74‑1 (MQ=255) caaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGcc > 1:679611/1‑76 (MQ=255) aaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAg > 2:1532596/1‑76 (MQ=255) aaGAAAAAACCCGCCGTTCGCTGGCGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAg > 1:1653233/1‑76 (MQ=255) | ATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAG > REL606/4140934‑4141082 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |