Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,224 G→T 19.9% R73R (CGC→CGA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,2240GT19.9% 206.0 / 27.6 80R73R (CGC→CGAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (31/33);  new base T (5/11);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

CTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGC  >  minE/1934155‑1934292
                                                                     |                                                                    
cTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATg                                                                      <  1:3424822/70‑1 (MQ=255)
cTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATg                                                                      <  1:3422759/70‑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)
             ccAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGc                                                        >  1:3403473/1‑71 (MQ=255)
             ccAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGc                                                        <  1:2043395/71‑1 (MQ=255)
               aaCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCa                                                                >  1:1720525/1‑61 (MQ=255)
                aCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGc                                                                     >  1:2494688/1‑55 (MQ=255)
                aCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGc                                                                     >  1:3071995/1‑55 (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: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)
                         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: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:2433148/1‑50 (MQ=255)
                                           ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGAt                                               >  1:2104523/1‑50 (MQ=255)
                                           ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttcttc                          >  1:3348784/1‑71 (MQ=255)
                                               ccATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttctt                           <  1:3335348/66‑1 (MQ=255)
                                                cATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTcaa                                                >  1:2869328/1‑42 (MQ=255)
                                                cATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTcaa                                                >  1:1452122/1‑42 (MQ=255)
                                                cATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTc                                                  <  1:2128789/42‑1 (MQ=255)
                                                cATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttcttc                          <  1:1460137/66‑1 (MQ=255)
                                                  tgataCGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGTTCCTGGTtatcc                   <  1:2818886/69‑2 (MQ=255)
                                                  taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                   <  1:425590/71‑1 (MQ=255)
                                                  taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                   <  1:2329692/71‑1 (MQ=255)
                                                  taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                   <  1:170121/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)
                                                    ataCGACGGGTCATCATTCGGCCAGCAACCGAGCCTTCAATGTTTAAGGATTCCAGTTCCTGGTtatcct                  <  1:860711/70‑3 (MQ=255)
                                                    ataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCa                                          >  1:372350/1‑46 (MQ=255)
                                                     taCGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGTTCCTGGTTATCctt                 <  1:2698314/70‑1 (MQ=255)
                                                     taCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                  <  1:190178/69‑1 (MQ=255)
                                                      acgacgGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGttcct                           <  1:27264/59‑3 (MQ=255)
                                                      acgacgGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGttcct                           <  1:461225/59‑3 (MQ=255)
                                                        gacgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTcc                                  <  1:2382506/50‑1 (MQ=255)
                                                         acgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                  <  1:760319/65‑1 (MQ=255)
                                                         acgGGTCATCATGCGGCCAGCAACGGAGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                  <  1:2024128/65‑1 (MQ=255)
                                                          cgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCttt                 >  1:3508625/1‑65 (MQ=255)
                                                           gggTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCg          <  1:3623730/71‑1 (MQ=255)
                                                             gTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa        <  1:1904250/71‑1 (MQ=255)
                                                              tcatcaTGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTtctc                    >  1:509796/1‑58 (MQ=255)
                                                              tcatcaTGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt       >  1:2079312/1‑71 (MQ=255)
                                                               catcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGt                        >  1:3870703/1‑53 (MQ=255)
                                                               catcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt       >  1:2449958/1‑70 (MQ=255)
                                                                atcatTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGttc                             <  1:1635228/47‑1 (MQ=255)
                                                                atcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                   <  1:2672363/57‑1 (MQ=255)
                                                                  catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCttt                 >  1:1086959/1‑57 (MQ=255)
                                                                  catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCt    <  1:3895473/70‑1 (MQ=255)
                                                                  catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTg   >  1:1004369/1‑71 (MQ=255)
                                                                   atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGc  >  1:3580720/1‑71 (MQ=255)
                                                                    tGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTcc                                  <  1:475559/38‑1 (MQ=255)
                                                                     |                                                                    
CTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGC  >  minE/1934155‑1934292

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: