Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,191 A→G 24.7% S84S (TCT→TCC lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,1910AG24.7% 145.0 / 50.9 89S84S (TCT→TCClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (38/29);  new base G (13/9);  total (51/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

TCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCC  >  minE/1934126‑1934260
                                                                 |                                                                     
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:1506555/68‑1 (MQ=255)
            tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                         <  1:1500069/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:3670634/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)
               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:3486407/1‑59 (MQ=255)
                  gCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTccc                                                            >  1:2435601/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)
                       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)
                          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: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:5491071/1‑55 (MQ=255)
                                  tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgacg                                                >  1:3691850/1‑55 (MQ=255)
                                  tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgac                                                 <  1:759894/54‑1 (MQ=255)
                                  tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgac                                                 <  1:63102/54‑1 (MQ=255)
                                    cGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAg                               >  1:4926516/1‑70 (MQ=255)
                                    cGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAg                               >  1:3759086/1‑70 (MQ=255)
                                    cGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAg                               >  1:1058439/1‑70 (MQ=255)
                                         gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                            >  1:4279799/1‑68 (MQ=255)
                                         gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                            >  1:5591560/1‑68 (MQ=255)
                                         gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                            >  1:932022/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:234300/63‑2 (MQ=255)
                                               atCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAAcc                          <  1:5623845/63‑2 (MQ=255)
                                               atCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAAcc                          <  1:4339776/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)
                                                             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)
                                                               aaaGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTc   >  1:3316697/1‑71 (MQ=255)
                                                                aaGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTcc  <  1:2878285/71‑1 (MQ=255)
                                                                 |                                                                     
TCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCC  >  minE/1934126‑1934260

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: