Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,188 G→A 22.9% F85F (TTC→TTT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,1880GA22.9% 178.9 / 48.3 96F85F (TTC→TTTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base A (13/9);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

TCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCT  >  minE/1934126‑1934257
                                                              |                                                                     
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:4254397/68‑1 (MQ=255)
            tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                      <  1:1500069/68‑1 (MQ=255)
            tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                      <  1:1506555/68‑1 (MQ=255)
             cACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGa                                                                  >  1:3729867/1‑55 (MQ=255)
              aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                                 >  1:4089260/1‑71 (MQ=255)
              aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                                 >  1:5192515/1‑71 (MQ=255)
              aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                                 >  1:2305793/1‑71 (MQ=255)
              aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                                 >  1:5635555/1‑71 (MQ=255)
              aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                                 >  1:3670634/1‑71 (MQ=255)
              aCGGGCAACGTACAGATGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                                                    >  1:2103692/1‑68 (MQ=255)
               cttgCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATac                                                 >  1:4010170/4‑70 (MQ=255)
               cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      <  1:1083991/49‑1 (MQ=255)
               cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      <  1:3553945/49‑1 (MQ=255)
               cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      <  1:5608964/49‑1 (MQ=255)
               cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      <  1:5206876/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:1821742/1‑71 (MQ=255)
                     aCGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGt                                          >  1:3189115/1‑71 (MQ=255)
                     aCGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGt                                          <  1:1442921/71‑1 (MQ=255)
                      cGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGt                                          >  1:2923742/1‑70 (MQ=255)
                       gTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                                                    <  1:4147612/43‑1 (MQ=255)
                       gTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                                                    <  1:2922051/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:2369524/1‑71 (MQ=255)
                          cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatca                                     >  1:2739834/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:2287907/1‑71 (MQ=255)
                               gAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCgg                                >  1:2526588/1‑71 (MQ=255)
                               gAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCg                                 >  1:1549295/1‑70 (MQ=255)
                               gAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCg                                 <  1:3245899/70‑1 (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:5474961/48‑1 (MQ=255)
                                  tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                                                    <  1:2251878/48‑1 (MQ=255)
                                  tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgacg                                             >  1:5491071/1‑55 (MQ=255)
                                  tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgacg                                             >  1:3691850/1‑55 (MQ=255)
                                  tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgac                                              <  1:63102/54‑1 (MQ=255)
                                  tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgac                                              <  1:759894/54‑1 (MQ=255)
                                    cGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAg                            >  1:1058439/1‑70 (MQ=255)
                                    cGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAg                            >  1:4926516/1‑70 (MQ=255)
                                    cGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAg                            >  1:3759086/1‑70 (MQ=255)
                                         gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                         >  1:932022/1‑68 (MQ=255)
                                         gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                         >  1:5591560/1‑68 (MQ=255)
                                         gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                         >  1:4279799/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:5157556/1‑65 (MQ=255)
                                              cGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACg                       >  1:5169788/1‑65 (MQ=255)
                                               atCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAAcc                       <  1:4339776/63‑2 (MQ=255)
                                               atCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAAcc                       <  1:234300/63‑2 (MQ=255)
                                               atCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAAcc                       <  1:5623845/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:5176971/1‑49 (MQ=255)
                                                          tAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGc                           >  1:4153184/1‑49 (MQ=255)
                                                           aaCGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAgcg    >  1:4978297/1‑71 (MQ=255)
                                                             cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGc                                  >  1:3304856/1‑39 (MQ=255)
                                                             cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                         >  1:1099062/1‑48 (MQ=255)
                                                             cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                         >  1:810261/1‑48 (MQ=255)
                                                              gAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCt  >  1:2343918/1‑70 (MQ=255)
                                                              |                                                                     
TCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCT  >  minE/1934126‑1934257

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: