Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,425,816 G→T 6.0% D92E (GAC→GAA) ‡ rffG ← dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,425,8160GT 6.0% 329.4 / 4.0 99D92E (GAC→GAA) ‡rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base G (45/48);  new base T (3/3);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAAT  >  minE/2425747‑2425885
                                                                     |                                                                     
gTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCg                                                                       <  1:3314546/70‑1 (MQ=255)
  aaGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTc                                                                     <  1:1631177/70‑1 (MQ=255)
   aGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCaa                                                                   >  1:687279/1‑71 (MQ=255)
    gCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCaa                                                                   <  1:2441393/70‑1 (MQ=255)
     cccGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCaa                                                                   >  1:529760/1‑69 (MQ=255)
      ccGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAg                                                                <  1:1885767/71‑1 (MQ=255)
       cgcgCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATa                                                                 >  1:115415/1‑69 (MQ=255)
       cgcgCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAg                                                                >  1:2762674/1‑70 (MQ=255)
        gcgcCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGaa                                                              >  1:264388/1‑71 (MQ=255)
          gccgcTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACg                                                            <  1:1797915/71‑1 (MQ=255)
          gccgcTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACg                                                            <  1:699644/71‑1 (MQ=255)
          gccgcTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTCCCGGGCCGTCAATAGAACg                                                            >  1:1662591/1‑71 (MQ=255)
           ccgcTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACg                                                            <  1:1733003/70‑1 (MQ=255)
             gcTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTc                                                         <  1:3654842/71‑1 (MQ=255)
                tCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAAt                                                                  <  1:1894911/59‑1 (MQ=255)
                tCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCaa                                                       >  1:834530/1‑70 (MQ=255)
                 cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTc                                                                     <  1:1397775/55‑1 (MQ=255)
                 cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                                                     >  1:3807353/1‑71 (MQ=255)
                 cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                                                     >  1:1956674/1‑71 (MQ=255)
                 cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                                                     >  1:1823733/1‑71 (MQ=255)
                  aaGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                                                     <  1:2355245/70‑1 (MQ=255)
                  aaGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACAt                                                    <  1:2614328/71‑1 (MQ=255)
                   aGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACAt                                                    <  1:594339/70‑1 (MQ=255)
                      aaTGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACAt                                                    <  1:3639418/67‑1 (MQ=255)
                      aaTGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACAt                                                    <  1:1425877/67‑1 (MQ=255)
                       aTGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGt                                                                      <  1:2374994/48‑1 (MQ=255)
                       aTGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCtt                                               >  1:1387052/1‑71 (MQ=255)
                        tGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCttt                                              <  1:1103193/71‑1 (MQ=255)
                        tGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCttt                                              <  1:2269879/71‑1 (MQ=255)
                         gTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTc                                             >  1:2844453/1‑71 (MQ=255)
                          tataAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTc                                             <  1:2995065/70‑1 (MQ=255)
                          cataAGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCttt                                              >  1:3209261/2‑69 (MQ=255)
                          cataAGTACCAACAATATTGGTTTCAATAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTcag                                           <  1:861050/70‑3 (MQ=255)
                             aaGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCgg                                           <  1:762616/69‑1 (MQ=255)
                             aaGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGc                                          >  1:1957527/1‑70 (MQ=255)
                              aGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTg                                        <  1:172685/71‑1 (MQ=255)
                              aGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACg                                                            >  1:2668448/1‑51 (MQ=255)
                               gTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTc                                                         <  1:3487039/53‑1 (MQ=255)
                               gTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGc                                       >  1:2278821/1‑71 (MQ=255)
                               gTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGc                                       >  1:59207/1‑71 (MQ=255)
                                tCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACAt                                                    >  1:3646852/1‑57 (MQ=255)
                                tCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGcc                                      <  1:2126022/71‑1 (MQ=255)
                                  ccGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGc                                       >  1:1535743/1‑68 (MQ=255)
                                  ccGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCa                                     <  1:851055/70‑1 (MQ=255)
                                  ccGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAg                                    <  1:3323253/71‑1 (MQ=255)
                                    gACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTc                                                         >  1:765827/1‑48 (MQ=255)
                                     aCAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAGGTg                                 <  1:764629/71‑1 (MQ=255)
                                      cAATGTTGGTTTCAATAAACTCTGCCGGGCCGTCAATAGAACGGTCAAc                                                      >  1:3013270/1‑49 (MQ=255)
                                      cAATGTTGGTTTCAATAAACGCTGCCTGGCCGTCAATAGAACGGTCAACAt                                                    <  1:145190/51‑1 (MQ=255)
                                      cAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATg                                 >  1:1830248/1‑70 (MQ=255)
                                      cAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATg                                 >  1:2659163/1‑70 (MQ=255)
                                      cAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGc                                <  1:154459/71‑1 (MQ=255)
                                      cAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGc                                <  1:2277127/71‑1 (MQ=255)
                                       aaTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCa                               <  1:3150288/71‑1 (MQ=255)
                                        aTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCAt                              >  1:3006506/1‑71 (MQ=255)
                                           ttGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTcag                                           <  1:2282845/55‑3 (MQ=255)
                                           ttGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCttt                                              <  1:1346486/52‑1 (MQ=255)
                                           ttGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAg                                    <  1:1890879/62‑1 (MQ=255)
                                             ggTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGAt                                  <  1:2775507/62‑1 (MQ=255)
                                              gTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGacac                         <  1:1047876/70‑1 (MQ=255)
                                               tttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGcc                                      <  1:3461569/56‑1 (MQ=255)
                                                ttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGc                                <  1:3480487/61‑1 (MQ=255)
                                                 tCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTc                     >  1:2886854/1‑71 (MQ=255)
                                                 tCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTc                     >  1:3920947/1‑71 (MQ=255)
                                                 tCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTc                     >  1:3741151/1‑71 (MQ=255)
                                                 tCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTc                     >  1:1634468/1‑71 (MQ=255)
                                                 tCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTc                     >  1:359633/1‑71 (MQ=255)
                                                  cAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCgg                                           >  1:734671/1‑48 (MQ=255)
                                                  cAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATTACACAGTCt                    >  1:531413/1‑71 (MQ=255)
                                                   aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCgg                                           <  1:405941/47‑1 (MQ=255)
                                                     tAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCg                                            <  1:1968123/44‑1 (MQ=255)
                                                     tAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGt                      >  1:2637057/1‑66 (MQ=255)
                                                     tAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGc                 >  1:2040312/1‑71 (MQ=255)
                                                      aaaCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTgg                  >  1:2017319/1‑69 (MQ=255)
                                                      aaaCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTgg                  >  1:3948973/1‑69 (MQ=255)
                                                      aaaCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTg                   >  1:2850337/1‑68 (MQ=255)
                                                      aaaCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCt                >  1:3240864/1‑71 (MQ=255)
                                                      aaaCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCt                >  1:3074137/1‑71 (MQ=255)
                                                      aaaCGCTGCCGGGCCGTCAATAGAACGGTCAAAATGGCTTTCGGCTGCCAGATGCATGACACAGTCTgg                  >  1:982845/1‑69 (MQ=255)
                                                       aaCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCtt                                               >  1:731789/1‑39 (MQ=255)
                                                        aCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGAt                                  >  1:882756/1‑51 (MQ=255)
                                                        aCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGAt                                  >  1:1190247/1‑51 (MQ=255)
                                                          gCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt                    <  1:1771928/63‑1 (MQ=255)
                                                              ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGAt             >  1:3786143/1‑66 (MQ=255)
                                                              ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCt          <  1:1513662/69‑1 (MQ=255)
                                                              ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCt          >  1:2919334/1‑69 (MQ=255)
                                                              ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCt          >  1:2370404/1‑69 (MQ=255)
                                                              ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCt          <  1:1420360/69‑1 (MQ=255)
                                                              ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTc         <  1:3039585/70‑1 (MQ=255)
                                                              ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTc         <  1:1531255/70‑1 (MQ=255)
                                                              ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTc         <  1:290285/70‑1 (MQ=255)
                                                                gggCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTc         >  1:172826/1‑68 (MQ=255)
                                                                 ggCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTg     <  1:3205400/71‑1 (MQ=255)
                                                                  gCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGa    <  1:2918446/71‑1 (MQ=255)
                                                                   ccGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCa        >  1:2026529/1‑66 (MQ=255)
                                                                    cGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGaa   >  1:1681617/1‑70 (MQ=255)
                                                                    cGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGa    <  1:967302/69‑1 (MQ=255)
                                                                    cGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAAt  <  1:1474585/71‑1 (MQ=255)
                                                                    cGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAAt  <  1:1227328/71‑1 (MQ=255)
                                                                    cGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAAt  <  1:1491050/71‑1 (MQ=255)
                                                                     |                                                                     
GTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAAT  >  minE/2425747‑2425885

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: