Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,280,704 | 0 | A | C | 37.5% | 28.5 / 16.1 | 24 | T25P (ACC→CCC) | radD | putative DNA repair helicase RadD |
Reads supporting (aligned to +/- strand): ref base A (14/1); new base C (0/9); total (14/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.65e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACG > NC_000913/2280640‑2280768 | tACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATaaaa < 1:1911403/65‑1 (MQ=255) cacTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAcc > 1:2643827/1‑65 (MQ=255) cacTCCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAcc < 1:920669/65‑1 (MQ=255) acTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATTAAAccc > 1:4156438/1‑65 (MQ=255) aatacCAGAAAGAAGCCGTGGATGCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGtt < 1:1811382/64‑1 (MQ=255) tatCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTAt > 1:1989868/1‑65 (MQ=255) tCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCg > 1:3006659/1‑65 (MQ=255) cAGCAAGAGGCGGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGt < 1:2865516/65‑1 (MQ=255) gCAAGAGGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgc > 1:1488384/1‑65 (MQ=255) aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgctg > 1:4077089/1‑65 (MQ=255) aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgctg > 1:4492416/1‑65 (MQ=255) agaagCCGTGGATCCCCCGCTCACCCTTTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc < 1:1855735/65‑1 (MQ=255) aaaagCCGTGAATGCCACCCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctgc < 1:4695878/63‑1 (MQ=255) aagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGccc > 1:4604327/1‑65 (MQ=255) aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCCCCGGCGCa > 1:523522/1‑65 (MQ=255) aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGgctgt > 1:3407391/1‑62 (MQ=255) aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGgctgt > 1:1100033/1‑62 (MQ=255) aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa > 1:3688935/1‑65 (MQ=255) aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa > 1:3603521/1‑65 (MQ=255) gCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAgg < 1:2832946/65‑1 (MQ=255) gCCACGCTAAACCATTTTCCCCGTAATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAgg < 1:1440239/65‑1 (MQ=255) gCCACGCCCACCCATTTCCGTCGCAATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAgg < 1:750915/65‑1 (MQ=255) cGCTCAACCATTTTCGTCGTCATAAAACCCCTCCCGTTATCCTGCTGCCCACCGGCGCAGGTaaa > 1:2754847/1‑65 (MQ=255) cTCACCCATTTTCGTCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg < 1:1157881/65‑1 (MQ=255) cTCAACAATTTTCTTCGTCAGAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg < 1:2347920/65‑1 (MQ=255) cacccATTTTCGTCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGcc < 1:3601669/62‑1 (MQ=255) cAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGcc < 1:261232/65‑1 (MQ=255) ccATTTTCGCCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTgg < 1:3217362/65‑1 (MQ=255) ccATTTCCGTCTTCATAAACCCCCTGCGGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTgg < 1:3457566/65‑1 (MQ=255) ccATTTCCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTgg < 1:1159691/65‑1 (MQ=255) cATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGt < 1:34580/65‑1 (MQ=255) tttCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGAt > 1:2402465/1‑65 (MQ=255) cgtcATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGa < 1:1862579/65‑1 (MQ=255) cgccATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGa < 1:4311365/62‑1 (MQ=255) catcATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGa < 1:4308241/63‑1 (MQ=255) gccATAACACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGaa < 1:3624804/63‑1 (MQ=255) tAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTgg < 1:4261612/65‑1 (MQ=255) aaaacCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTgg < 1:3283158/64‑1 (MQ=255) aCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACg > 1:339920/1‑65 (MQ=255) | TACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACG > NC_000913/2280640‑2280768 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |