| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,248,116:1 | +G | 12.0% | coding (1829/1896 nt) | thiC → | thiamin (pyrimidine moiety) biosynthesis protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,248,116 | 1 | . | G | 12.0% | 223.3 / 37.3 | 75 | coding (1829/1896 nt) | thiC | thiamin (pyrimidine moiety) biosynthesis protein |
| Reads supporting (aligned to +/- strand): ref base . (38/28); new base G (7/2); total (45/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.01e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01 | |||||||||||
ATTCTGCTCGATGAAAATCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGCGTAAGGAGG > minE/2248048‑2248175 | aTTCTGCTCGATGAAAATCAGCCAGGAAGTGCGTGATTACGCCGCCCCGCAAACTATTGAAATGGGAAT‑gg > 1:2174100/1‑71 (MQ=255)aTTCTGCTCGATGAAAATCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑gg > 1:773490/1‑71 (MQ=255)aTTCTGCTCGATGAAAATCAGCCAGGAAGTGCGTGAATACGCCGCCACGCAAACTATTGAAATGGGAAT‑gg > 1:1694877/1‑71 (MQ=255) tCTGCTCGATGAAAATCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCg > 1:397357/1‑71 (MQ=255) tCTGCTCGATGAAAATCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCg > 1:3479005/1‑71 (MQ=255) cTCGATGAAAATCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGa < 1:2053010/69‑1 (MQ=255) aTGAAAATCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑ggg > 1:991812/1‑61 (MQ=255) aTGAAAATCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑ggg > 1:1440515/1‑61 (MQ=255) aTGAAAATCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑ggg > 1:3258475/1‑61 (MQ=255) aTGAAAATCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑ggg > 1:2206624/1‑61 (MQ=255) aaaTCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGgag < 1:3382677/71‑1 (MQ=255) aaaTCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGgag < 1:3716740/71‑1 (MQ=255) aaaTCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGgag < 1:3149213/71‑1 (MQ=255) aaaTCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGgag < 1:1850616/71‑1 (MQ=255) cAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCg < 1:4067524/71‑1 (MQ=255) cAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCg < 1:2674871/71‑1 (MQ=255) aaGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGc > 1:1486838/1‑70 (MQ=255) aaGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGc > 1:1089458/1‑70 (MQ=255) aaGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGc > 1:969450/1‑70 (MQ=255) aaGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGc > 1:3974421/1‑70 (MQ=255) gTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGc > 1:241195/1‑68 (MQ=255) tGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGa < 1:3065825/46‑1 (MQ=255) tGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCaga < 1:353846/71‑1 (MQ=255) cGTGATTACGCCGCCACGCAAACTATTGAAATGGGAATGGGCGGATATGTCGGAGAACTTCCGTGCCagag > 1:3739253/1‑71 (MQ=255) cGTGATTACGCCGCCACGCAAACTATTGAAATGGGAATGGGCGGATATGTCGGAGAACTTCCGTGCCagag > 1:3606141/1‑71 (MQ=255) cGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAgg > 1:2638891/1‑71 (MQ=255) gTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGa < 1:3898269/43‑1 (MQ=255) gTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGa < 1:943003/43‑1 (MQ=255) gTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGa < 1:914943/43‑1 (MQ=255) tGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCg < 1:470106/71‑1 (MQ=255) aTTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCg > 1:1808428/1‑69 (MQ=255) aTTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCg > 1:280270/1‑69 (MQ=255) aTTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCg > 1:1424000/1‑69 (MQ=255) ttACGCCGCCACGCAAACTATTGAAATGGGAATGGGCGGATATGTCGgag < 1:691334/50‑1 (MQ=255) ttACGCCGCCACGCAAACTATTGAAATGGGAATGGGCGGATATGTCGgag < 1:3220114/50‑1 (MQ=255) tACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGgag > 1:342521/1‑70 (MQ=255) cgccgcCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGaa < 1:356477/70‑1 (MQ=255) cgccgcCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGaa < 1:133388/70‑1 (MQ=255) cgccgcCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGaa < 1:21143/70‑1 (MQ=255) cgccgcCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGaa < 1:3390779/70‑1 (MQ=255) cgccgcCACGCAAACAATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGaa < 1:3506494/70‑1 (MQ=255) cgccACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCt < 1:421637/71‑1 (MQ=255) cgccACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCt < 1:4093436/71‑1 (MQ=255) cgccACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCt < 1:1757207/71‑1 (MQ=255) cgccACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCt < 1:3001451/71‑1 (MQ=255) cgccACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCt < 1:2996637/71‑1 (MQ=255) ccACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTAc < 1:96265/71‑1 (MQ=255) cGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTAcc > 1:811097/1‑69 (MQ=255) cGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTAcc > 1:2774728/1‑69 (MQ=255) cGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTAcc > 1:2385730/1‑69 (MQ=255) cAAACTATTGAAATGGGAATGGGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGc > 1:1243963/1‑71 (MQ=255) cAAACTATTGAAATGGGAATGGGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGc > 1:1631931/1‑71 (MQ=255) cAAACTATTGAAATGGGAATGGGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGc > 1:2309040/1‑71 (MQ=255) cAAACTATTGAAATGGGAAT‑GGCGTATATGTCGGAGAAc > 1:1683098/1‑39 (MQ=255) cAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAAc > 1:547771/1‑39 (MQ=255) cAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAAc > 1:1735568/1‑39 (MQ=255) cAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAAc > 1:393863/1‑39 (MQ=255) cAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAAc > 1:2718342/1‑39 (MQ=255) cAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAAc > 1:3305944/1‑39 (MQ=255) cAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGc > 1:694727/1‑70 (MQ=255) cAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGCg > 1:673735/1‑71 (MQ=255) aaCTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCt < 1:92413/62‑1 (MQ=255) aaCTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCt < 1:2083534/62‑1 (MQ=255) gAAATGGGAATGGGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGCGTAaggag > 1:1971250/1‑70 (MQ=255) gAAATGGGAATGGGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGCGTAaggag > 1:1975887/1‑70 (MQ=255) gAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGc > 1:4135369/1‑38 (MQ=255) gAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGc > 1:3041571/1‑38 (MQ=255) gAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGc > 1:3526591/1‑38 (MQ=255) gAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGc > 1:553045/1‑38 (MQ=255) gAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGc > 1:4022643/1‑38 (MQ=255) gAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTAc > 1:2142108/1‑57 (MQ=255) gAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGCGTAaggagg > 1:3643070/1‑70 (MQ=255) gAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGCGTAaggagg > 1:3465544/1‑70 (MQ=255) gAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGCGTAaggagg > 1:2491651/1‑70 (MQ=255) aaTGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGaa < 1:3825097/49‑1 (MQ=255) tGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGaaa < 1:301375/48‑1 (MQ=255) gggAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTa > 1:3716576/1‑51 (MQ=255) ggAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATc > 1:3462727/1‑48 (MQ=255) ggAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATc > 1:1043231/1‑48 (MQ=255) ggAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATc > 1:1687915/1‑48 (MQ=255) ggAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATc > 1:983146/1‑48 (MQ=255) ggAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATc > 1:4112723/1‑48 (MQ=255) | ATTCTGCTCGATGAAAATCAGCCAGGAAGTGCGTGATTACGCCGCCACGCAAACTATTGAAATGGGAAT‑GGCGGATATGTCGGAGAACTTCCGTGCCAGAGGCGGAGAAATCTACCTGCGTAAGGAGG > minE/2248048‑2248175 |
| 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: ‑ |