| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,235,108 | A→T | 44.8% | Q85L (CAG→CTG) | purH → | fused IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,235,108 | 0 | A | T | 44.8% | 14.2 / 61.7 | 58 | Q85L (CAG→CTG) | purH | fused IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (12/20); new base T (18/8); total (30/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.98e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATGGATGGACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGATCCAGCCTATCGATATGGTGGTTGTTAACCTGTAT > minE/2235044‑2235172 | atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCTGGAc < 1:808927/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCTGGAc < 1:114464/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCTGGAc < 1:406733/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCTGGAc < 1:2268275/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCTGGAc < 1:1893201/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCTGGAc < 1:1849342/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCTGGAc < 1:1676202/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGAc < 1:243928/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGAc < 1:3227901/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGAc < 1:2815036/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGAc < 1:238960/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGAc < 1:1720368/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGAc < 1:1645091/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGAc < 1:1361213/69‑1 (MQ=255)atggatggACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGCCCGTCGCGGCCTGGAc < 1:733427/69‑1 (MQ=255) tggatggaCGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGAc < 1:2381179/68‑1 (MQ=255) gcgTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGAc < 1:128679/59‑1 (MQ=255) tGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGAc < 1:1552723/56‑1 (MQ=255) ccTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAAc > 1:3354582/1‑68 (MQ=255) ccTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACa < 1:2719364/69‑1 (MQ=255) ccTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACa < 1:444604/69‑1 (MQ=255) tCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATg > 1:1190294/1‑47 (MQ=255) tCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGa > 1:3427521/1‑68 (MQ=255) tCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGa > 1:3319635/1‑68 (MQ=255) tCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGa > 1:212606/1‑68 (MQ=255) tCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGa > 1:2825076/1‑68 (MQ=255) cGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCAt < 1:2962169/49‑1 (MQ=255) cGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCAt < 1:2850796/49‑1 (MQ=255) cGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCAt < 1:511187/49‑1 (MQ=255) cGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCAt < 1:1656695/49‑1 (MQ=255) cGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCAt < 1:1055949/49‑1 (MQ=255) cGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCAt < 1:2056373/49‑1 (MQ=255) cGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCAt < 1:265356/49‑1 (MQ=255) tggtggCATTCTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGATCCAGCCTATCGa < 1:3139439/69‑1 (MQ=255) cTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGATCCa > 1:697318/1‑50 (MQ=255) cTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGATCCAGCCTATCGATATGGTGgtt > 1:2063869/1‑69 (MQ=255) cTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGATCCAGCCTATCGATATGGTGgtt > 1:179886/1‑69 (MQ=255) cTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGATCCAGCCTATCGATATGGTGgtt > 1:2969136/1‑69 (MQ=255) cTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGATCCAGCCTATCGATATGGTGgtt > 1:1278526/1‑69 (MQ=255) cTGGGCCGTAGCGGCCAGGACGATGCCATTATGGAAGAACATCAGATCCAGCCTATCGATATGGTGgtt > 1:1427460/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATAtgg > 1:2225153/1‑52 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:39049/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:3282680/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:1570410/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:1824457/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:526987/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:1705958/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:793778/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:1648895/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:349438/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:1855070/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:2979636/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:2069157/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:2869494/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:1053057/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:2813546/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:2750991/1‑69 (MQ=255) ggCCTGGACGATGCCATTATGGAAGAACATCAGATCCTGCCTATCGATATGGTGGTTGTTAACCTGTAt > 1:2737667/1‑69 (MQ=255) | ATGGATGGACGCGTGAAGACCCTGCATCCGAAAGTACATGGTGGCATTCTGGGCCGTCGCGGCCAGGACGATGCCATTATGGAAGAACATCAGATCCAGCCTATCGATATGGTGGTTGTTAACCTGTAT > minE/2235044‑2235172 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |