Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,016 A→G 21.0% L143L (TTG→CTG)  lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,0160AG21.0% 173.6 / 41.2 92L143L (TTG→CTG) lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (28/44);  new base G (9/10);  total (38/54)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGA  >  minE/1933947‑1934081
                                                                     |                                                                 
gcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAAc                                                                  >  1:3602980/1‑71 (MQ=255)
 gacTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCt                                                                 <  1:3181762/69‑1 (MQ=255)
 gacTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCt                                                                 <  1:4835933/69‑1 (MQ=255)
 cgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACt                                                                 >  1:1744131/1‑71 (MQ=255)
  gcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTc                                                                <  1:3085453/71‑1 (MQ=255)
   cTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCACCt                                                                 >  1:1213546/1‑69 (MQ=255)
      ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTcag                                                              <  1:1144142/69‑3 (MQ=255)
      ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGta                                                            >  1:5676856/1‑70 (MQ=255)
       cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTca                                                               >  1:4077343/1‑66 (MQ=255)
       cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGtac                                                           <  1:1974970/71‑3 (MQ=255)
       cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGtac                                                           <  1:4136137/71‑3 (MQ=255)
       cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGta                                                            <  1:4003635/70‑2 (MQ=255)
        tGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTg                                                            >  1:285824/1‑69 (MQ=255)
          aTCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                                                        <  1:4825226/71‑1 (MQ=255)
          aTCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGGGCAGt                                                        >  1:1274071/1‑71 (MQ=255)
          aTCCTGCAAGCCGTGGAATTTATACGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGt                                                        >  1:2917538/1‑71 (MQ=255)
           tCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                                                        >  1:4916423/1‑70 (MQ=255)
           tCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                                                        >  1:3631352/1‑70 (MQ=255)
           tCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                                                        >  1:2325026/1‑70 (MQ=255)
           tCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTg                                                       <  1:4174745/71‑1 (MQ=255)
           tCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTg                                                       <  1:3500613/71‑1 (MQ=255)
           tCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTg                                                       <  1:1761591/71‑1 (MQ=255)
           tCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGAGCAGTg                                                       <  1:735382/71‑1 (MQ=255)
                    ccGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                              <  1:4668860/71‑1 (MQ=255)
                     catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                                             >  1:4588118/3‑70 (MQ=255)
                     catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                                             >  1:2573148/3‑70 (MQ=255)
                     catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                                             >  1:1855730/3‑70 (MQ=255)
                     cGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACa                                               <  1:4822578/69‑1 (MQ=255)
                      atGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTg                                                       <  1:1258198/59‑1 (MQ=255)
                       tGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAGCTc                                           <  1:3823357/71‑1 (MQ=255)
                       tGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCg                                                               <  1:2863152/51‑1 (MQ=255)
                        ggAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGa                                                 >  1:2926615/1‑64 (MQ=255)
                         gAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGc                                         <  1:191712/71‑1 (MQ=255)
                          aaTTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTc                                                                <  1:386645/47‑1 (MQ=255)
                            tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACt                                                                 <  1:5386285/44‑1 (MQ=255)
                            tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACt                                                                 <  1:3768983/44‑1 (MQ=255)
                            tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACt                                                                 <  1:2585115/44‑1 (MQ=255)
                            tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTgg                                                      <  1:3767021/55‑1 (MQ=255)
                            tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCgg                                      >  1:2880493/1‑71 (MQ=255)
                             ttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGAc                                                <  1:2890008/60‑1 (MQ=255)
                                  cGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttgg                                <  1:5295937/71‑1 (MQ=255)
                                  cGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttgg                                <  1:3564349/71‑1 (MQ=255)
                                    gCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTc                                           <  1:4178431/58‑1 (MQ=255)
                                    gCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttggtt                              >  1:1499819/1‑71 (MQ=255)
                                    gCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttggtt                              >  1:5415451/1‑71 (MQ=255)
                                      aGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggtttt                                  >  1:4613305/1‑65 (MQ=255)
                                       gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                                                 <  1:1941728/49‑1 (MQ=255)
                                       gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                                                 <  1:4663411/49‑1 (MQ=255)
                                       gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                                                 <  1:861621/49‑1 (MQ=255)
                                       gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGgttttggttttg                           <  1:215619/71‑1 (MQ=255)
                                       gCGGACGCAGTGCTTTGGGCAGCAGACGCAACTCGGTGCAGTGGATAGa                                                 <  1:3394994/49‑1 (MQ=255)
                                        cGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGa                          <  1:5582937/71‑1 (MQ=255)
                                          gACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAAc                        <  1:4641451/71‑1 (MQ=255)
                                           aCGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGaa                         >  1:2462441/1‑69 (MQ=255)
                                              cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                              <  1:3062072/45‑1 (MQ=255)
                                              cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                              <  1:4849127/45‑1 (MQ=255)
                                              cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                              <  1:1403694/45‑1 (MQ=255)
                                              cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                              <  1:3005811/45‑1 (MQ=255)
                                                gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGAc                                                <  1:3277571/41‑1 (MQ=255)
                                                gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGt                                             <  1:5187435/44‑1 (MQ=255)
                                                gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGt                                             <  1:5414204/44‑1 (MQ=255)
                                                gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGc                     <  1:5055446/68‑1 (MQ=255)
                                                gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGc                     <  1:1988341/68‑1 (MQ=255)
                                                 tGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggtt                                    <  1:2294863/52‑1 (MQ=255)
                                                  gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc                >  1:4893703/1‑71 (MQ=255)
                                                  gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc                >  1:854649/1‑71 (MQ=255)
                                                  gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc                >  1:47924/1‑71 (MQ=255)
                                                  gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc                >  1:1074442/1‑71 (MQ=255)
                                                  gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc                >  1:4105097/1‑71 (MQ=255)
                                                  gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAc                 >  1:1799409/1‑70 (MQ=255)
                                                   cTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                              >  1:1796875/1‑40 (MQ=255)
                                                   cTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                              >  1:5202966/1‑40 (MQ=255)
                                                   cTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCt               >  1:5323365/1‑71 (MQ=255)
                                                     ttGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                                                 <  1:2630518/35‑1 (MQ=255)
                                                      tGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACa                       >  1:591886/1‑60 (MQ=255)
                                                      tGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCttt             <  1:2978676/70‑1 (MQ=255)
                                                      tGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCtt              <  1:4657059/69‑1 (MQ=255)
                                                      tGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTc            >  1:2457260/1‑71 (MQ=255)
                                                      tGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTc            >  1:5371789/1‑71 (MQ=255)
                                                       ggTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTcg           >  1:1529964/1‑71 (MQ=255)
                                                           agcagACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGcc       >  1:5550570/1‑71 (MQ=255)
                                                           agcagACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGcc       <  1:5575288/71‑1 (MQ=255)
                                                             cagTCGCAACTCGGTGCAGTGGATAGACAGTTCGCCGgttttggttttg                           <  1:848158/49‑1 (MQ=255)
                                                             cagACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGgttttggttttg                           <  1:969624/49‑1 (MQ=255)
                                                             cagACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCtt                   <  1:4008047/57‑1 (MQ=255)
                                                              agACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg    >  1:3497594/1‑71 (MQ=255)
                                                              agACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg    <  1:4847325/71‑1 (MQ=255)
                                                              agACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg    >  1:5254303/1‑71 (MQ=255)
                                                              agACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg    >  1:4204476/1‑71 (MQ=255)
                                                               gACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttt                                   <  1:3367423/39‑1 (MQ=255)
                                                               gACGCAACTCGGTGCAGTGGATAGACAGTTAGCCggttt                                   <  1:4013623/39‑1 (MQ=255)
                                                                  gCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGa  >  1:1397935/1‑69 (MQ=255)
                                                                     |                                                                 
GCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGA  >  minE/1933947‑1934081

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: