Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,200 A→C 15.5% G81G (GGT→GGG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,2000AC15.5% 161.1 / 33.8 84G81G (GGT→GGGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base C (6/7);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

GATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTC  >  minE/1934133‑1934268
                                                                   |                                                                    
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)
  tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                                 <  1:3872507/71‑1 (MQ=255)
  tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                                 <  1:654561/71‑1 (MQ=255)
  tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                                                                 <  1:46477/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)
     tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                                                               >  1:874430/1‑70 (MQ=255)
     tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                                                               >  1:2383178/1‑70 (MQ=255)
       aCGTGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                                            >  1:1721183/1‑71 (MQ=255)
       aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                                            <  1:2946840/71‑1 (MQ=255)
       aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                                                            <  1:3568748/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)
               cGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGtc                                                    <  1:3945102/71‑1 (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:1020139/70‑1 (MQ=255)
                   cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatc                                                 <  1:146099/70‑1 (MQ=255)
                    aGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCTATAATACGACGGGTcatcat                                               <  1:1842042/71‑1 (MQ=255)
                      ttGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTcatcattc                                             >  1:2098882/2‑69 (MQ=255)
                      ttGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTcatcattc                                             >  1:1474238/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:3403473/1‑71 (MQ=255)
                                   ccAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGc                                <  1:2043395/71‑1 (MQ=255)
                                     aaCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCa                                        >  1:1720525/1‑61 (MQ=255)
                                      aCGTCCTGCAGGGTAACGAAAGACGCTTTACCCataata                                                             >  1:1275403/1‑39 (MQ=255)
                                      aCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGc                                             >  1:3071995/1‑55 (MQ=255)
                                      aCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGc                                             >  1:2494688/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:1776962/70‑2 (MQ=255)
                                           cTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTca                         <  1:3842180/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)
                                                         aaGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGAc                             >  1:1816200/1‑52 (MQ=255)
                                                         aaGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTc           >  1:66336/1‑70 (MQ=255)
                                                           gaggCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAg        >  1:2529101/4‑71 (MQ=255)
                                                              gCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGttc     <  1:3302059/71‑1 (MQ=255)
                                                              gCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttc     <  1:3534421/71‑1 (MQ=255)
                                                               cTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttct    <  1:3938307/71‑1 (MQ=255)
                                                                 ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGAt                       >  1:2104523/1‑50 (MQ=255)
                                                                 ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGAt                       >  1:2433148/1‑50 (MQ=255)
                                                                 ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttcttc  >  1:3348784/1‑71 (MQ=255)
                                                                   |                                                                    
GATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTC  >  minE/1934133‑1934268

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: