Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,330 G→A 22.7% I371I (ATC→ATT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,3300GA22.7% 249.6 / 67.2 132I371I (ATC→ATTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (55/47);  new base A (15/15);  total (70/62)
Fisher's exact test for biased strand distribution p-value = 8.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 136 total aligned reads displayed.

CCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGA  >  minE/1933260‑1933400
                                                                      |                                                                      
  ggATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAAt                                                                      <  1:2151697/71‑1 (MQ=255)
  ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGAt                                                                      <  1:4125749/71‑1 (MQ=255)
    atatTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATAc                                                                    >  1:5397678/1‑71 (MQ=255)
    atatTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATAc                                                                    <  1:1503976/71‑1 (MQ=255)
    atatTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATAc                                                                    >  1:2914735/1‑71 (MQ=255)
       ttCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGAt                                                                      <  1:3361860/66‑1 (MQ=255)
        tCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGcc                                                                >  1:3587386/1‑71 (MQ=255)
         ccgTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTcc                                                                >  1:228746/3‑70 (MQ=255)
         cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGcc                                                                >  1:1606530/1‑70 (MQ=255)
          cgTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTccc                                                               <  1:3444668/69‑1 (MQ=255)
             aaTGAAGCTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAc                                                           >  1:3738040/1‑71 (MQ=255)
              aTGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcc                                                          >  1:3130815/1‑71 (MQ=255)
              aTGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcc                                                          >  1:5092796/1‑71 (MQ=255)
              aTGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcc                                                          >  1:2666226/1‑71 (MQ=255)
               tGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                       <  1:4639703/57‑1 (MQ=255)
               tGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                       <  1:3930688/57‑1 (MQ=255)
               tGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                       <  1:171847/57‑1 (MQ=255)
               tGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                       <  1:1088435/57‑1 (MQ=255)
               tGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGcc                                                                <  1:3814991/64‑1 (MQ=255)
                 aaGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa                                                                     >  1:4037722/1‑57 (MQ=255)
                 aaGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa                                                                     >  1:450541/1‑57 (MQ=255)
                 aaGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa                                                                     >  1:2208617/1‑57 (MQ=255)
                 aaGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCa                                                              >  1:25551/1‑64 (MQ=255)
                 aaGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc                                                         >  1:1942019/1‑69 (MQ=255)
                 aaGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc                                                         >  1:5525834/1‑69 (MQ=255)
                  aGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAACCCCCAg                                                      <  1:375750/71‑1 (MQ=255)
                     tCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                       >  1:3669043/1‑51 (MQ=255)
                       ggCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAc                                                           <  1:58183/61‑1 (MQ=255)
                       ggCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAc                                                           <  1:4098982/61‑1 (MQ=255)
                       ggCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAc                                                           <  1:2444210/61‑1 (MQ=255)
                       ggCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGctct                                                  <  1:3906577/70‑1 (MQ=255)
                       ggCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAACCCCCa                                                       >  1:3610316/1‑65 (MQ=255)
                        gCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                       <  1:3830996/48‑1 (MQ=255)
                           gAATCAGATGTGCTTCTGCCACTTCTTCGAATATCTCGGTAAAGa                                                                       <  1:3153150/45‑1 (MQ=38)
                           gAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAACCCCCAGcttt                                                  >  1:3166104/1‑64 (MQ=255)
                            aaTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcccc                                                        >  1:1300250/1‑59 (MQ=255)
                            aaTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTtctc                                              <  1:3421519/69‑1 (MQ=255)
                            aaTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCa                                             <  1:3710205/70‑1 (MQ=255)
                            aaTCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAACCCCCAGCTTTtctcta                                            >  1:1322916/1‑69 (MQ=255)
                            aaTCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAACCCCCAGCTTTtctc                                              <  1:1853951/69‑1 (MQ=255)
                            aaTCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAACCCCCAGCTTTtctc                                              <  1:1493536/69‑1 (MQ=255)
                             aTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcc                                                          <  1:3894561/56‑1 (MQ=255)
                             aTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTtctc                                              <  1:4173209/68‑1 (MQ=255)
                              tCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                       <  1:215251/42‑1 (MQ=255)
                              tCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcccc                                                        >  1:2997193/1‑57 (MQ=255)
                               cAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGccc                                                               >  1:3603968/1‑49 (MQ=255)
                               cAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCaa                                            >  1:1590314/1‑68 (MQ=255)
                               cAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAACCCCCAGCTTTtctc                                              >  1:50082/1‑66 (MQ=255)
                                aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGcc                                                                <  1:133440/47‑1 (MQ=255)
                                aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                                        <  1:3297181/71‑1 (MQ=255)
                                aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                                        >  1:4659103/1‑71 (MQ=255)
                                aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                                        <  1:2546909/71‑1 (MQ=255)
                                 gATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTgg                                       <  1:4687913/71‑1 (MQ=255)
                                   tgtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGAt                                     <  1:5369291/71‑1 (MQ=255)
                                   tgtgCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCa                                                              <  1:1279444/46‑1 (MQ=19)
                                    gtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa                                                                     >  1:662941/1‑38 (MQ=255)
                                    gtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa                                                                     >  1:471681/1‑38 (MQ=255)
                                    gtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCa                                                              >  1:1040531/1‑45 (MQ=255)
                                    gtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCtcttc                                                >  1:4011272/1‑59 (MQ=255)
                                    gtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATg                                    >  1:2899991/1‑71 (MQ=255)
                                    gtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATg                                    >  1:2393523/1‑71 (MQ=255)
                                     tgCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAACCCCCAGCTTTTCTCTAcc                                          <  1:3678908/64‑2 (MQ=255)
                                       cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCtcttc                                                >  1:4239722/1‑56 (MQ=255)
                                       cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTgg                                       >  1:347108/1‑65 (MQ=255)
                                       cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGcc                                  >  1:1034932/1‑70 (MQ=255)
                                       cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGcc                                  >  1:1536290/1‑70 (MQ=255)
                                       cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGcc                                  >  1:1875216/1‑70 (MQ=255)
                                       cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGcc                                  >  1:3481753/1‑70 (MQ=255)
                                       cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGcc                                  >  1:2069040/1‑70 (MQ=255)
                                         tCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCa                                                              >  1:2966395/1‑40 (MQ=255)
                                         tCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCa                                                              >  1:1984331/1‑40 (MQ=255)
                                          cTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGCCCCCAGCtcttc                                                >  1:2372027/1‑53 (MQ=255)
                                          cTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGa                                <  1:2313055/69‑1 (MQ=255)
                                            gCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTtctc                                              >  1:3801942/1‑53 (MQ=255)
                                               aCTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCtctt                                                 >  1:5169978/1‑47 (MQ=255)
                                                cttcttCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGAtt                          >  1:5553026/1‑69 (MQ=255)
                                                cttcttCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGAtt                          >  1:548779/1‑69 (MQ=255)
                                                  tcttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAg                       >  1:2793419/1‑70 (MQ=255)
                                                  tcttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAg                       >  1:296842/1‑70 (MQ=255)
                                                   cttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGc                      <  1:697422/70‑1 (MQ=255)
                                                   cttcGAAGATCTAGGTAACGATACTGCCAAGACCCCAGct                                                    >  1:4581491/1‑40 (MQ=255)
                                                    ttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGAt                                     <  1:1923287/54‑1 (MQ=255)
                                                    ttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGAt                                     <  1:2489382/54‑1 (MQ=255)
                                                     tcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGAt                           <  1:2928213/63‑1 (MQ=255)
                                                     tcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAAt                   >  1:1137033/1‑71 (MQ=255)
                                                     tcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAAt                   >  1:693485/1‑71 (MQ=255)
                                                        aaGATCTCTGTGACAATACGTCCCAACCCCCAGCTTTTCTCTACCGTAATACCGATAGATTCAGCTAATGc                <  1:5280438/71‑1 (MQ=255)
                                                         aGATCTCGGTAACGATACGGCCCAGACCCCAGCtctt                                                 >  1:2710813/1‑37 (MQ=255)
                                                          gATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCtt              <  1:4023790/71‑1 (MQ=255)
                                                           aTCTCTGTGACAATACGTCCCAACCCCCAGCTTTTCTCTACCGTAATACCGATAGATTCAGCTAATGCtt              >  1:2622580/1‑70 (MQ=255)
                                                             ctctGTGACAATACGTCCCAACCCCCAGCTTTTCTCTACCGTAATACCGATAGATTCAGCTAATGCttta            >  1:2088103/1‑69 (MQ=255)
                                                                ggTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTc            <  1:476645/67‑1 (MQ=255)
                                                                ggTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTc            <  1:1343296/67‑1 (MQ=255)
                                                                ggTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCg           >  1:2012272/1‑68 (MQ=255)
                                                                 gTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCaga       <  1:5367815/71‑1 (MQ=255)
                                                                  tAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                                        <  1:4947452/37‑1 (MQ=255)
                                                                  tAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCagag      >  1:3408136/1‑71 (MQ=255)
                                                                      gATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGc                      <  1:2944620/51‑1 (MQ=255)
                                                                      gATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGc                      <  1:3530729/51‑1 (MQ=255)
                                                                      gATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGa  <  1:5642619/71‑1 (MQ=255)
                                                                      |                                                                      
CCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCGA  >  minE/1933260‑1933400

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: