Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,696 A→T 20.3% V249V (GTT→GTA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,6960AT20.3% 129.0 / 29.7 59V249V (GTT→GTAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base T (10/2);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 9.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGT  >  minE/1933628‑1933766
                                                                    |                                                                      
ttATGACGTACGGAAATACCTTCGTTTCGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCaccaac                                                                      <  1:530284/71‑1 (MQ=255)
  aTGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCaccaacca                                                                    >  1:989482/1‑71 (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:2878222/71‑1 (MQ=255)
             aaaTACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTg                                                         <  1:3909333/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)
                     tCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAg                                                 >  1:1538647/1‑71 (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:255040/1‑71 (MQ=255)
                          aCGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCg                                            >  1:3490400/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:2900841/1‑70 (MQ=255)
                               aGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcgc                                        >  1:1448985/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)
                                                      cGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGa                                      <  1:3037342/49‑1 (MQ=255)
                                                      cGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATc                >  1:1476339/1‑71 (MQ=255)
                                                         tCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAg             >  1:1432559/1‑71 (MQ=255)
                                                           gAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTc                      <  1:3342719/60‑1 (MQ=255)
                                                             aGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGcgcg         >  1:4067940/1‑71 (MQ=255)
                                                              gCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCTCGATACGCAGGTACATGTCGAGATCCAGCGCGt        <  1:47521/71‑1 (MQ=255)
                                                              gCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt        >  1:2937929/1‑71 (MQ=255)
                                                              gCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt        <  1:3471807/71‑1 (MQ=255)
                                                               ccaccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                            <  1:1699890/50‑1 (MQ=255)
                                                               ccaccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                            <  1:3563919/50‑1 (MQ=255)
                                                               ccaccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGtt       <  1:2867274/71‑1 (MQ=255)
                                                                caccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                            <  1:2427662/49‑1 (MQ=255)
                                                                caccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                            <  1:2265559/49‑1 (MQ=255)
                                                                 accaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt     <  1:3777217/71‑1 (MQ=255)
                                                                 accaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt     <  1:940860/71‑1 (MQ=255)
                                                                  ccaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAgg                              <  1:2162682/45‑1 (MQ=255)
                                                                  ccaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAgcgc          <  1:1302375/65‑1 (MQ=255)
                                                                  ccaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtg    <  1:1778782/71‑1 (MQ=255)
                                                                    aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                            >  1:3974055/1‑45 (MQ=255)
                                                                    aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggt  >  1:2263193/1‑71 (MQ=255)
                                                                    aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggt  <  1:4006010/71‑1 (MQ=255)
                                                                    |                                                                      
TTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGT  >  minE/1933628‑1933766

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: