Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,320,260 G→A 51.6% D429D (GAC→GAT gnd ← gluconate‑6‑phosphate dehydrogenase, decarboxylating

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,320,2600GA51.6% ‑1.9 / 85.0 93D429D (GAC→GATgndgluconate‑6‑phosphate dehydrogenase, decarboxylating
Reads supporting (aligned to +/- strand):  ref base G (11/34);  new base A (27/21);  total (38/55)
Fisher's exact test for biased strand distribution p-value = 2.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAAC  >  minE/1320197‑1320325
                                                               |                                                                 
cgcACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGt                                                               >  1:2972426/1‑68 (MQ=255)
cgcACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGt                                                               >  1:1888759/1‑68 (MQ=255)
cgcACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGt                                                               >  1:2094226/1‑68 (MQ=255)
cgcACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGt                                                               >  1:2286165/1‑68 (MQ=255)
cgcACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGt                                                               >  1:2300729/1‑68 (MQ=255)
cgcACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGt                                                               >  1:163925/1‑68 (MQ=255)
cgcACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGt                                                               >  1:2656958/1‑68 (MQ=255)
cgcACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGt                                                               >  1:2686574/1‑68 (MQ=255)
cgcACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGt                                                               >  1:1325899/1‑68 (MQ=255)
cgcACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGt                                                               >  1:1103941/1‑68 (MQ=255)
              cACGCTGTGCCTGGATCATGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:472564/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:367967/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:3162677/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:1075354/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:289854/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:278231/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:586664/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:689486/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:818907/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:839867/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:2103559/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:2067971/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:905348/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:1628689/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:1110614/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:1147302/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:1810862/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:1304321/69‑3 (MQ=255)
              cACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:1321262/69‑3 (MQ=255)
               aCGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:1356613/68‑3 (MQ=255)
               aCGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGcag                                                <  1:250438/68‑3 (MQ=255)
                   tgtgCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGgag                                           <  1:2073029/69‑1 (MQ=255)
                             tCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGTTCGGa                                  <  1:1428196/68‑1 (MQ=255)
                             tCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGa                                  <  1:616646/68‑1 (MQ=255)
                             tCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGa                                  <  1:2968685/68‑1 (MQ=255)
                             tCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGa                                  <  1:2919878/68‑1 (MQ=255)
                             tCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGa                                  <  1:24010/68‑1 (MQ=255)
                                 gTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:2978700/69‑1 (MQ=255)
                                 gTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:1698225/69‑1 (MQ=255)
                                 gTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:2401741/69‑1 (MQ=255)
                                 gTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:1619640/69‑1 (MQ=255)
                                 gTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:744317/69‑1 (MQ=255)
                                 gTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:994640/69‑1 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCATCCGCAGCGGCGAAGGTCGGAAc                                >  1:1393658/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCACCCGCAGCGGCGAAGGTCGGAAc                                >  1:2496537/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:2028392/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:3212671/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:3243763/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:357020/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:2082954/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:376779/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:460937/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:131263/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:478607/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:3302517/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:1217372/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:2010223/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:75876/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:835634/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:1954267/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:969489/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:1940766/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:2690538/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:2866260/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:2908016/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:2798566/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:2142981/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:2987697/1‑59 (MQ=255)
                                        ggcagAACAGCGGCACGGTACCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:623894/1‑59 (MQ=255)
                                        ggcagAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCg                      <  1:874966/69‑1 (MQ=255)
                                        ggcagAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCg                      <  1:3431120/69‑1 (MQ=255)
                                        ggcagAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCg                      <  1:1844395/69‑1 (MQ=255)
                                        ggcagAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCg                      <  1:2471536/69‑1 (MQ=255)
                                        ggcagAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCg                      <  1:1836601/69‑1 (MQ=255)
                                        ggcacAACAGCGGCACGGTAGCTATCGTAATAGGCAACCGCAGCGGCGAAGGTCGGAAc                                >  1:3300472/1‑59 (MQ=255)
                                             aaCAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTg                 <  1:2655768/69‑1 (MQ=255)
                                             aaCAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTg                 <  1:1594726/69‑1 (MQ=255)
                                             aaCAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTg                 <  1:3227149/69‑1 (MQ=255)
                                             aaCAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTg                 <  1:1042067/69‑1 (MQ=255)
                                               cagcagCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:275536/55‑1 (MQ=255)
                                               cagcagCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:716386/55‑1 (MQ=255)
                                               cagcagCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:256488/55‑1 (MQ=255)
                                               cagcagCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:2367945/55‑1 (MQ=255)
                                               cagcagCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:425691/55‑1 (MQ=255)
                                               cagcagCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:170788/55‑1 (MQ=255)
                                               cagcagCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCgg                             <  1:3332856/55‑1 (MQ=255)
                                                         gTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGc     >  1:3366656/1‑69 (MQ=255)
                                                            gCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAAc  <  1:2798337/69‑1 (MQ=255)
                                                            gCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAAc  <  1:2009430/69‑1 (MQ=255)
                                                            gCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAAc  <  1:3214456/69‑1 (MQ=255)
                                                            gCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAAc  <  1:314369/69‑1 (MQ=255)
                                                            gCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAAc  <  1:1782148/69‑1 (MQ=255)
                                                              tGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAAc  <  1:2888531/67‑1 (MQ=255)
                                                               |                                                                 
CGCACCAAAATAGTCACGCTGTGCCTGGATCAGGTTCGCAGGCAGAACAGCAGCACGGTAGCTGTCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAAC  >  minE/1320197‑1320325

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: