Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,141,009 | 0 | G | A | 21.3% | 159.8 / 12.9 | 61 | E28K (GAA→AAA) | yijC | DNA‑binding transcriptional repressor |
Reads supporting (aligned to +/- strand): ref base G (39/9); new base A (0/13); total (39/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.58e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTG > REL606/4140934‑4141080 | aTTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg < 1:500665/76‑1 (MQ=255) tctcTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaa < 1:1873094/76‑1 (MQ=255) ctctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaa > 1:74851/1‑76 (MQ=255) ctctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaa < 2:1375619/76‑1 (MQ=255) tctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaa > 1:396727/1‑76 (MQ=255) ctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaa > 1:1873614/1‑76 (MQ=255) ctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaa > 2:972947/1‑76 (MQ=255) tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 2:1035995/76‑1 (MQ=255) tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 1:632221/76‑1 (MQ=255) ggTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc < 1:624480/76‑1 (MQ=255) tataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAccc > 2:1145764/1‑76 (MQ=255) ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg > 1:884329/1‑76 (MQ=255) ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg > 1:993264/1‑76 (MQ=255) ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg > 1:249869/1‑76 (MQ=255) ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg > 1:274781/1‑76 (MQ=255) ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg > 2:689998/1‑76 (MQ=255) tGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGt > 1:387258/1‑76 (MQ=255) tGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATAGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGt > 1:470478/1‑76 (MQ=255) gcagTTTTTTGGGAAAGAATTAAAACATGGTGAGGGGCGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCt < 2:1089449/76‑1 (MQ=255) cagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtg > 2:975216/1‑76 (MQ=255) agTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtgg > 1:1219951/1‑76 (MQ=255) agTACTTTTGGAAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtgg < 2:301434/76‑1 (MQ=255) aCTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACACAAAGAAAAAACCAGCCGTTCGCgggtgg > 2:1465305/1‑76 (MQ=255) aCTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtgg > 1:10224/1‑76 (MQ=255) ttttgGGCAAGGTTTCAGACTTCGTGAGGGGCGTAGGAGCGAAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGa < 2:1284855/75‑1 (MQ=255) cTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGa > 2:884514/1‑76 (MQ=255) cTTTTGGCAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGAAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGa < 2:200425/76‑1 (MQ=255) cTTTGGGCAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGa < 2:1418350/76‑1 (MQ=255) ttttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGaa > 2:548950/1‑76 (MQ=255) tttgGGCAAGGATTCAGACATGGGGATGGGGGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGaa < 1:1506610/76‑1 (MQ=255) tttgGGAAAGGATCAAAACATCGTGAGGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGaa < 1:1708887/76‑1 (MQ=255) ttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc > 1:1172030/1‑76 (MQ=255) tGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGcc > 1:156430/1‑76 (MQ=255) tGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGcc > 1:93232/1‑76 (MQ=255) ggaaaGGATCAAGACATGGTGAGGGGCAAAAGACCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCg < 2:693392/73‑1 (MQ=255) ggCAAGGATTCAGACATGGTGAGGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCg < 2:74381/76‑1 (MQ=255) gATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTa > 2:556237/1‑76 (MQ=255) aTTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAg > 2:1453360/1‑76 (MQ=255) ttCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc > 1:458472/1‑76 (MQ=255) cAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa > 1:1408196/1‑76 (MQ=255) aGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa > 2:621042/1‑76 (MQ=255) aGACATCGTGAGGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa < 1:685739/76‑1 (MQ=255) aGACAGCGGGATGGCCGAAAGAGCGCAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa < 2:1851384/76‑1 (MQ=255) tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGc > 2:429580/1‑76 (MQ=255) gaggggCGTAAGAGCGCAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCt < 2:122848/73‑1 (MQ=255) gaggggCGAAAGAGCACAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCt < 2:7632/73‑1 (MQ=255) aTGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGATg > 1:603596/1‑76 (MQ=255) tGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGa > 2:543450/1‑76 (MQ=255) gCGAAAGACCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACg < 1:574424/76‑1 (MQ=255) cgaaaGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc < 2:1193036/73‑1 (MQ=255) tAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAg > 2:254166/1‑76 (MQ=255) tAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAg > 1:1477544/1‑76 (MQ=255) aaaGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAg < 2:548845/75‑1 (MQ=255) aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGc > 2:1800251/1‑76 (MQ=255) aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGc > 2:710884/1‑76 (MQ=255) aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGc > 1:962220/1‑76 (MQ=255) agagCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCt > 2:1625283/1‑76 (MQ=255) gagCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt < 2:1709610/76‑1 (MQ=255) aagCACAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt < 2:1138343/75‑1 (MQ=255) agCACAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc < 2:252447/76‑1 (MQ=255) caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCa > 2:1555701/1‑76 (MQ=255) caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCa > 1:90361/1‑76 (MQ=255) aacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAg > 1:318215/1‑76 (MQ=255) aacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAg > 2:1555371/1‑76 (MQ=255) caaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGcc > 2:823360/1‑76 (MQ=255) aaaaaaaaCACCCCGCCGTTCGGTGTTGGAAGCCGCATTTAGCCAATTAAGTGCTGACCGCAGCTTCGCCAGACTg < 2:1392442/76‑1 (MQ=255) | ATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTG > REL606/4140934‑4141080 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |