Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,566,491 | C→G | 28.9% | T187T (ACG→ACC) | 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,491 | 0 | C | G | 28.9% | 109.7 / 51.2 | 83 | T187T (ACG→ACC) | mtlD | mannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base C (53/6); new base G (9/15); total (62/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.56e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01 |
GCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTGCCGAAG > minE/2566432‑2566558 | gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc > 1:2665490/1‑68 (MQ=255) gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc > 1:1579838/1‑68 (MQ=255) gCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATc > 1:152780/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:1463242/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:743692/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:754062/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:1601123/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:2486076/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:1430636/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:1393854/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:91391/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:2706564/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:1216927/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:3065205/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:3125401/50‑1 (MQ=255) tAATTCCATGCCTGGGATGTTCGGCAGCGCGCCTTTGAACTGGGTTTTAt < 1:100187/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:2286244/68‑1 (MQ=255) tGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCg < 1:1580575/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:2955218/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:1259490/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:82009/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:1610260/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:34754/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:1919067/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:1914400/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:316714/1‑69 (MQ=255) aTGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACt > 1:1778927/1‑69 (MQ=255) aTGTTAGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAAt > 1:2635553/1‑41 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTccg > 1:2399553/1‑68 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTccg > 1:1189838/1‑68 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:2589782/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:3060920/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:409354/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:3202237/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:2529406/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:3256386/1‑69 (MQ=255) tCGGCAGCGCGCCTTTGAACTGGGTTTTATCAACAATCCATTCGCTGAATGTTTCTACCGTCACTTCCa > 1:181888/1‑69 (MQ=255) cAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCgg > 1:2532689/1‑69 (MQ=255) cAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCgg > 1:1662977/1‑69 (MQ=255) gcgcCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGa > 1:539373/1‑40 (MQ=255) cgcCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATc > 1:464729/1‑67 (MQ=255) cgcCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATc > 1:658245/1‑67 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc > 1:3232922/1‑42 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc > 1:3272637/1‑42 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc > 1:331176/1‑42 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc > 1:1775047/1‑42 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc > 1:3211309/1‑42 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc > 1:3205733/1‑42 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc > 1:131942/1‑42 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTc > 1:1954469/1‑42 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGAt > 1:780417/1‑55 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGt > 1:2236224/1‑62 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:500118/1‑69 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:956632/1‑69 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:325013/1‑69 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:3249748/1‑69 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:2669118/1‑69 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:2663814/1‑69 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:2194021/1‑69 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:1746584/1‑69 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:1713584/1‑69 (MQ=255) tGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:1546786/1‑69 (MQ=255) gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaag > 1:3081363/1‑69 (MQ=255) gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaag > 1:2157509/1‑69 (MQ=255) gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaag > 1:1844802/1‑69 (MQ=255) gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaag > 1:1337802/1‑69 (MQ=255) gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaag > 1:1187491/1‑69 (MQ=255) gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:2930879/1‑68 (MQ=255) gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTgccgaa > 1:1039215/1‑68 (MQ=255) gCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTACCAGCGGATCGTTAGTTgccgaag > 1:1649968/1‑69 (MQ=255) | GCCATCAGGTTGTCGGTTAACTCCATGCCTGGGATGTTCGGCAGTGCGCCTTTGAACTGCGTTTTATCGACAATCCATTCGCTGAAGGTTTCTACCGTCACTTCCAGCGGATCGTTAGTTGCCGAAG > minE/2566432‑2566558 |
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: ‑ |