Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,681 A→C 18.2% R254R (CGT→CGG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,6810AC18.2% 126.8 / 40.8 77R254R (CGT→CGGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (35/28);  new base C (12/2);  total (47/30)
Fisher's exact test for biased strand distribution p-value = 6.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

ATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGA  >  minE/1933613‑1933748
                                                                    |                                                                   
atGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATacac                                                                    <  1:3805667/70‑1 (MQ=255)
  ggTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGc                                                                  <  1:62105/70‑1 (MQ=255)
   gTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCgtt                                                                 >  1:3387577/1‑68 (MQ=255)
    tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCTAATACACGCTc                                                                >  1:853191/1‑70 (MQ=255)
    tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTc                                                                >  1:3319114/1‑70 (MQ=255)
    tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTc                                                                >  1:4016534/1‑70 (MQ=255)
    tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                                                               <  1:716624/71‑1 (MQ=255)
        ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                           >  1:1974117/1‑71 (MQ=255)
        ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                           >  1:1120992/1‑71 (MQ=255)
        ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                           <  1:3821955/71‑1 (MQ=255)
          ctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGcca                                                         >  1:406838/1‑71 (MQ=255)
             ggTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGccacc                                                       <  1:2882579/70‑1 (MQ=255)
               ttATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAAccgccta                                                     >  1:261118/1‑68 (MQ=255)
               ttATGACGTACGGAAATACCTTCGTTTCGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCaccaac                                                    <  1:530284/71‑1 (MQ=255)
                 aTGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCaccaacca                                                  >  1:989482/1‑71 (MQ=255)
                   gACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                                                               >  1:913760/1‑56 (MQ=255)
                   gACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAg                                                            >  1:523210/1‑59 (MQ=255)
                   gACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAg                                                            >  1:242739/1‑59 (MQ=255)
                   gACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAg                                                            >  1:2801985/1‑59 (MQ=255)
                    aCGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAAccacca                                               <  1:2439472/71‑1 (MQ=255)
                     cGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac                                                 >  1:3973233/1‑68 (MQ=255)
                       tACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac                                                 <  1:217576/66‑1 (MQ=255)
                        aCGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACg                                           >  1:108448/1‑71 (MQ=255)
                         cGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGc                                          >  1:1033318/1‑71 (MQ=255)
                          ggAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGc                                          >  1:3762003/1‑70 (MQ=255)
                           gAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCaccaa                                                     >  1:4026208/1‑58 (MQ=255)
                            aaaTACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTg                                       <  1:3909333/71‑1 (MQ=255)
                            aaaTACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTg                                       <  1:2878222/71‑1 (MQ=255)
                            aaaTACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGttt                                        >  1:3601423/1‑70 (MQ=255)
                             aaTACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGa                                      <  1:10027/71‑1 (MQ=255)
                             aaTACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGttt                                        >  1:3770173/1‑69 (MQ=255)
                             aaTACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGttt                                        >  1:3167801/1‑69 (MQ=255)
                               tACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGttt                                        >  1:56374/1‑67 (MQ=255)
                                aCCTTCGTAACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCaccaac                                                    >  1:3004721/1‑54 (MQ=255)
                                   ttCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                                                             <  1:1915596/42‑1 (MQ=255)
                                   ttCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                                                             <  1:1493998/42‑1 (MQ=255)
                                   ttCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                                                             <  1:2806749/42‑1 (MQ=255)
                                   ttCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                                                             <  1:341326/42‑1 (MQ=255)
                                    tCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGcca                                                         >  1:1316429/1‑45 (MQ=255)
                                    tCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGcca                                                         >  1:1037450/1‑45 (MQ=255)
                                    tCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAg                               >  1:1538647/1‑71 (MQ=255)
                                     cGTTACGGAAGTTACGGTTGATTTCGAATACACGc                                                                  >  1:2884728/1‑35 (MQ=255)
                                      gTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGtt                             >  1:1418443/1‑71 (MQ=255)
                                       ttACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAg                                              <  1:4071597/53‑1 (MQ=255)
                                       ttACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTc                            <  1:312358/71‑1 (MQ=255)
                                         aCGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAga                                    >  1:2041964/1‑60 (MQ=255)
                                         aCGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCg                          >  1:53451/1‑71 (MQ=255)
                                         aCGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCg                          >  1:3490400/1‑71 (MQ=255)
                                         aCGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCg                          >  1:255040/1‑71 (MQ=255)
                                          cGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCgg                         >  1:203821/1‑71 (MQ=255)
                                          cGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCg                          <  1:2881085/70‑1 (MQ=255)
                                             aaGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcgc                      >  1:2449018/1‑71 (MQ=255)
                                              aGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcgc                      >  1:1448985/1‑70 (MQ=255)
                                              aGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcgc                      >  1:2900841/1‑70 (MQ=255)
                                               gTTACGGTTGATTTCGAATACACGCTCGAAGCCaccaac                                                    <  1:2432471/39‑1 (MQ=255)
                                               gTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGa                    <  1:3098223/71‑1 (MQ=255)
                                                  aCGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATAc                 >  1:3611396/1‑71 (MQ=255)
                                                   cGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGa                    <  1:451250/67‑1 (MQ=255)
                                                   cGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACg                <  1:803030/71‑1 (MQ=255)
                                                   cGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACg                <  1:807925/71‑1 (MQ=255)
                                                    ggTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCAGGCGCGATACg                <  1:4092885/70‑1 (MQ=255)
                                                    ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc               <  1:3029277/71‑1 (MQ=255)
                                                    ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc               <  1:3872398/71‑1 (MQ=255)
                                                    ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc               <  1:193866/71‑1 (MQ=255)
                                                      ttGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACa                                >  1:1208894/1‑52 (MQ=255)
                                                       tGATTTCGAATACACGCTCGAAGCCACCAAccacca                                               >  1:741788/1‑36 (MQ=255)
                                                         aTTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa          >  1:3631984/1‑71 (MQ=255)
                                                          tttCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGa                    >  1:3823431/1‑60 (MQ=255)
                                                           ttCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTc                            >  1:776922/1‑51 (MQ=255)
                                                           ttCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACa        >  1:3673052/1‑71 (MQ=255)
                                                            tCGAATACACGCTCGACGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACAt       >  1:1742955/1‑71 (MQ=255)
                                                             cGAATACACGCTCGAAGCCACCAACCAACAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACAt       <  1:2267886/70‑1 (MQ=255)
                                                              gAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACAtat     >  1:4106411/1‑69 (MQ=255)
                                                              gAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACAta      >  1:3495070/1‑69 (MQ=255)
                                                              gAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACAta      >  1:3100084/1‑69 (MQ=255)
                                                                 tACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTAc                                 <  1:1927913/40‑1 (MQ=255)
                                                                 tACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTAc                                 <  1:1808327/40‑1 (MQ=255)
                                                                 tACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCga  <  1:3602721/71‑1 (MQ=255)
                                                                   cacGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACa        >  1:2419638/1‑63 (MQ=255)
                                                                   cacGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTc    >  1:744424/1‑67 (MQ=255)
                                                                    |                                                                   
ATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGA  >  minE/1933613‑1933748

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: