Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 117,332 G→A 51.8% A326A (GCG→GCA acnB → bifunctional aconitate hydratase 2 and 2‑methylisocitrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE117,3320GA51.8% ‑3.4 / 47.1 54A326A (GCG→GCAacnBbifunctional aconitate hydratase 2 and 2‑methylisocitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (5/21);  new base A (18/10);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

GGCGACGTGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAGTGCGTGCTGGTGGCCGTATTCCG  >  minE/117264‑117392
                                                                    |                                                            
ggCGACGTGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCg                                                              >  1:2298003/1‑69 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:10862/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:2520157/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:2755614/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:2397618/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:220634/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:2089277/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:2044244/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:1873373/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:3186796/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:1746714/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:3269382/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:1633749/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:1699686/69‑1 (MQ=255)
       tGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCg                                                       <  1:413944/69‑1 (MQ=255)
            gACGTTTACCCGTATAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:2978915/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGa                                                      >  1:2942413/1‑65 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:2594312/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:973022/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:1041813/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:824088/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:655810/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:618364/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:2876998/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:2898345/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:488164/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:3158091/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:3345935/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:2488664/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:1174130/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:1983068/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:1826355/1‑68 (MQ=255)
            gACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCAACCTTCGAACt                                                   >  1:1501352/1‑68 (MQ=255)
                  tACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCGAACTGAAAAcc                                            >  1:1271015/1‑69 (MQ=255)
                               gTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCGAACTGAAAACCGACGTGCtgattg                               >  1:306745/1‑69 (MQ=255)
                                           aCCACGAAACCGGCGAACTGCTGGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAgtgcgtg                   >  1:3026210/1‑69 (MQ=255)
                                               cGAAACCGGCGAACTGCTGGCAACCTTCGAACTGAAAACCGACg                                        <  1:2522477/44‑1 (MQ=255)
                                               cGAAACCGGCGAACTGCTGGCAACCTTCGAACTGAAAACCGACg                                        <  1:112077/44‑1 (MQ=255)
                                               cGAAACCGGCGAACTGCTGGCAACCTTCGAACTGAAAACCGACg                                        <  1:928725/44‑1 (MQ=255)
                                               cGAAACCGGCGAACTGCTGGCAACCTTCGAACTGAAAACCGACg                                        <  1:1251996/44‑1 (MQ=255)
                                               cGAAACCGGCGAACTGCTGGCAACCTTCGAACTGAAAACCGACg                                        <  1:632081/44‑1 (MQ=255)
                                               cGAAACCGGCGAACTGCTGGCAACCTTCGAACTGAAAACCGACg                                        <  1:1263511/44‑1 (MQ=255)
                                               cGAAACCGGCGAACTGCTGGCAACCTTCGAACTGAAAACCGACg                                        <  1:334026/44‑1 (MQ=255)
                                               cGAAACCGGCGAACTGCTGGCAACCTTCGAACTGAAAACCGACg                                        <  1:2183865/44‑1 (MQ=255)
                                               cGAAACCGGCGAACTGCTGGCAACCTTCGAACTGAAAACCGACg                                        <  1:2673070/44‑1 (MQ=255)
                                                     cGGCGAACTGCTGGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAGTGCGTTCTGGTGGCCg         <  1:499849/69‑1 (MQ=255)
                                                     cGGCGAACTGCTGGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAGTGCGTGCTGGTGGCCg         <  1:402081/69‑1 (MQ=255)
                                                     cGGCGAACTGCTGGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAGTGCGTGCTGGTGGCCg         <  1:528379/69‑1 (MQ=255)
                                                     cGGCGAACTGCTGGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAGTGCGTGCTGGTGGCCg         <  1:2100391/69‑1 (MQ=255)
                                                     cGGCGAACTGCTGGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAGTGCGTGCTGGTGGCCg         <  1:2784781/69‑1 (MQ=255)
                                                     cGGCGAACTGCTGGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAGTGCGTGCTGGTGGCCg         <  1:2534991/69‑1 (MQ=255)
                                                       gCGAACTGCTGGCAACCTTCGAACTGAAAACCGACg                                        <  1:1465589/36‑1 (MQ=255)
                                                            ctgctgGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAGTGCGTGCTGGTGGCCGTATTCCg  >  1:1663861/1‑69 (MQ=255)
                                                            ctgctgGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAGTGCGTGCTGGTGGCCGTATTCCg  >  1:1936645/1‑69 (MQ=255)
                                                                tgGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAg                         <  1:2075456/42‑1 (MQ=255)
                                                                    |                                                            
GGCGACGTGATTGACGTTTACCCGTACAAAGGTGAAGTGCGTAACCACGAAACCGGCGAACTGCTGGCGACCTTCGAACTGAAAACCGACGTGCTGATTGATGAAGTGCGTGCTGGTGGCCGTATTCCG  >  minE/117264‑117392

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: