Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,167 G→A 22.5% G92G (GGC→GGT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,1670GA22.5% 189.4 / 51.3 102G92G (GGC→GGTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (37/42);  new base A (19/4);  total (56/46)
Fisher's exact test for biased strand distribution p-value = 3.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 107 total aligned reads displayed.

ATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAA  >  minE/1934097‑1934234
                                                                      |                                                                   
aTTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCg                                                                     >  1:2809910/1‑71 (MQ=255)
aTTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCg                                                                     >  1:2943895/1‑71 (MQ=255)
  tttttGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGc                                                                    >  1:418436/1‑70 (MQ=255)
  tttttAAACTGATCGTTATAAACACCTTCTGGCAGGCTATCTCTTGCAACGTACAGTTGAATACGGCCAcc                                                                   >  1:530027/1‑71 (MQ=255)
   ttttAAACTGATCGTTATAAACACCTTCTGGCAGGCTATCTCTTGCAACGTACAGTTGAATACGGCCAcc                                                                   >  1:3316314/1‑70 (MQ=255)
   ttttAAACTGATCGTTATAAACACCTTCTGGCAGGCTATCTCTTGCAACGTACAGTTGAATACGGCCAcc                                                                   >  1:3243048/1‑70 (MQ=255)
    tttGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCa                                                                  <  1:2318093/70‑1 (MQ=255)
      taaaCTGATCGTTATAAACACCTTCTGGCAGGCTATCTCTTGCAACGTACAGTTGAATACGGCCACCGAca                                                               >  1:3441751/3‑70 (MQ=255)
      taaaCTGATCGTTATAAACACCTTCTGGCAGGCTATCTCTTGCAACGTACAGTTGAATACGGCCACCGAca                                                               >  1:2601234/3‑70 (MQ=255)
       gAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACg                                                               <  1:50296/70‑1 (MQ=255)
       gAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACg                                                               <  1:1604905/70‑1 (MQ=255)
        aaCTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTc                                                             <  1:2636232/71‑1 (MQ=255)
         aCTGATCGTTATAAACACCTTCTGGCAGGCTATCTCTTGCAACGTACAGTTGAATACGGCCACCGACATcc                                                            >  1:3362005/1‑71 (MQ=255)
             cTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTg                                                          <  1:2561669/69‑1 (MQ=255)
             cTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCa                                                        <  1:3604473/71‑1 (MQ=255)
               cGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAgg                                                      >  1:4000042/1‑71 (MQ=255)
                gTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAggg                                                     <  1:1469185/71‑1 (MQ=255)
                gTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAggg                                                     >  1:3317548/1‑71 (MQ=255)
                  tataAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGc                                                         <  1:4041612/65‑1 (MQ=255)
                  tataAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGt                                                    >  1:218029/1‑70 (MQ=255)
                  tataAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGt                                                    >  1:4648808/1‑70 (MQ=255)
                   ataAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTaa                                                  >  1:4515001/1‑71 (MQ=255)
                   ataAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTaa                                                  >  1:3361187/1‑71 (MQ=255)
                      aaCGCCTTCCGGGAGATCGTCACGGGCAACGTCCAGCTGAATGCGACCGCCa                                                                  >  1:1493144/1‑52 (MQ=255)
                        cGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTc                                                             >  1:1558446/1‑55 (MQ=255)
                            ttCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGcc                                                                   <  1:1479649/45‑1 (MQ=255)
                            ttCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGc                                                         <  1:3751635/55‑1 (MQ=255)
                             tCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGc                                        >  1:2419192/1‑71 (MQ=255)
                              ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTTCAGGGTAACGAAAGACGc                                        <  1:4037873/70‑1 (MQ=255)
                              ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCt                                       >  1:3594207/1‑71 (MQ=255)
                              ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCt                                       >  1:43532/1‑71 (MQ=255)
                              ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCt                                       >  1:830366/1‑71 (MQ=255)
                                gggAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCttt                                     <  1:2600185/71‑1 (MQ=255)
                                 ggAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTa                                    <  1:4818060/71‑1 (MQ=255)
                                    gaTCGTCACGGGCAACGTACAGCTGAATGCGACCGcc                                                                   >  1:3085563/1‑37 (MQ=255)
                                    gaTCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAAc                                                                >  1:5595980/1‑40 (MQ=255)
                                      tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                               <  1:1369846/71‑1 (MQ=255)
                                      tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                               <  1:3923450/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)
                                          cACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAgg                                                      >  1:2967745/1‑44 (MQ=255)
                                          cACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAgg                                                      >  1:5028445/1‑44 (MQ=255)
                                          cACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGa                                           >  1:3729867/1‑55 (MQ=255)
                                           aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                          >  1:5192515/1‑71 (MQ=255)
                                           aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                          >  1:4089260/1‑71 (MQ=255)
                                           aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                          >  1:5635555/1‑71 (MQ=255)
                                           aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                          >  1:2305793/1‑71 (MQ=255)
                                           aCGGGCAACGTACAGATGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                             >  1:2103692/1‑68 (MQ=255)
                                            cttgCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATac                          >  1:4010170/4‑70 (MQ=255)
                                            cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                               <  1:5608964/49‑1 (MQ=255)
                                            cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                               <  1:1083991/49‑1 (MQ=255)
                                            cGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                               <  1:3553945/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)
                                               gCAACGTACAGTTGAATCCGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTccc                                  >  1:2586352/1‑59 (MQ=255)
                                               gCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTccc                                  >  1:3486407/1‑59 (MQ=255)
                                               gCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTccc                                  >  1:2435601/1‑59 (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)
                                                       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: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: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:63102/54‑1 (MQ=255)
                                                               tGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgac                       <  1:759894/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:5591560/1‑68 (MQ=255)
                                                                      gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa  >  1:4279799/1‑68 (MQ=255)
                                                                      gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa  >  1:932022/1‑68 (MQ=255)
                                                                      |                                                                   
ATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAA  >  minE/1934097‑1934234

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: