Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,206 A→G 21.5% I79I (ATT→ATC lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,2060AG21.5% 135.0 / 34.2 71I79I (ATT→ATClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (23/32);  new base G (8/7);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

TCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTT  >  minE/1934138‑1934276
                                                                    |                                                                      
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:2562071/1‑70 (MQ=255)
             aCAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcat                                                          >  1:3091605/1‑70 (MQ=255)
              cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatc                                                         <  1:1020139/70‑1 (MQ=255)
              cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatc                                                         <  1:1719053/70‑1 (MQ=255)
              cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatc                                                         <  1:146099/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: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:2322979/71‑1 (MQ=255)
                                    tCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTc                                  <  1:2141904/71‑1 (MQ=255)
                                    tCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTc                                  <  1:2844694/71‑1 (MQ=255)
                                      cTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCg                                 <  1:2069430/70‑1 (MQ=255)
                                      cTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCg                                 <  1:1274658/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:598419/1‑71 (MQ=255)
                                          aGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGtt                            >  1:258544/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:3248566/1‑71 (MQ=255)
                                                 cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCt                     <  1:3849325/71‑1 (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)
                                                                ccATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttctt           <  1:3335348/66‑1 (MQ=255)
                                                                 cATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTcaa                                >  1:1452122/1‑42 (MQ=255)
                                                                 cATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTcaa                                >  1:2869328/1‑42 (MQ=255)
                                                                 cATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTc                                  <  1:2128789/42‑1 (MQ=255)
                                                                 cATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttcttc          <  1:1460137/66‑1 (MQ=255)
                                                                   taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct   <  1:170121/71‑1 (MQ=255)
                                                                   taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct   <  1:425590/71‑1 (MQ=255)
                                                                   taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct   <  1:2329692/71‑1 (MQ=255)
                                                                    aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt  >  1:3785940/1‑71 (MQ=255)
                                                                    aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt  >  1:37839/1‑71 (MQ=255)
                                                                    aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt  >  1:3586116/1‑71 (MQ=255)
                                                                    |                                                                      
TCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTT  >  minE/1934138‑1934276

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: