Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,555 C→T 10.0% S296S (TCG→TCA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,5550CT10.0% 359.0 / 21.4 130S296S (TCG→TCAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (55/62);  new base T (5/8);  total (60/70)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 131 total aligned reads displayed.

GTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTG  >  minE/1933486‑1933625
                                                                     |                                                                      
gTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCg                                                                       <  1:4424617/71‑1 (MQ=255)
 tCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGa                                                                      >  1:3530621/1‑71 (MQ=255)
     gCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  >  1:616503/1‑71 (MQ=255)
       gTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  >  1:4477600/1‑69 (MQ=255)
        tCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTc                                                               <  1:3177385/71‑1 (MQ=255)
          aCGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  >  1:861474/1‑66 (MQ=255)
          aCGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  >  1:3219102/1‑66 (MQ=255)
          aCGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  >  1:1030118/1‑66 (MQ=255)
            gTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTc                                                                   <  1:5471403/63‑1 (MQ=255)
            gTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTc                                                                   <  1:4816158/63‑1 (MQ=255)
            gTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  <  1:2111644/64‑1 (MQ=255)
            gTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCt                                                           >  1:3284125/1‑71 (MQ=255)
             tCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTc                                                          >  1:5632444/1‑71 (MQ=255)
              cGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCt                                                           >  1:5144250/1‑69 (MQ=255)
               gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                        >  1:2898334/1‑71 (MQ=255)
               gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                        >  1:5327732/1‑71 (MQ=255)
               gCCGTAGGTCACTTCCGTCTTACCGAGAACATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                        <  1:5013603/71‑1 (MQ=255)
                 cGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCa                                                              >  1:2358104/1‑63 (MQ=255)
                 cGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATc                                                      <  1:546588/71‑1 (MQ=255)
                  gTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                       >  1:2143972/1‑69 (MQ=255)
                   tAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAg                                                    <  1:2949391/71‑1 (MQ=255)
                     ggTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGa                                                   <  1:5094003/70‑1 (MQ=255)
                      gTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCa                                                              <  1:529255/58‑1 (MQ=255)
                      gTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCa                                                              <  1:4394475/58‑1 (MQ=255)
                        cACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCg                                                                       >  1:3494102/1‑47 (MQ=255)
                        cACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  >  1:494790/1‑52 (MQ=255)
                        cACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  >  1:283935/1‑52 (MQ=255)
                        cACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  >  1:607336/1‑52 (MQ=255)
                         aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCg                                                                       <  1:1604421/46‑1 (MQ=255)
                         aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                              <  1:5454919/71‑1 (MQ=255)
                         aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                              >  1:2426161/1‑71 (MQ=255)
                         aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                              >  1:1174200/1‑71 (MQ=255)
                           ttCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGc                                                            <  1:3616657/55‑1 (MQ=255)
                           ttCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGc                                                            <  1:3274494/55‑1 (MQ=255)
                             ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATc                                                      >  1:208252/1‑59 (MQ=255)
                             ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                          >  1:3082879/1‑71 (MQ=255)
                             ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                          >  1:4152859/1‑71 (MQ=255)
                             ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                           <  1:4282210/70‑1 (MQ=255)
                             ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                           <  1:4825934/70‑1 (MQ=255)
                              cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                          <  1:2970151/70‑1 (MQ=255)
                              cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                          <  1:525005/70‑1 (MQ=255)
                              cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                         <  1:5629324/71‑1 (MQ=255)
                              cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                         <  1:4032619/71‑1 (MQ=255)
                              cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                         <  1:1125447/71‑1 (MQ=255)
                                tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGt                                                                <  1:3726233/46‑1 (MQ=255)
                                tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGt                                                                <  1:1097762/46‑1 (MQ=255)
                                tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCtt                                               <  1:1070019/63‑1 (MQ=255)
                                tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                       <  1:1550941/71‑1 (MQ=255)
                                  ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGc                                                            <  1:851229/48‑1 (MQ=255)
                                  ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGc                                     >  1:2824540/1‑71 (MQ=255)
                                    aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                             <  1:750842/61‑1 (MQ=255)
                                    aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                             <  1:2674842/61‑1 (MQ=255)
                                        agaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCa                                                              <  1:1806145/40‑1 (MQ=255)
                                        agaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                               >  1:251775/1‑71 (MQ=255)
                                          aaTATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                             <  1:3865103/71‑1 (MQ=255)
                                          aaTATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                             >  1:5144597/1‑71 (MQ=255)
                                             atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                       <  1:1886964/58‑1 (MQ=255)
                                             atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                           >  1:4320873/1‑70 (MQ=255)
                                              tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCg                                                         <  1:4707523/39‑1 (MQ=255)
                                              tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                       <  1:3635940/41‑1 (MQ=255)
                                              tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                           >  1:3359077/1‑69 (MQ=255)
                                              tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                           >  1:2809981/1‑69 (MQ=255)
                                              tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                         <  1:4120994/71‑1 (MQ=255)
                                                cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                        <  1:2918867/70‑1 (MQ=255)
                                                cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                        <  1:4458211/70‑1 (MQ=255)
                                                cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                         >  1:2636665/1‑69 (MQ=255)
                                                 tgtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTagag                      >  1:489102/1‑71 (MQ=255)
                                                  gtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGt                     >  1:4180999/1‑71 (MQ=255)
                                                  gtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGt                     >  1:2505897/1‑71 (MQ=255)
                                                     ccAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTc                   <  1:2002871/70‑1 (MQ=255)
                                                       agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                              >  1:5168495/1‑41 (MQ=255)
                                                       agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                             >  1:4668741/1‑58 (MQ=255)
                                                       agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                <  1:4582410/71‑1 (MQ=255)
                                                       agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                <  1:4214173/71‑1 (MQ=255)
                                                        gagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                <  1:3095157/70‑1 (MQ=255)
                                                         agTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcatc               >  1:1944328/1‑70 (MQ=255)
                                                         agTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCAGCGTa                                  >  1:2250359/1‑51 (MQ=255)
                                                           tACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                              <  1:3446813/37‑1 (MQ=255)
                                                           tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                             <  1:3404465/38‑1 (MQ=255)
                                                           tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                             <  1:4564623/38‑1 (MQ=255)
                                                           tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                             <  1:5620026/38‑1 (MQ=255)
                                                           tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATg            <  1:3973897/71‑1 (MQ=255)
                                                            aCGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                             <  1:3490659/37‑1 (MQ=255)
                                                             cGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGt          <  1:4610910/71‑1 (MQ=255)
                                                             cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt          >  1:3125452/1‑71 (MQ=255)
                                                             cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt          >  1:3218272/1‑71 (MQ=255)
                                                              ggAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg         <  1:2298891/71‑1 (MQ=255)
                                                              ggAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg         <  1:1453187/71‑1 (MQ=255)
                                                              ggAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                         <  1:4345748/55‑1 (MQ=255)
                                                              ggAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg         <  1:1628587/71‑1 (MQ=255)
                                                               gAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg         <  1:3540865/70‑1 (MQ=255)
                                                               gAACAGTGACTCTGTCAGTTCAATCAAATCGTAGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg         <  1:207200/70‑1 (MQ=255)
                                                               gAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg         >  1:4119688/1‑70 (MQ=255)
                                                                aaCAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGaa       >  1:4333181/1‑71 (MQ=255)
                                                                  cAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAct     >  1:573568/1‑71 (MQ=255)
                                                                    gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct   >  1:3805927/1‑71 (MQ=255)
                                                                    gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct   >  1:2719605/1‑71 (MQ=255)
                                                                    gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct   >  1:4744690/1‑71 (MQ=255)
                                                                    gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActc    >  1:4926756/1‑70 (MQ=255)
                                                                     cGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTg  >  1:4254158/1‑71 (MQ=255)
                                                                     |                                                                      
GTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTG  >  minE/1933486‑1933625

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: