Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,011 A→G 22.4% R144R (CGT→CGC lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,0110AG22.4% 168.9 / 46.6 95R144R (CGT→CGClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (29/44);  new base G (11/10);  total (40/55)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

GCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAG  >  minE/1933945‑1934079
                                                                  |                                                                    
gcgcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGc                                                                   <  1:3250009/70‑1 (MQ=255)
  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)
      ttCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGgcgc                                                                   <  1:3471179/64‑1 (MQ=255)
      ttCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGgcgc                                                                   <  1:4853049/64‑1 (MQ=255)
        ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTcag                                                            <  1:1144142/69‑3 (MQ=255)
        ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGta                                                          >  1:5676856/1‑70 (MQ=255)
        ccTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGAc                                                                     >  1:94270/1‑60 (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)
          tGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCa                                                                  >  1:2312386/1‑61 (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:1855730/3‑70 (MQ=255)
                       catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                                           >  1:2573148/3‑70 (MQ=255)
                       cGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACa                                             <  1:4822578/69‑1 (MQ=255)
                        atGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTg                                                     <  1:1258198/59‑1 (MQ=255)
                         tGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGgcgc                                                                   <  1:3170555/45‑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)
                                tATCCCGCAGCGGACGCAGTGCTTTGGTCAGCAGACg                                                                    >  1:4659781/1‑37 (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:5415451/1‑71 (MQ=255)
                                      gCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttggtt                            >  1:1499819/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:1403694/45‑1 (MQ=255)
                                                cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                            <  1:3005811/45‑1 (MQ=255)
                                                cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                            <  1:4849127/45‑1 (MQ=255)
                                                cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                            <  1:3062072/45‑1 (MQ=255)
                                                  gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGAc                                              <  1:3277571/41‑1 (MQ=255)
                                                  gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGt                                           <  1:5414204/44‑1 (MQ=255)
                                                  gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGt                                           <  1:5187435/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:4105097/1‑71 (MQ=255)
                                                    gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc              >  1:4893703/1‑71 (MQ=255)
                                                    gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc              >  1:47924/1‑71 (MQ=255)
                                                    gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc              >  1:1074442/1‑71 (MQ=255)
                                                    gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc              >  1:854649/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:5371789/1‑71 (MQ=255)
                                                        tGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTc          >  1:2457260/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:5254303/1‑71 (MQ=255)
                                                                agACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg  <  1:4847325/71‑1 (MQ=255)
                                                                agACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg  >  1:4204476/1‑71 (MQ=255)
                                                                agACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg  >  1:3497594/1‑71 (MQ=255)
                                                                 gACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttt                                 <  1:3367423/39‑1 (MQ=255)
                                                                 gACGCAACTCGGTGCAGTGGATAGACAGTTAGCCggttt                                 <  1:4013623/39‑1 (MQ=255)
                                                                  |                                                                    
GCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAG  >  minE/1933945‑1934079

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: