| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,974,977 | C→T | 40.5% | intergenic (+123/‑301) | metK → / → galP | methionine adenosyltransferase 1/D‑galactose transporter |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,974,977 | 0 | C | T | 40.5% | 27.7 / 42.6 | 64 | intergenic (+123/‑301) | metK/galP | methionine adenosyltransferase 1/D‑galactose transporter |
| Reads supporting (aligned to +/- strand): ref base C (27/11); new base T (13/13); total (40/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.17e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATT > minE/1974910‑1975043 | aGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCa < 1:2072529/69‑1 (MQ=255)aGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCa < 1:1379740/69‑1 (MQ=255) tCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAg < 1:2274980/69‑1 (MQ=255) cACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAg > 1:1349586/1‑68 (MQ=255) cACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAg > 1:2686639/1‑68 (MQ=255) cACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAg > 1:2549387/1‑68 (MQ=255) cACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAg < 1:2260621/68‑1 (MQ=255) cACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAg > 1:677313/1‑68 (MQ=255) cgcATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatc < 1:1407572/69‑1 (MQ=255) taTCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatc < 1:526658/66‑1 (MQ=255) cGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatc > 1:2332681/1‑63 (MQ=255) cGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatc > 1:466689/1‑63 (MQ=255) cGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatc < 1:3104281/63‑1 (MQ=255) cGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatc < 1:2245530/63‑1 (MQ=255) cGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatc < 1:3263625/63‑1 (MQ=255) cGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatc < 1:3401154/63‑1 (MQ=255) cGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatc < 1:2708563/63‑1 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaaa > 1:1983160/1‑69 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaaa > 1:1061845/1‑69 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaaa > 1:3112938/1‑69 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaaa > 1:1903750/1‑69 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaaa > 1:1664505/1‑69 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaaa > 1:2769901/1‑69 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaaa > 1:2393530/1‑69 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaaa > 1:2544919/1‑69 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaaa > 1:938324/1‑69 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaa > 1:101047/1‑68 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaa > 1:2626457/1‑68 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGaa > 1:2895079/1‑68 (MQ=255) aaaaCCGCCCGCACAATAACATCATTCTTCCTAAACACGTTTTACCGCAAATTATCATCACTACTGaaa > 1:1789957/1‑69 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:893935/69‑1 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:1060156/69‑1 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:678467/69‑1 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:3222793/69‑1 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:470576/69‑1 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:24527/69‑1 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:1607541/69‑1 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:1810658/69‑1 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:1472786/69‑1 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:2159775/69‑1 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:2614908/69‑1 (MQ=255) aTAACATCATTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:2239265/69‑1 (MQ=255) catTCTTCCTCAACACGTTTTACCGCAAATTATCATCACTACTGAAATCGATTACaccaacc < 1:2110184/62‑1 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGAt > 1:3281987/1‑51 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:3432145/1‑68 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:2965475/1‑68 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:605419/1‑68 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:1843897/1‑68 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:2000441/1‑68 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:460813/1‑68 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:2270154/1‑68 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:3330079/1‑68 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:3291119/1‑68 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:2188825/1‑68 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:2204499/1‑68 (MQ=255) atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAg > 1:856574/1‑68 (MQ=255) cACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCCACCACAACAGACAAAGATTTGt > 1:1437719/1‑68 (MQ=255) cACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGt > 1:96273/1‑68 (MQ=255) cACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGt > 1:951399/1‑68 (MQ=255) cACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGt > 1:484129/1‑68 (MQ=255) cACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGt > 1:445384/1‑68 (MQ=255) cACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGt > 1:3173152/1‑68 (MQ=255) cACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGt > 1:3025029/1‑68 (MQ=255) tCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATAtt > 1:748417/1‑68 (MQ=255) | AGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATT > minE/1974910‑1975043 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |