Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,425,822 A→T 5.9% S90S (TCT→TCA rffG ← dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,425,8220AT 5.9% 291.5 / 8.2 102S90S (TCT→TCArffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base A (49/47);  new base T (3/3);  total (52/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

Only 100 of 103 total aligned reads displayed.

CCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATACGCGT  >  minE/2425753‑2425891
                                                                     |                                                                     
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)
          tCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCaa                                                             >  1:834530/1‑70 (MQ=255)
           cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                                                           >  1:3807353/1‑71 (MQ=255)
           cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                                                           >  1:1823733/1‑71 (MQ=255)
           cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                                                           >  1:1956674/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)
                 aTGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCtt                                                     >  1:1387052/1‑71 (MQ=255)
                  tGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCttt                                                    <  1:2269879/71‑1 (MQ=255)
                  tGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCttt                                                    <  1:1103193/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:59207/1‑71 (MQ=255)
                         gTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGc                                             >  1:2278821/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:2659163/1‑70 (MQ=255)
                                cAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATg                                       >  1:1830248/1‑70 (MQ=255)
                                cAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGc                                      <  1:2277127/71‑1 (MQ=255)
                                cAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGc                                      <  1:154459/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:3920947/1‑71 (MQ=255)
                                           tCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTc                           >  1:3741151/1‑71 (MQ=255)
                                           tCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTc                           >  1:2886854/1‑71 (MQ=255)
                                           tCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTc                           >  1:359633/1‑71 (MQ=255)
                                           tCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTc                           >  1:1634468/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)
                                               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:3074137/1‑71 (MQ=255)
                                                aaaCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCt                      >  1:3240864/1‑71 (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:1420360/69‑1 (MQ=255)
                                                        ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCt                >  1:2919334/1‑69 (MQ=255)
                                                        ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCt                >  1:2370404/1‑69 (MQ=255)
                                                        ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTc               <  1:1531255/70‑1 (MQ=255)
                                                        ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTc               <  1:290285/70‑1 (MQ=255)
                                                        ccGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTc               <  1:3039585/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:1491050/71‑1 (MQ=255)
                                                              cGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAAt        <  1:1474585/71‑1 (MQ=255)
                                                              cGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAAt        <  1:1227328/71‑1 (MQ=255)
                                                                tCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCt                <  1:2861656/61‑1 (MQ=255)
                                                                tCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATAc      >  1:686029/1‑71 (MQ=255)
                                                                 cAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATAcg     >  1:1501768/1‑71 (MQ=255)
                                                                  aaTAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATAcgc    >  1:2982325/1‑71 (MQ=255)
                                                                  aaTAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATAcgc    >  1:2442891/1‑71 (MQ=255)
                                                                  aaTAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATAcgc    >  1:1796504/1‑71 (MQ=255)
                                                                  aaTAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATAc      <  1:1064342/69‑1 (MQ=255)
                                                                   aTAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGa                                  >  1:2709219/1‑40 (MQ=255)
                                                                   aTAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGa                                  >  1:1189468/1‑40 (MQ=255)
                                                                    tAGAACGGTCAACATGGCTTTCGGCTGCCAGATGc                                      <  1:2141930/35‑1 (MQ=255)
                                                                    tAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATACGCGt  <  1:3792981/71‑1 (MQ=255)
                                                                     |                                                                     
CCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGATGCTCAGTGAATACGCGT  >  minE/2425753‑2425891

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: