Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,639 G→A 14.7% S268S (TCC→TCT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,6390GA14.7% 276.5 / 32.5 109S268S (TCC→TCTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (50/43);  new base A (9/7);  total (59/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 = 5.39e-01

Only 100 of 116 total aligned reads displayed.

GATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCT  >  minE/1933570‑1933709
                                                                     |                                                                      
gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACgg                                                                       <  1:3446948/71‑1 (MQ=255)
gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACg                                                                        >  1:3189364/1‑70 (MQ=255)
gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACg                                                                        >  1:3195083/1‑70 (MQ=255)
gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACg                                                                        <  1:1206485/70‑1 (MQ=255)
gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACg                                                                        >  1:1757474/1‑70 (MQ=255)
 aTCAGATCTTTGTAATCTGCGTAAGCCATGAAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                                                                      <  1:1078116/71‑1 (MQ=255)
 aTCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGa                                                                      <  1:3543291/71‑1 (MQ=255)
  tCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaa                                                                     <  1:666948/71‑1 (MQ=255)
  tCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                                                                      >  1:3618664/1‑70 (MQ=255)
   cAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTGCATCATGGTGAACTCTGGGTTATGACGTACGGaa                                                                     <  1:1139601/70‑1 (MQ=255)
    aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                                                                    >  1:1231389/1‑70 (MQ=255)
      aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                                                                      <  1:2270945/66‑1 (MQ=255)
        cTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                                                                    >  1:3423575/1‑66 (MQ=255)
         tttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                                                                      <  1:2408485/63‑1 (MQ=255)
         tttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCtt                                                              >  1:1431638/1‑71 (MQ=255)
         tttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCtt                                                              >  1:1810435/1‑71 (MQ=255)
          ttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTc                                                             >  1:3648812/1‑71 (MQ=255)
          ttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTc                                                             >  1:3713919/1‑71 (MQ=255)
           tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCg                                                            >  1:2566959/1‑71 (MQ=255)
           tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCg                                                            >  1:3843048/1‑71 (MQ=255)
             tAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGtt                                                          <  1:3654142/71‑1 (MQ=255)
               aTCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTAc                                                        >  1:1800792/1‑71 (MQ=255)
                tCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGtt                                                          <  1:1603635/68‑1 (MQ=255)
                  tGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACgg                                                      >  1:1432986/1‑70 (MQ=255)
                    cGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACgg                                                      <  1:3763661/68‑1 (MQ=255)
                    cGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAg                                                   <  1:1151388/71‑1 (MQ=255)
                      tAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACgg                                                      >  1:1722944/1‑66 (MQ=255)
                         gCCGTGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACg                                              >  1:898899/1‑71 (MQ=255)
                         gCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTa                                                <  1:3566231/69‑1 (MQ=255)
                          ccATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACgg                                             <  1:31784/71‑1 (MQ=255)
                           cATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGt                                            >  1:2432281/1‑71 (MQ=255)
                           cATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGt                                            <  1:606162/71‑1 (MQ=255)
                            aTGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                                                                      <  1:2149913/44‑1 (MQ=255)
                            aTGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGtt                                           <  1:3346568/71‑1 (MQ=255)
                            aTGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGtt                                           <  1:3031687/71‑1 (MQ=255)
                             tGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                                                                      <  1:1944385/43‑1 (MQ=255)
                             tGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGt                                            >  1:883771/1‑69 (MQ=255)
                              gTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTAc                                               <  1:272517/65‑1 (MQ=255)
                               tAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGa                                                                      <  1:3144676/41‑1 (MQ=255)
                               tAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGAt                                        <  1:351408/71‑1 (MQ=255)
                                agagTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACTGAAGTTACGGTTGAtt                                       >  1:2010645/1‑71 (MQ=255)
                                agagTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGAtt                                       >  1:3744116/1‑71 (MQ=255)
                                 gagTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGAttt                                      >  1:2922537/1‑71 (MQ=255)
                                 gagTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGAtt                                       <  1:2123761/70‑1 (MQ=255)
                                 gagTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGAtt                                       <  1:3854885/70‑1 (MQ=255)
                                 gagTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGAtt                                       <  1:3868686/70‑1 (MQ=255)
                                  agTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACgg                                             <  1:2750015/63‑1 (MQ=255)
                                  agTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACgg                                             <  1:2640157/63‑1 (MQ=255)
                                  agTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTc                                     >  1:1255903/1‑71 (MQ=255)
                                  agTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTc                                     <  1:2960355/71‑1 (MQ=255)
                                  agTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTc                                     <  1:2596932/71‑1 (MQ=255)
                                   gTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTc                                     >  1:1558291/1‑70 (MQ=255)
                                   gTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGAt                                        >  1:1976200/1‑67 (MQ=255)
                                   gTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCg                                    <  1:922304/71‑1 (MQ=255)
                                   gTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCg                                    <  1:2235807/71‑1 (MQ=255)
                                   gTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCg                                    <  1:3881830/71‑1 (MQ=255)
                                    ttCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGtt                                           >  1:3804349/1‑63 (MQ=255)
                                    ttCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCg                                    >  1:1794347/1‑70 (MQ=255)
                                    ttCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCg                                    >  1:3181857/1‑70 (MQ=255)
                                     tCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGaa                                  >  1:466830/1‑71 (MQ=255)
                                      ccATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTAc                                                        >  1:92746/1‑48 (MQ=255)
                                      ccATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGaa                                  <  1:3432322/70‑1 (MQ=255)
                                        atcatTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTAccg                                             >  1:1370776/1‑55 (MQ=255)
                                        atcatGGTGAACTCTGGGTTATGACGTACGGAAATACCTTc                                                             >  1:2138235/1‑41 (MQ=255)
                                        atcatGGTGAACTCTGGGTTATGACGTACGGAAATACCTTc                                                             >  1:767046/1‑41 (MQ=255)
                                         tcatGGTGAACTCTGGGTTATGACGTACGGAAATACCtt                                                              <  1:3080424/39‑1 (MQ=255)
                                          catTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATaccc                             >  1:1453502/1‑69 (MQ=255)
                                          catGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGaa                                  >  1:442877/1‑66 (MQ=255)
                                           atGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTa                                                         >  1:2546357/1‑42 (MQ=255)
                                            ttgTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCgt                           <  1:2639753/69‑2 (MQ=255)
                                            tGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGtt                                           >  1:1989752/1‑55 (MQ=255)
                                              gTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTc                         >  1:1442235/1‑71 (MQ=255)
                                               tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGaa                                  >  1:3324809/1‑61 (MQ=255)
                                               tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                        <  1:3220104/71‑1 (MQ=255)
                                               tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                        <  1:2740201/71‑1 (MQ=255)
                                               tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                        <  1:249991/71‑1 (MQ=255)
                                                 aaCTCAGGATTATGGCGAAGAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGaa                                  >  1:1545477/1‑59 (MQ=255)
                                                   ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                    <  1:843694/71‑1 (MQ=255)
                                                   ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                    >  1:3670150/1‑71 (MQ=255)
                                                   ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                    >  1:3504748/1‑71 (MQ=255)
                                                     ctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGcca                  >  1:1724942/1‑71 (MQ=255)
                                                      tGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGccac                 >  1:3473895/1‑71 (MQ=255)
                                                         attATGGCGAACAGAAATACCTTCATTACGGAAGtt                                                 <  1:1261407/35‑1 (MQ=16)
                                                          ttATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCgcctac             <  1:922462/71‑1 (MQ=255)
                                                          ttATGACGTACGGAAATACCTTCGTTAGGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCaccaac             <  1:1063126/71‑1 (MQ=255)
                                                           tATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGa                                   >  1:374608/1‑48 (MQ=37)
                                                           tATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATac                               >  1:2349790/1‑52 (MQ=255)
                                                             tGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac          <  1:1300805/71‑1 (MQ=255)
                                                              gACGTACGGATATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                      >  1:1279412/1‑58 (MQ=255)
                                                              gACGTACGGATATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                      >  1:1274816/1‑58 (MQ=255)
                                                               aCGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccacc         <  1:3086059/70‑1 (MQ=255)
                                                                cGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac          >  1:1376230/1‑68 (MQ=255)
                                                                cGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac          >  1:3630647/1‑68 (MQ=255)
                                                                 gTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGcca                  >  1:189450/1‑59 (MQ=255)
                                                                  tACGGAAATACCTTCGTTACGGAAGTTACGGTTGAtt                                       >  1:1625328/1‑37 (MQ=255)
                                                                  tACGGAAATACCTTCGTTACGGAAGTTACGGTTGAtt                                       >  1:3913098/1‑37 (MQ=255)
                                                                   aCGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACg    >  1:3649797/1‑71 (MQ=255)
                                                                    cGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGc   >  1:2166419/1‑71 (MQ=255)
                                                                     ggAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCt  >  1:1701368/1‑71 (MQ=255)
                                                                     ggAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCt  >  1:2088152/1‑71 (MQ=255)
                                                                     |                                                                      
GATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCT  >  minE/1933570‑1933709

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: