Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,231,892 | 0 | A | C | 36.4% | 99.7 / 17.2 | 55 | T25P (ACC→CCC) | yejH | predicted ATP‑dependet helicase |
Reads supporting (aligned to +/- strand): ref base A (30/5); new base C (0/20); total (30/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.05e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACCATGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTCGTG > REL606/2231817‑2231967 | aCCATGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATaaaa < 1:103446/76‑1 (MQ=255) aTGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAccc > 2:1488022/1‑76 (MQ=255) aTGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAccc < 2:1133261/76‑1 (MQ=255) gATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCt < 1:1763005/76‑1 (MQ=255) aTTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTg > 1:1493308/1‑76 (MQ=255) tttttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGc > 1:140746/1‑76 (MQ=255) tttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTTGTCGTCATAAAACCCCTGCCg > 1:1579952/1‑76 (MQ=255) tttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCg < 2:1738757/76‑1 (MQ=255) tttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCg > 1:974400/1‑76 (MQ=255) tttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCg > 1:252819/1‑76 (MQ=255) ttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGt < 1:1081804/76‑1 (MQ=255) tACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGtt > 2:979972/1‑76 (MQ=255) tACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGtt > 1:1976328/1‑76 (MQ=255) acacTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTa < 2:2045248/76‑1 (MQ=255) cccATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCAtgctg < 1:960267/76‑1 (MQ=255) ccATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgctgc > 2:675041/1‑76 (MQ=255) cATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGcc > 2:224375/1‑76 (MQ=255) atatCAGCAAGAAGCCGTGGATGCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGccc < 2:986346/76‑1 (MQ=255) tatCAGCAAGAAGCCGTGGATGCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCa < 2:226202/76‑1 (MQ=255) tatCAGCAAGAAGCCGTGGATGCCCCGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCa < 1:125309/76‑1 (MQ=255) tatCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCa > 2:1424492/1‑76 (MQ=255) tatCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCa > 1:433690/1‑76 (MQ=255) tatCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCa > 1:2132863/1‑76 (MQ=255) tCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTTGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCAcc > 1:1606331/1‑76 (MQ=255) ccAGCAAGAAGCCGTGGATGCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCAcc < 1:1852037/75‑1 (MQ=255) aagaagCCGTGGATGCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGgcgc < 1:1312692/76‑1 (MQ=255) aagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAgg > 1:1898478/1‑76 (MQ=255) agCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGt < 2:1055732/76‑1 (MQ=255) agCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGt < 1:589232/76‑1 (MQ=255) gCCGTGGATGCCACGCTCAACCATTTCCGTCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTa < 2:554298/76‑1 (MQ=255) ccGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaa < 2:1931664/76‑1 (MQ=255) cGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaaa > 2:1349426/1‑76 (MQ=255) cGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGCaaa > 2:669276/1‑76 (MQ=255) ggggATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaaaa < 2:550301/74‑1 (MQ=255) gTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTCCCGTTATCGTGCTGCCCACCGGCGCAGGTaaaa > 2:1139260/1‑76 (MQ=255) tGGATGCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg < 1:1984729/76‑1 (MQ=255) tGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg > 1:933426/1‑76 (MQ=255) tGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg > 2:582843/1‑76 (MQ=255) tGGATGCCACGCTCAACCATTTNCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg > 2:296919/1‑76 (MQ=255) tGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCGGCCCACCGGCGCAGGTAAAAGCCTg < 1:293059/76‑1 (MQ=255) ccACGCCCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGt < 1:1079362/76‑1 (MQ=255) cACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTg > 2:1691356/1‑76 (MQ=255) cACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTg > 2:1448295/1‑76 (MQ=255) cccTCAACCATTTTCGTCGCCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGAt < 2:547717/74‑1 (MQ=255) cGCTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGAt < 2:586822/76‑1 (MQ=255) cGCTCACCCATTTCCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGAt < 1:576048/76‑1 (MQ=255) cGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGAt < 1:176188/76‑1 (MQ=255) gCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATa < 1:1374704/76‑1 (MQ=255) cTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAg < 1:2021579/76‑1 (MQ=255) aaCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCgg < 2:836143/76‑1 (MQ=255) aCCTTTTCCGTCGTCAAAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGa < 1:522203/76‑1 (MQ=255) ccATTTTCGTCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGaa < 1:1008781/76‑1 (MQ=255) ccATTTTCCTCCACAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGaa < 1:1319498/76‑1 (MQ=255) aTTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACt < 1:1187076/76‑1 (MQ=255) tcgtcgTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCa > 1:131955/1‑76 (MQ=255) cgtcgtCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCAc > 2:1474130/1‑76 (MQ=255) tcgtcATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACgg < 1:49164/76‑1 (MQ=255) cgtcATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGCGATAGCGGAAATGGCACggc > 2:1779436/1‑76 (MQ=255) cgtcATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGCTAAAAGCCTGGTGATAGCGGAACTGGCACggc > 1:544567/1‑76 (MQ=255) gtcATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACggct > 2:1643420/1‑76 (MQ=255) aTAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACggctggc > 1:1309122/1‑76 (MQ=255) tAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACggctggct > 2:1697768/1‑76 (MQ=255) aaaaCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTc > 1:1876667/1‑76 (MQ=255) aaaaCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTc > 2:2084578/1‑76 (MQ=255) aaCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTCgg > 2:181228/1‑75 (MQ=255) aaCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTCGt > 2:1564684/1‑76 (MQ=255) aCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTCGTg > 2:1930436/1‑76 (MQ=255) | ACCATGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTCGTG > REL606/2231817‑2231967 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |