Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,699 C→A 19.9% V248V (GTG→GTT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,6990CA19.9% 171.9 / 19.4 60V248V (GTG→GTTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (23/25);  new base A (10/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 4.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

ATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGG  >  minE/1933630‑1933769
                                                                     |                                                                      
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)
           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)
                   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: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:2900841/1‑70 (MQ=255)
                             aGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcgc                                           >  1:1448985/1‑70 (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:807925/71‑1 (MQ=255)
                                  cGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACg                                     <  1:803030/71‑1 (MQ=255)
                                   ggTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCAGGCGCGATACg                                     <  1:4092885/70‑1 (MQ=255)
                                   ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                    <  1:193866/71‑1 (MQ=255)
                                   ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                    <  1:3029277/71‑1 (MQ=255)
                                   ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                    <  1:3872398/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:3100084/1‑69 (MQ=255)
                                             gAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACAta                           >  1:3495070/1‑69 (MQ=255)
                                                tACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTAc                                                      <  1:1808327/40‑1 (MQ=255)
                                                tACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTAc                                                      <  1:1927913/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:3471807/71‑1 (MQ=255)
                                                            gCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt           >  1:2937929/1‑71 (MQ=255)
                                                             ccaccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                               <  1:3563919/50‑1 (MQ=255)
                                                             ccaccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                               <  1:1699890/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:4006010/71‑1 (MQ=255)
                                                                  aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggt     >  1:2263193/1‑71 (MQ=255)
                                                                   accaccAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtg    >  1:860261/1‑71 (MQ=255)
                                                                   accaccAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggt     >  1:4068966/1‑70 (MQ=255)
                                                                     caccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtggg  >  1:1790564/1‑71 (MQ=255)
                                                                     |                                                                      
ATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGG  >  minE/1933630‑1933769

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: