Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,188 G→A 11.2% F85F (TTC→TTT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,1880GA11.2% 271.1 / 18.8 98F85F (TTC→TTTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (42/45);  new base A (5/6);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCT  >  minE/1934119‑1934257
                                                                     |                                                                     
aaCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      <  1:414129/71‑1 (MQ=255)
aaCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      <  1:1110733/71‑1 (MQ=255)
aaCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      <  1:1251231/71‑1 (MQ=255)
 aCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      >  1:2757308/1‑70 (MQ=255)
  cGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      >  1:3482992/1‑69 (MQ=255)
  cGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      >  1:2351674/1‑69 (MQ=255)
  cGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      >  1:1195468/1‑69 (MQ=255)
     cTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGAc                                                                 >  1:3822468/1‑71 (MQ=255)
        ccGGGAGATCGTCACGGGCAACGTACAGCTGTATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGc                                                               <  1:1339389/70‑1 (MQ=255)
        ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACg                                                                <  1:1693437/69‑1 (MQ=255)
        ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCt                                                              >  1:1586997/1‑71 (MQ=255)
          gggAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      <  1:3897061/61‑1 (MQ=255)
          gggAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCtt                                                             >  1:2888995/1‑70 (MQ=255)
              gaTCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTAcc                                                         >  1:3931835/1‑70 (MQ=255)
              gaTCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTAcc                                                         >  1:1706607/1‑70 (MQ=255)
              gaTCGTCACGGGCAACGTACAGCTGAAAGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTAccc                                                        <  1:436605/71‑1 (MQ=255)
               aTCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTAccc                                                        <  1:2413681/70‑1 (MQ=255)
               aTCGTCACGGGCAACGTACAGCTGAATACGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCa                                                       >  1:3835260/1‑71 (MQ=255)
                tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                                                    <  1:3817657/57‑1 (MQ=255)
                tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                                                    <  1:1893285/57‑1 (MQ=255)
                tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                                                    <  1:1682024/57‑1 (MQ=255)
                tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                      <  1:654561/71‑1 (MQ=255)
                tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                      <  1:46477/71‑1 (MQ=255)
                tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                      <  1:3872507/71‑1 (MQ=255)
                 cGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                      >  1:1019782/1‑70 (MQ=255)
                 cGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                      >  1:2058786/1‑70 (MQ=255)
                 cGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                      >  1:1667779/1‑70 (MQ=255)
                  gTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTAccc                                                        >  1:1479895/1‑67 (MQ=255)
                   tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                                                      <  1:2197201/52‑1 (MQ=255)
                   tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                                                    >  1:2383178/1‑70 (MQ=255)
                   tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                                                    >  1:874430/1‑70 (MQ=255)
                    cACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaa                                                                     >  1:3636070/1‑52 (MQ=255)
                    cACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaa                                                                     >  1:3682362/1‑52 (MQ=255)
                     aCGTGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                                 >  1:1721183/1‑71 (MQ=255)
                     aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                                                    >  1:314147/1‑52 (MQ=255)
                     aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                                 <  1:3568748/71‑1 (MQ=255)
                     aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                                 <  1:2946840/71‑1 (MQ=255)
                      cttgCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATac                                                 >  1:1342083/4‑70 (MQ=255)
                       gggCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCata                                                     >  1:2664461/1‑65 (MQ=255)
                       gggCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATTATacga                                               <  1:3212337/71‑1 (MQ=255)
                         gCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgacg                                             <  1:1001202/71‑1 (MQ=255)
                          cAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                                                    <  1:1176573/47‑1 (MQ=255)
                             cGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGtc                                         <  1:3945102/71‑1 (MQ=255)
                             cGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCt                                                              >  1:2175408/1‑50 (MQ=255)
                                aCAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGc                                                               >  1:1685103/1‑46 (MQ=255)
                                aCAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcat                                       >  1:3091605/1‑70 (MQ=255)
                                aCAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcat                                       >  1:2562071/1‑70 (MQ=255)
                                 cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatc                                      <  1:1719053/70‑1 (MQ=255)
                                 cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatc                                      <  1:146099/70‑1 (MQ=255)
                                 cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatc                                      <  1:1020139/70‑1 (MQ=255)
                                  aGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCTATAATACGACGGGTcatcat                                    <  1:1842042/71‑1 (MQ=255)
                                    ttGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTcatcattc                                  >  1:1474238/2‑69 (MQ=255)
                                    ttGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTcatcattc                                  >  1:2098882/2‑69 (MQ=255)
                                    ttGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTcatcat                                    <  1:3442194/68‑1 (MQ=255)
                                    cTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACggg                                           <  1:2773660/62‑1 (MQ=255)
                                    cTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATg                                   <  1:3424822/70‑1 (MQ=255)
                                    cTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATg                                   <  1:3422759/70‑1 (MQ=255)
                                     tGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgac                                              <  1:2419705/58‑1 (MQ=255)
                                      gAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCgg                                >  1:613890/1‑71 (MQ=255)
                                      gAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCgg                                >  1:2435959/1‑71 (MQ=255)
                                       aaTGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCgg                                <  1:1966525/70‑1 (MQ=255)
                                       aaTACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCgg                                <  1:1549006/70‑1 (MQ=255)
                                        aTGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGc                               <  1:3845564/70‑1 (MQ=255)
                                             aCCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCaa                         <  1:1676818/71‑1 (MQ=255)
                                               cGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACg                       <  1:3731146/71‑1 (MQ=255)
                                                gcaaaCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACgg                      >  1:1476048/4‑71 (MQ=255)
                                                gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAAAgg                      >  1:1291757/1‑71 (MQ=255)
                                                gCCAACGTCCTGCAGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGc                     <  1:1437688/71‑1 (MQ=255)
                                                 ccAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataata                                                  <  1:3426055/42‑1 (MQ=255)
                                                 ccAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGc                     <  1:2043395/71‑1 (MQ=255)
                                                 ccAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGc                     >  1:3403473/1‑71 (MQ=255)
                                                   aaCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCa                             >  1:1720525/1‑61 (MQ=255)
                                                    aCGTCCTGCAGGGTAACGAAAGACGCTTTACCCataata                                                  >  1:1275403/1‑39 (MQ=255)
                                                    aCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGc                                  >  1:2494688/1‑55 (MQ=255)
                                                    aCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGc                                  >  1:3071995/1‑55 (MQ=255)
                                                    aCGTCCTGCAGGGTAACGAAAGACGATTTACCCataata                                                  >  1:879284/1‑39 (MQ=255)
                                                       tCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTc               <  1:2141904/71‑1 (MQ=255)
                                                       tCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTc               <  1:2322979/71‑1 (MQ=255)
                                                       tCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTc               <  1:2844694/71‑1 (MQ=255)
                                                         cTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCg              <  1:1274658/70‑1 (MQ=255)
                                                         cTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCg              <  1:2069430/70‑1 (MQ=255)
                                                         cTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGa             <  1:2120878/71‑1 (MQ=255)
                                                         cTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTca              <  1:3842180/70‑2 (MQ=255)
                                                         cTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTca              <  1:1776962/70‑2 (MQ=255)
                                                         cTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACAGAGACTTcaa             >  1:1040718/1‑69 (MQ=255)
                                                          tGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACtt                <  1:1167715/67‑1 (MQ=255)
                                                            cAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTcatcat                                    >  1:327039/1‑45 (MQ=255)
                                                             aGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGtt         >  1:258544/1‑71 (MQ=255)
                                                             aGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGtt         >  1:598419/1‑71 (MQ=255)
                                                               ggTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCa       <  1:3160277/71‑1 (MQ=255)
                                                                 tAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCg              <  1:1337919/62‑1 (MQ=255)
                                                                  aaCGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACg                       <  1:78977/52‑1 (MQ=255)
                                                                  aaCGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAgcg    >  1:3533358/1‑71 (MQ=255)
                                                                   aCGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGAt            >  1:3665650/1‑62 (MQ=255)
                                                                    cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGtt         >  1:696461/1‑64 (MQ=255)
                                                                    cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCa       >  1:3529402/1‑66 (MQ=255)
                                                                    cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCt  <  1:3849325/71‑1 (MQ=255)
                                                                    cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCt  >  1:3248566/1‑71 (MQ=255)
                                                                     |                                                                     
AACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCT  >  minE/1934119‑1934257

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: