Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,720 T→C 11.7% E241E (GAA→GAG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,7200TC11.7% 205.6 / 17.8 77E241E (GAA→GAGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base C (7/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGG  >  minE/1933651‑1933790
                                                                     |                                                                      
gTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGtt                                                                       >  1:1418443/1‑71 (MQ=255)
 ttACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTc                                                                      <  1:312358/71‑1 (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:3872398/71‑1 (MQ=255)
              ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                                         <  1:193866/71‑1 (MQ=255)
              ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                                         <  1:3029277/71‑1 (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)
                           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: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)
                                               caaccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTAtgatgg                        >  1:3355044/3‑71 (MQ=255)
                                               caaccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTAtgatgg                        >  1:2820590/3‑71 (MQ=255)
                                                caccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtggg                       >  1:1790564/1‑71 (MQ=255)
                                                 accaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                                                    <  1:188923/41‑1 (MQ=255)
                                                 accaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCa                                      >  1:723509/1‑55 (MQ=255)
                                                  ccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtg                         <  1:3342585/67‑1 (MQ=255)
                                                  ccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtgggtg                     <  1:3637099/71‑1 (MQ=255)
                                                  ccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtgggt                      >  1:2704845/1‑70 (MQ=255)
                                                    aGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtgggtg                     >  1:1325526/1‑69 (MQ=255)
                                                     gACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCa                                      >  1:118621/1‑51 (MQ=255)
                                                     gACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                  <  1:146971/71‑1 (MQ=255)
                                                      aCGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGtaataa                 <  1:708954/71‑1 (MQ=255)
                                                      aCGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaa                 <  1:1254616/71‑1 (MQ=255)
                                                       cGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCa                                      >  1:727859/1‑49 (MQ=255)
                                                        gCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggt                          >  1:3728706/1‑60 (MQ=255)
                                                          ttGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtg                         >  1:2975726/1‑59 (MQ=255)
                                                            gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTa                              >  1:2845528/1‑52 (MQ=255)
                                                            gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtg                            >  1:4072739/1‑54 (MQ=255)
                                                                  cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt                                >  1:2624653/1‑44 (MQ=255)
                                                                  cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCa      <  1:1590013/70‑1 (MQ=255)
                                                                  cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg     >  1:2046856/1‑71 (MQ=255)
                                                                   aGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaaa                <  1:3268151/59‑1 (MQ=255)
                                                                   aGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGc    <  1:2405991/71‑1 (MQ=255)
                                                                    gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaaa                <  1:3847449/58‑1 (MQ=255)
                                                                     ttCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCgg  >  1:662680/1‑71 (MQ=255)
                                                                     |                                                                      
GTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGG  >  minE/1933651‑1933790

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: