Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,185 A→T 23.9% V86V (GTT→GTA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,1850AT23.9% 174.8 / 57.3 92V86V (GTT→GTAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (34/36);  new base T (13/9);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCG  >  minE/1934116‑1934255
                                                                     |                                                                      
ataAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTaa                                                                       >  1:3361187/1‑71 (MQ=255)
ataAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTaa                                                                       >  1:4515001/1‑71 (MQ=255)
          tCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGc                                                             >  1:2419192/1‑71 (MQ=255)
           ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTTCAGGGTAACGAAAGACGc                                                             <  1:4037873/70‑1 (MQ=255)
           ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGc                                                             <  1:315032/70‑1 (MQ=255)
           ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCt                                                            >  1:830366/1‑71 (MQ=255)
           ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCt                                                            >  1:3594207/1‑71 (MQ=255)
           ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCt                                                            >  1:43532/1‑71 (MQ=255)
             gggAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCttt                                                          <  1:2600185/71‑1 (MQ=255)
              ggAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTa                                                         <  1:4818060/71‑1 (MQ=255)
                   tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                    <  1:3923450/71‑1 (MQ=255)
                   tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                    <  1:1369846/71‑1 (MQ=255)
                      tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataat                                                 <  1:4953092/71‑1 (MQ=255)
                      tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                    <  1:1506555/68‑1 (MQ=255)
                      tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                    <  1:1500069/68‑1 (MQ=255)
                      tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                    <  1:4254397/68‑1 (MQ=255)
                       cACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGa                                                                >  1:3729867/1‑55 (MQ=255)
                        aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                               >  1:5192515/1‑71 (MQ=255)
                        aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                               >  1:3670634/1‑71 (MQ=255)
                        aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                               >  1:2305793/1‑71 (MQ=255)
                        aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                               >  1:4089260/1‑71 (MQ=255)
                        aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                               >  1:5635555/1‑71 (MQ=255)
                        aCGGGCAACGTACAGATGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                                                  >  1:2103692/1‑68 (MQ=255)
                         cttgCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATac                                               >  1:4010170/4‑70 (MQ=255)
                         cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                    <  1:3553945/49‑1 (MQ=255)
                         cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                    <  1:1083991/49‑1 (MQ=255)
                         cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                    <  1:5206876/49‑1 (MQ=255)
                         cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                    <  1:5608964/49‑1 (MQ=255)
                         cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                    <  1:3261681/65‑1 (MQ=255)
                         cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                               >  1:5431537/1‑70 (MQ=255)
                          ttgCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATacga                                             >  1:1882957/3‑71 (MQ=255)
                          gggCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacga                                             <  1:4257391/71‑1 (MQ=255)
                            gCAACGTACAGTTGAATCCGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTccc                                                       >  1:2586352/1‑59 (MQ=255)
                            gCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTccc                                                       >  1:2435601/1‑59 (MQ=255)
                            gCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTccc                                                       >  1:3486407/1‑59 (MQ=255)
                            gCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATacga                                             >  1:97317/1‑69 (MQ=255)
                             cAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACgg                                          >  1:1178934/1‑71 (MQ=255)
                              aaCGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCa                                                     <  1:1086633/59‑1 (MQ=255)
                               aCGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGt                                        >  1:3189115/1‑71 (MQ=255)
                               aCGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGt                                        >  1:1821742/1‑71 (MQ=255)
                               aCGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGt                                        <  1:1442921/71‑1 (MQ=255)
                                cGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGt                                        >  1:2923742/1‑70 (MQ=255)
                                 gTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                                                  <  1:2922051/43‑1 (MQ=255)
                                 gTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                                                  <  1:4147612/43‑1 (MQ=255)
                                 gTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                    <  1:5027833/41‑1 (MQ=255)
                                 gTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtca                                      <  1:2699140/71‑1 (MQ=255)
                                  tACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGtcat                                     <  1:459105/71‑1 (MQ=255)
                                  tACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACg                                                              <  1:3633276/46‑1 (MQ=255)
                                  tACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcat                                     >  1:2695632/1‑71 (MQ=255)
                                  tACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcat                                     <  1:192176/71‑1 (MQ=255)
                                   aCAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacg                                              >  1:1351298/1‑61 (MQ=255)
                                    cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaa                                                                   >  1:256010/1‑39 (MQ=255)
                                    cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCt                                                            >  1:1692253/1‑46 (MQ=255)
                                    cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatca                                   >  1:2739834/1‑71 (MQ=255)
                                    cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatca                                   >  1:2369524/1‑71 (MQ=255)
                                    cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatc                                    <  1:5568952/70‑1 (MQ=255)
                                     aGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTcatcat                                  <  1:4756993/71‑1 (MQ=255)
                                        tGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCg                               >  1:65613/1‑71 (MQ=255)
                                         gAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCttt                                                          <  1:2181078/43‑1 (MQ=255)
                                         gAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCgg                              >  1:2526588/1‑71 (MQ=255)
                                         gAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCgg                              >  1:2287907/1‑71 (MQ=255)
                                         gAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCg                               <  1:3245899/70‑1 (MQ=255)
                                         gAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCg                               >  1:1549295/1‑70 (MQ=255)
                                          aaTGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCa                                                     <  1:3626458/47‑1 (MQ=255)
                                          aaTACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCgg                              <  1:4012404/70‑1 (MQ=255)
                                          aaTACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCg                               <  1:4738092/69‑1 (MQ=255)
                                            tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                                                  <  1:2251878/48‑1 (MQ=255)
                                            tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                                                  <  1:5474961/48‑1 (MQ=255)
                                            tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgacg                                           >  1:3691850/1‑55 (MQ=255)
                                            tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgacg                                           >  1:5491071/1‑55 (MQ=255)
                                            tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgac                                            <  1:759894/54‑1 (MQ=255)
                                            tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgac                                            <  1:63102/54‑1 (MQ=255)
                                              cGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAg                          >  1:1058439/1‑70 (MQ=255)
                                              cGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAg                          >  1:3759086/1‑70 (MQ=255)
                                              cGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAg                          >  1:4926516/1‑70 (MQ=255)
                                                   gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                       >  1:932022/1‑68 (MQ=255)
                                                   gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                       >  1:4279799/1‑68 (MQ=255)
                                                   gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                       >  1:5591560/1‑68 (MQ=255)
                                                    ccAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                                                  >  1:3067599/1‑40 (MQ=255)
                                                    ccAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACgg                    >  1:1209557/1‑70 (MQ=255)
                                                    ccAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGc                   <  1:587097/71‑1 (MQ=255)
                                                      aaCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGc                         <  1:3508765/63‑1 (MQ=255)
                                                        cGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACg                     >  1:5169788/1‑65 (MQ=255)
                                                        cGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACg                     >  1:5157556/1‑65 (MQ=255)
                                                         atCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAAcc                     <  1:4339776/63‑2 (MQ=255)
                                                         atCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAAcc                     <  1:5623845/63‑2 (MQ=255)
                                                         atCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAAcc                     <  1:234300/63‑2 (MQ=255)
                                                           ccTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAAcc                     <  1:3715333/62‑2 (MQ=255)
                                                            cTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTca            <  1:3961549/70‑2 (MQ=255)
                                                             tGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGa           >  1:5539454/1‑70 (MQ=255)
                                                                 gggTAACGAAAGACGCTTTACCCATAATACGACGGGtcat                                     <  1:5414647/40‑1 (MQ=255)
                                                                  ggTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCa     <  1:5377813/71‑1 (MQ=255)
                                                                    tAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGc                         >  1:4153184/1‑49 (MQ=255)
                                                                    tAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGc                         >  1:5176971/1‑49 (MQ=255)
                                                                     aaCGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAgcg  >  1:4978297/1‑71 (MQ=255)
                                                                     |                                                                      
ATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCG  >  minE/1934116‑1934255

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: