Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,280,704 | 0 | A | C | 62.1% | 3.5 / 27.3 | 29 | T25P (ACC→CCC) | radD | putative DNA repair helicase RadD |
Reads supporting (aligned to +/- strand): ref base A (10/1); new base C (0/18); total (10/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.49e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACG > NC_000913/2280640‑2280768 | tACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATaaaa > 1:1360153/1‑65 (MQ=255) cccTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAcc < 1:3062641/63‑1 (MQ=255) acTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAccc > 1:1138980/1‑65 (MQ=255) ccccATATCAGCAAGAGGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGcc < 1:2532382/64‑1 (MQ=255) cccATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCg > 1:1852610/1‑65 (MQ=255) ccAAATCAGAAAAAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGt < 1:2435765/65‑1 (MQ=255) tatCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTAt > 1:2620133/1‑65 (MQ=255) cAGAAAGAAGCGGGGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGt < 1:2671036/65‑1 (MQ=255) cAGAAAAAAGCCGTGGATGCCCCGCTCACCCATTTCCGTCGTCATAAACCCCCTGCCGTTATCGt < 1:2635877/65‑1 (MQ=255) aaGCAAAAAGCCGTGGATGCCACGCCCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGt < 1:3767906/64‑1 (MQ=255) cAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgct < 1:2726054/65‑1 (MQ=255) aaagaggCCGTGGATGCCACGCTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgct < 1:2914152/64‑1 (MQ=255) aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTCTCgtggtg > 1:3063412/1‑65 (MQ=255) aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgctg > 1:592115/1‑65 (MQ=255) aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgctg > 1:3707277/1‑65 (MQ=255) ggaagCCGTGGATGCCACGCTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc < 1:2650078/64‑1 (MQ=255) agaagCCGTGGATCCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc < 1:2379903/65‑1 (MQ=255) aaaagCCGTGGATGCCACGCTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc < 1:3360789/63‑1 (MQ=255) tGGATGCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTNATCGTGCTGCCCACCGgc < 1:2785398/65‑1 (MQ=255) aTGCCACGCTCAACCATTTTCGTCGTTATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa > 1:3331811/1‑65 (MQ=255) aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa > 1:1136812/1‑65 (MQ=255) aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa > 1:3405780/1‑65 (MQ=255) gCCACGCTCAACCATTTTCGTCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGNAGCAgg < 1:2204918/65‑1 (MQ=255) ccACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGt < 1:2354708/65‑1 (MQ=255) aCGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTac > 1:3338281/1‑64 (MQ=255) cccAACCATTTTCGCCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg < 1:2458167/63‑1 (MQ=255) cTCACCCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg < 1:1893587/65‑1 (MQ=255) cTCAACCTTTTTCTGCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg < 1:3043073/65‑1 (MQ=255) cTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg > 1:1513104/1‑65 (MQ=255) cTCAACCATTTTCGGCGTAATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg < 1:2736584/65‑1 (MQ=255) cTCAACCATTTGCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg < 1:445961/65‑1 (MQ=255) accAACCATTTTCTTCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg < 1:2681428/63‑1 (MQ=255) cctttttCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTgg < 1:1172453/62‑1 (MQ=255) ctttttCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGt < 1:2619437/63‑1 (MQ=255) cgtcATAACCCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGa < 1:3291029/65‑1 (MQ=255) cgtcATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGa < 1:2118573/65‑1 (MQ=255) tAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTgg < 1:2264507/65‑1 (MQ=255) tAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGATGATAGCGGAACTgg < 1:3198621/65‑1 (MQ=255) aCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACg > 1:1090439/1‑65 (MQ=255) | TACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACG > NC_000913/2280640‑2280768 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |