| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,014,794 | G→C | 22.0% | A425A (GCC→GCG) | fumA ← | aerobic Class I fumarate hydratase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,014,794 | 0 | G | C | 22.0% | 110.6 / 28.6 | 59 | A425A (GCC→GCG) | fumA | aerobic Class I fumarate hydratase |
| Reads supporting (aligned to +/- strand): ref base G (31/15); new base C (8/5); total (39/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.46e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 | |||||||||||
ACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGC > minE/1014726‑1014863 | aCATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGc < 1:2207995/71‑1 (MQ=255) aTAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCg > 1:1045288/1‑71 (MQ=255) aaGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGa > 1:3023545/1‑71 (MQ=255) atccatccGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGAccc > 1:2789702/1‑70 (MQ=255) cgtcccGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCgtagta < 1:3608595/71‑1 (MQ=255) ctgcggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtag > 1:2669235/3‑65 (MQ=255) tgcggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtag > 1:3835326/2‑64 (MQ=255) ggcggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCg < 1:1963437/71‑1 (MQ=255) cggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCggg > 1:3979551/1‑70 (MQ=255) ggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgat < 1:2927941/71‑1 (MQ=255) gTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgatg < 1:1281739/71‑1 (MQ=255) gggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCttt > 1:2622623/1‑71 (MQ=255) gggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCtt > 1:82283/1‑70 (MQ=255) ggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt < 1:3307703/70‑1 (MQ=255) ggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt < 1:549856/70‑1 (MQ=255) ccAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTaaa < 1:830660/55‑1 (MQ=255) ccAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTaaa < 1:1415740/55‑1 (MQ=255) cAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATg < 1:5629680/71‑1 (MQ=255) aaGTGAACCTGATGGATAACCGGCAGGTGTTTTCGCCGGACCcgc < 1:1037267/45‑1 (MQ=255) aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATg > 1:2967144/1‑70 (MQ=255) aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATg > 1:2722137/1‑70 (MQ=255) aaGAGAACGGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATc < 1:961661/63‑1 (MQ=255) gtgaACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtag < 1:1254013/45‑1 (MQ=255) gtgaACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgta > 1:448847/3‑46 (MQ=255) gagaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGga > 1:4586304/1‑56 (MQ=255) gagaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTAc > 1:47173/1‑71 (MQ=255) agaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt < 1:770335/63‑1 (MQ=255) gaACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTg > 1:5638293/1‑71 (MQ=255) aaCCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTg < 1:4151197/62‑1 (MQ=255) ggAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa > 1:590103/1‑59 (MQ=255) ggAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa > 1:86848/1‑59 (MQ=255) gAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAg > 1:808093/1‑71 (MQ=255) aTAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTaa > 1:3987559/1‑37 (MQ=255) aaCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgatga > 1:314085/1‑44 (MQ=255) aaCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCg > 1:3455482/1‑60 (MQ=255) aaCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCg > 1:5326622/1‑60 (MQ=255) aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa > 1:1864890/1‑50 (MQ=255) aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCa > 1:3814459/1‑70 (MQ=255) aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCa > 1:5414016/1‑70 (MQ=255) aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAc > 1:101591/1‑71 (MQ=255) ccTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAcc > 1:4527079/1‑71 (MQ=255) tCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCtt > 1:549322/1‑43 (MQ=255) tCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAcc > 1:2011655/1‑68 (MQ=255) tCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt > 1:3026398/1‑70 (MQ=255) cGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGccc < 1:2412323/60‑1 (MQ=255) ggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATc < 1:2039985/35‑1 (MQ=25) ttttGGCCGGACGCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTatccatc > 1:4626157/1‑71 (MQ=255) ttttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTg < 1:235991/38‑1 (MQ=255) ttttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt < 1:1013865/63‑1 (MQ=255) ttttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTatccatc > 1:3414559/1‑71 (MQ=255) ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCg < 1:2721111/48‑1 (MQ=37) ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGAt > 1:1479560/1‑38 (MQ=255) ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGAt > 1:3959589/1‑38 (MQ=255) ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTg > 1:2857754/1‑44 (MQ=255) ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTatcca > 1:4590974/1‑67 (MQ=255) ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg > 1:2099960/1‑71 (MQ=255) ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg > 1:3490139/1‑71 (MQ=255) ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg > 1:1742759/1‑71 (MQ=255) ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg > 1:1345685/1‑71 (MQ=255) ggCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGc < 1:5612886/70‑1 (MQ=255) | ACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGC > minE/1014726‑1014863 |
| 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: ‑ |