Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,566,476 | T→C | 40.7% | A192A (GCA→GCG) | mtlD ← | mannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,566,476 | 0 | T | C | 40.7% | 19.6 / 41.0 | 59 | A192A (GCA→GCG) | mtlD | mannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base T (24/11); new base C (9/15); total (33/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.19e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01 |
TTACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGG > minE/2566417‑2566541 | ttACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTg < 1:1145248/69‑1 (MQ=255) ttACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTg < 1:890001/69‑1 (MQ=255) tACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGa > 1:2911975/1‑69 (MQ=255) tACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGa < 1:1982950/69‑1 (MQ=255) tACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGa > 1:2290288/1‑69 (MQ=255) tACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGa < 1:2433804/69‑1 (MQ=255) ttCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTg < 1:2675103/69‑1 (MQ=255) gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc > 1:2665490/1‑68 (MQ=255) gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc > 1:152780/1‑68 (MQ=255) gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc > 1:1579838/1‑68 (MQ=255) aTCAGGTTGTCGGTTAACTCCACGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACa > 1:321522/1‑69 (MQ=255) gTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATcc > 1:1778188/1‑68 (MQ=255) gTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCa > 1:1792763/1‑69 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:2486076/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:3125401/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:3065205/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:2706564/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:743692/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:100187/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:754062/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:1216927/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:91391/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:1393854/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:1601123/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:1430636/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:1463242/50‑1 (MQ=255) ccTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:1438561/40‑1 (MQ=38) tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTAGCg < 1:2352130/68‑1 (MQ=255) tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg < 1:516122/68‑1 (MQ=255) tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg < 1:2566838/68‑1 (MQ=255) tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg < 1:1580575/68‑1 (MQ=255) tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg < 1:2286244/68‑1 (MQ=255) ggATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg < 1:291843/66‑1 (MQ=255) gATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:952849/35‑1 (MQ=37) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTa > 1:3283925/1‑61 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:316714/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:82009/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:34754/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:2955218/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:1076020/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:2792618/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:2893544/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:1259490/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:1610260/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:1778927/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:1914400/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:1919067/1‑69 (MQ=255) aTGTTAGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAAt > 1:2635553/1‑41 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTccg > 1:1189838/1‑68 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTccg > 1:2399553/1‑68 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:3256386/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:409354/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:2529406/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:181888/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:3202237/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:2589782/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:3060920/1‑69 (MQ=255) cAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCgg > 1:1662977/1‑69 (MQ=255) cAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCgg > 1:2532689/1‑69 (MQ=255) | TTACGTTCGACAAATGCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGG > minE/2566417‑2566541 |
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: ‑ |