Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,872,291 G→A 9.0% G181D (GGT→GAT) ‡ ulaF → L‑ribulose 5‑phosphate 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,872,2910GA 9.0% 251.1 / 10.6 89G181D (GGT→GAT) ‡ulaFL‑ribulose 5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (28/53);  new base A (6/2);  total (34/55)
Fisher's exact test for biased strand distribution p-value = 5.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 101 total aligned reads displayed.

TTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACA  >  minE/2872226‑2872357
                                                                 |                                                                  
ttGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtggt                                                                >  1:2306317/1‑70 (MQ=255)
ttGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtggt                                                                >  1:923546/1‑70 (MQ=255)
    ggtgTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATg                                                            <  1:2067678/70‑1 (MQ=255)
     gtgtATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATgg                                                           <  1:2211007/70‑1 (MQ=255)
        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtgg                                                                 >  1:4805118/1‑61 (MQ=255)
        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtgg                                                                 >  1:2581399/1‑61 (MQ=255)
        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtgg                                                                 >  1:3860845/1‑61 (MQ=255)
        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga                                                       <  1:4604790/71‑1 (MQ=255)
        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaag                                                        >  1:4285458/1‑70 (MQ=255)
          tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaagaa                                                      <  1:4547312/70‑1 (MQ=255)
            cccACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGgaa                                                         >  1:848685/3‑65 (MQ=255)
            cccACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGtc                                                   >  1:4479915/3‑70 (MQ=255)
            cccACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGtc                                                   >  1:2887968/3‑70 (MQ=255)
            cccACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGt                                                    >  1:1172157/3‑70 (MQ=255)
            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt                                                    <  1:4557762/70‑1 (MQ=255)
            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTg                                                   <  1:4289954/71‑1 (MQ=255)
            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTg                                                   >  1:374889/1‑71 (MQ=255)
            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTg                                                   <  1:3622580/71‑1 (MQ=255)
             ccACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGtc                                                   >  1:3668662/2‑69 (MQ=255)
              cACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGTCGc                                                 <  1:1034563/71‑1 (MQ=255)
               aCGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt                                                    >  1:3283436/1‑67 (MQ=255)
               aCGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt                                                    >  1:5679899/1‑67 (MQ=255)
               aCGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGTCGct                                                >  1:648511/1‑70 (MQ=255)
               aCGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGTCGct                                                <  1:744477/71‑2 (MQ=255)
                cGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCg                                                >  1:3421720/1‑70 (MQ=255)
                cGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGa                                               <  1:394521/71‑1 (MQ=255)
                cGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGa                                               <  1:5533142/71‑1 (MQ=255)
                 gggCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGa                                               >  1:80711/1‑70 (MQ=255)
                   gCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGAt                                                             >  1:4939775/1‑54 (MQ=255)
                   gCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                                            >  1:209393/1‑71 (MQ=255)
                   gCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                                            <  1:4340449/71‑1 (MQ=255)
                    ccGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa                                            <  1:4346794/70‑1 (MQ=255)
                      gTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATgg                                         >  1:1590732/1‑71 (MQ=255)
                       ttCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtggt                                                                <  1:3008695/47‑1 (MQ=255)
                         cGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCgg                                                                    >  1:3264103/1‑41 (MQ=255)
                         cGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCgg                                                                    >  1:3268981/1‑41 (MQ=255)
                         cGCCTGGGGGAAAGATGCTCAAGATGCGGTGCATAACGCGGTGGTGa                                                              >  1:3726122/1‑47 (MQ=255)
                          gCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTg                                     <  1:4897729/71‑1 (MQ=255)
                          gCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTg                                     <  1:4904463/71‑1 (MQ=255)
                          gCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTg                                     <  1:4288885/71‑1 (MQ=255)
                           ccTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga                                                       <  1:3710590/52‑1 (MQ=255)
                             tGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGAt                                  >  1:351066/1‑71 (MQ=255)
                             tGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGAt                                  >  1:3644140/1‑71 (MQ=255)
                             tGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGAt                                  >  1:5132382/1‑71 (MQ=255)
                              gggggAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt                                                    <  1:1929678/52‑1 (MQ=255)
                               ggggAAAGATGCTCACGATGCGGTGCATAACGCggtggt                                                                <  1:1243594/39‑1 (MQ=255)
                               ggggAAAGATGCTCACGATGCGGTGCATAACGCggtggt                                                                <  1:4909594/39‑1 (MQ=255)
                               ggggAAAGATGCTCACGATGCGGTGCATAACGCggtggt                                                                <  1:2546126/39‑1 (MQ=255)
                               ggggAAAGATGCTCACGATGCGGTGCATAACGCggtggt                                                                <  1:5049914/39‑1 (MQ=255)
                               ggggAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGAt                                  <  1:1433365/69‑1 (MQ=255)
                                   aaaGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt                                                    <  1:1849441/47‑1 (MQ=255)
                                    aaGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCcgc                           >  1:490196/1‑71 (MQ=255)
                                     aGATTCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt                                                    <  1:4170256/45‑1 (MQ=255)
                                     aGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt                                                    <  1:2794660/45‑1 (MQ=255)
                                     aGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt                                                    <  1:429674/45‑1 (MQ=255)
                                     aGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt                                                    <  1:4422703/45‑1 (MQ=255)
                                     aGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTgg                                                  >  1:4602987/1‑47 (MQ=255)
                                     aGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTgg                                                  >  1:2661147/1‑47 (MQ=255)
                                     aGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTgg                                                  >  1:3316616/1‑47 (MQ=255)
                                     aGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTgg                                                  >  1:148463/1‑47 (MQ=255)
                                     aGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaa                                              <  1:1545982/51‑1 (MQ=255)
                                     aGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaa                                              <  1:1144192/51‑1 (MQ=255)
                                     aGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaa                                              <  1:4427771/51‑1 (MQ=255)
                                     aGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaa                                              <  1:4716890/51‑1 (MQ=255)
                                       aTGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCgg                         >  1:4235227/1‑70 (MQ=255)
                                       aTGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCgg                         >  1:1123374/1‑70 (MQ=255)
                                         gCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAt                      >  1:2695070/1‑71 (MQ=255)
                                          cTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt                     <  1:2629783/71‑1 (MQ=255)
                                          cTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAt                      <  1:796782/70‑1 (MQ=255)
                                           tCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt                     >  1:2634633/1‑70 (MQ=255)
                                           tCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAt                      <  1:66643/69‑1 (MQ=255)
                                             aCGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTaa                   <  1:1990982/70‑1 (MQ=255)
                                             aCGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAAc                  >  1:803648/1‑71 (MQ=255)
                                             aCGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAAc                  >  1:3027528/1‑71 (MQ=255)
                                              cGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTg                                     <  1:1080159/51‑1 (MQ=255)
                                              cGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTg                                     <  1:1450193/51‑1 (MQ=255)
                                              cGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTg                                     <  1:278121/51‑1 (MQ=255)
                                                aTGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCca               >  1:3285736/1‑71 (MQ=255)
                                                aTGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCca               >  1:5442204/1‑71 (MQ=255)
                                                 tGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCcac              <  1:5400436/71‑1 (MQ=255)
                                                   cGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaa                                              <  1:3606104/37‑1 (MQ=255)
                                                   cGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCcgc                           <  1:479265/56‑1 (MQ=255)
                                                   cGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCcgc                           <  1:3874169/56‑1 (MQ=255)
                                                   aggTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaa                                              <  1:2204294/36‑1 (MQ=255)
                                                    ggTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAAc           <  1:4748852/71‑1 (MQ=255)
                                                     gTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt                     <  1:5330901/60‑1 (MQ=255)
                                                     gTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt                     <  1:603351/60‑1 (MQ=255)
                                                     gTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt                     <  1:593340/60‑1 (MQ=255)
                                                     gTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCcac              <  1:2913053/67‑1 (MQ=255)
                                                     gTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCcac              <  1:1541497/67‑1 (MQ=255)
                                                     gTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCcac              <  1:5342703/67‑1 (MQ=255)
                                                      tGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAATTGGCGTg                                     <  1:4648006/43‑1 (MQ=255)
                                                      tGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTg                                     <  1:3872119/43‑1 (MQ=255)
                                                      tGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTg                                     <  1:1786805/43‑1 (MQ=255)
                                                      tGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTg                                     <  1:80033/43‑1 (MQ=255)
                                                      tGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTg                                     <  1:4825754/43‑1 (MQ=255)
                                                        cATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa        <  1:4491733/70‑1 (MQ=255)
                                                            aCGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGccc                             >  1:4066239/1‑45 (MQ=255)
                                                             cgcgGTTGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATcaca  <  1:1232952/71‑1 (MQ=255)
                                                             cgcgGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt                     <  1:1778596/52‑1 (MQ=255)
                                                                 |                                                                  
TTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACA  >  minE/2872226‑2872357

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: