Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,141,009 | 0 | G | A | 26.3% | 147.5 / 10.1 | 57 | E28K (GAA→AAA) | yijC | DNA‑binding transcriptional repressor |
Reads supporting (aligned to +/- strand): ref base G (37/5); new base A (0/15); total (37/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.03e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGA > REL606/4140936‑4141081 | tctcTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaa > 1:125934/1‑76 (MQ=255) tctcTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaa > 2:1599922/1‑76 (MQ=255) tctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaa > 2:1285518/1‑76 (MQ=255) tctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaa > 1:1476858/1‑76 (MQ=255) tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 2:1486539/76‑1 (MQ=255) ggTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc > 1:374872/1‑76 (MQ=255) tataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAccc > 1:280088/1‑76 (MQ=255) tataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAccc < 1:299124/76‑1 (MQ=255) ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg < 2:1250872/76‑1 (MQ=255) ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg > 2:792340/1‑76 (MQ=255) ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg < 2:1634128/76‑1 (MQ=255) ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg > 2:1390742/1‑76 (MQ=255) aGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgcc > 1:489931/1‑76 (MQ=255) gTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgccg > 1:1360780/1‑76 (MQ=255) tGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGt > 1:1002623/1‑76 (MQ=255) ccAGCAGTTTTTTGGGCAAGGATTCAGACATGGTGAGGGGCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTc < 2:463086/76‑1 (MQ=255) cagcagTTTTTTGGGAAAGGATTCAGACATCGTGAGGGGCGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCg < 1:656068/76‑1 (MQ=255) cagcagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCg > 2:65833/1‑76 (MQ=255) cagcagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCg > 2:1233045/1‑76 (MQ=255) cagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCtg < 1:65833/76‑1 (MQ=255) agTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtgg > 2:754137/1‑76 (MQ=255) gTACTTTTGGCAAGGATTCAGACATCGTGAGGGGCGAAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggt < 2:162562/76‑1 (MQ=255) gTACTTTGGGCAAGGATTAAGACATCGTGATGGGCGTAAGAGCAAAAAAAAAAAAAAAAACCCGCCGTTCGCtggt < 2:162538/76‑1 (MQ=255) ttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc > 1:1201871/1‑76 (MQ=255) gCAAGTATTCAGACATCGTGAGGGGGGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGc < 1:39800/76‑1 (MQ=255) ccAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGc > 2:1706446/2‑76 (MQ=255) cAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCa > 1:688463/1‑76 (MQ=255) ttacAGACCGCGTGGTGGGGTTAAGAGCGCAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAg < 2:1454788/73‑1 (MQ=255) aTTCAGACGTCGTGAGGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAg < 2:273647/76‑1 (MQ=255) aTTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAg > 1:847317/1‑76 (MQ=255) aTTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAg > 1:338927/1‑76 (MQ=255) tCAAACATCGTGAGGGGCGAAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc < 2:1703976/76‑1 (MQ=255) cAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCTCATTTAGccc > 2:1439326/1‑75 (MQ=255) cAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa > 2:1070923/1‑76 (MQ=255) cAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa > 1:206215/1‑76 (MQ=255) cAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa > 2:338542/1‑76 (MQ=255) cAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa > 2:1439342/1‑76 (MQ=255) aGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa > 2:213556/1‑76 (MQ=255) aGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa > 2:1179116/1‑76 (MQ=255) gACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAt > 1:491422/1‑76 (MQ=255) cATCGTGAGGGGCGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTa < 1:520137/76‑1 (MQ=255) aaTCGAGATGGGCGAAAGAGCACAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTa < 2:1099417/75‑1 (MQ=255) tCGTGAGGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAg < 2:1299892/76‑1 (MQ=255) cGTGAGGGGCGAAAGACCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGt < 2:561591/76‑1 (MQ=255) gTGAGGGGCGTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTg < 2:1492280/76‑1 (MQ=255) ggggCGTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGa < 1:150505/75‑1 (MQ=255) aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGc > 2:919771/1‑76 (MQ=255) aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGc > 2:759232/1‑76 (MQ=255) aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGc > 2:883972/1‑76 (MQ=255) gagCGCAACAAAAAGAAAAAACCCGCCGTTCACTGGTGGAAGCCGCATTTAGCCAATTAAGTCCTGAACGCAGCtt > 2:1540594/1‑76 (MQ=255) ggcgaAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc < 1:1073351/75‑1 (MQ=255) agCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc < 1:1138285/76‑1 (MQ=255) agCCCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc < 1:1195004/76‑1 (MQ=255) gcgaAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCg < 2:215804/76‑1 (MQ=255) gcgaAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCg < 1:264989/76‑1 (MQ=255) cgcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGc > 2:729128/1‑76 (MQ=255) gcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGcc > 1:1192582/1‑76 (MQ=255) caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCa > 1:1289931/1‑76 (MQ=255) aacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAg > 2:474977/1‑76 (MQ=255) aacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAg > 1:457882/1‑76 (MQ=255) aaaaaaaaaaCCCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTTCTGAACGCAGCTTCGCCAGCCt < 2:1423303/76‑1 (MQ=255) aaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGACTGa > 1:961854/1‑76 (MQ=255) | TCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGA > REL606/4140936‑4141081 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |