Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,570 G→A 13.7% I291I (ATC→ATT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,5700GA13.7% 311.3 / 32.2 124I291I (ATC→ATTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (54/53);  new base A (7/10);  total (61/63)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 127 total aligned reads displayed.

GCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACG  >  minE/1933501‑1933639
                                                                     |                                                                     
gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                                      >  1:2898334/1‑71 (MQ=255)
gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                                      >  1:5327732/1‑71 (MQ=255)
gCCGTAGGTCACTTCCGTCTTACCGAGAACATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                                      <  1:5013603/71‑1 (MQ=255)
  cGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                                     <  1:25966/70‑1 (MQ=255)
  cGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATc                                                                    <  1:3403133/71‑1 (MQ=255)
  cGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATc                                                                    <  1:546588/71‑1 (MQ=255)
   gTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                                     >  1:2143972/1‑69 (MQ=255)
    tAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAg                                                                  <  1:2949391/71‑1 (MQ=255)
      ggTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                                     >  1:4607585/1‑66 (MQ=255)
          aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                                            <  1:5454919/71‑1 (MQ=255)
          aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                                            >  1:2426161/1‑71 (MQ=255)
          aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                                            >  1:1174200/1‑71 (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:4825934/70‑1 (MQ=255)
              ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                         <  1:3235233/70‑1 (MQ=255)
              ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                         <  1:4282210/70‑1 (MQ=255)
               cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                                        <  1:2970151/70‑1 (MQ=255)
               cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                                       <  1:1125447/71‑1 (MQ=255)
               cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                                       <  1:4032619/71‑1 (MQ=255)
               cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                                       <  1:5629324/71‑1 (MQ=255)
                 tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCtt                                                             <  1:1569337/63‑1 (MQ=255)
                 tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCtt                                                             <  1:1070019/63‑1 (MQ=255)
                 tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                                     <  1:4385361/71‑1 (MQ=255)
                 tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                                     <  1:1550941/71‑1 (MQ=255)
                   ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGt                                                          >  1:363106/1‑64 (MQ=255)
                   ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGc                                                   >  1:2824540/1‑71 (MQ=255)
                     aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                                           <  1:750842/61‑1 (MQ=255)
                     aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                                           <  1:2674842/61‑1 (MQ=255)
                         agaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                             >  1:251775/1‑71 (MQ=255)
                           aaTATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                                           >  1:5144597/1‑71 (MQ=255)
                              atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                                     <  1:1886964/58‑1 (MQ=255)
                               tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCg                                                                       <  1:4707523/39‑1 (MQ=255)
                               tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                                         >  1:2809981/1‑69 (MQ=255)
                               tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                       <  1:4120994/71‑1 (MQ=255)
                                 cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                                     <  1:4287362/39‑1 (MQ=255)
                                 cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                                      <  1:1164551/70‑1 (MQ=255)
                                 cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                                      <  1:4458211/70‑1 (MQ=255)
                                 cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                                      <  1:2918867/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)
                                      ccAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTc                                 <  1:2002871/70‑1 (MQ=255)
                                        agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                                            >  1:5168495/1‑41 (MQ=255)
                                        agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                                           >  1:4795297/1‑58 (MQ=255)
                                        agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                              <  1:4582410/71‑1 (MQ=255)
                                        agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                              <  1:4214173/71‑1 (MQ=255)
                                          agTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcatc                             >  1:1944328/1‑70 (MQ=255)
                                          agTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCAGCGTa                                                >  1:2250359/1‑51 (MQ=255)
                                           gTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                                        >  1:2766433/1‑42 (MQ=255)
                                            tACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                                            <  1:3446813/37‑1 (MQ=255)
                                            tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                                           <  1:5620026/38‑1 (MQ=255)
                                            tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                                           <  1:4564623/38‑1 (MQ=255)
                                            tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                                           <  1:3404465/38‑1 (MQ=255)
                                            tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATg                          <  1:3973897/71‑1 (MQ=255)
                                             aCGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                                           <  1:3490659/37‑1 (MQ=255)
                                              cGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGt                        <  1:5607190/71‑1 (MQ=255)
                                              cGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGt                        <  1:4610910/71‑1 (MQ=255)
                                              cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                        >  1:3125452/1‑71 (MQ=255)
                                              cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                        >  1:3174490/1‑71 (MQ=255)
                                              cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                        >  1:3218272/1‑71 (MQ=255)
                                               ggAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg                       <  1:1453187/71‑1 (MQ=255)
                                               ggAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg                       <  1:2298891/71‑1 (MQ=255)
                                               ggAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                       <  1:4345748/55‑1 (MQ=255)
                                                gAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg                       <  1:3540865/70‑1 (MQ=255)
                                                gAACAGTGACTCTGTCAGTTCAATCAAATCGTAGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg                       <  1:207200/70‑1 (MQ=255)
                                                 aaCAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGaa                     >  1:4333181/1‑71 (MQ=255)
                                                  aCAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                        <  1:4098631/67‑1 (MQ=255)
                                                   cAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAct                   >  1:573568/1‑71 (MQ=255)
                                                    aGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActc                  >  1:3809577/1‑71 (MQ=255)
                                                     gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                 >  1:5368644/1‑71 (MQ=255)
                                                     gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                 >  1:2719605/1‑71 (MQ=255)
                                                     gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                 >  1:3805927/1‑71 (MQ=255)
                                                     gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActc                  >  1:4926756/1‑70 (MQ=255)
                                                      cGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTg                >  1:4254158/1‑71 (MQ=255)
                                                        aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                                    >  1:4052027/1‑49 (MQ=255)
                                                        aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                                    >  1:3804555/1‑49 (MQ=255)
                                                        aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                                    >  1:2863214/1‑49 (MQ=255)
                                                        aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                                      <  1:1526961/47‑1 (MQ=255)
                                                        aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg               >  1:1752748/1‑70 (MQ=255)
                                                           cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATTGTGAAct                   <  1:4907958/63‑1 (MQ=255)
                                                           cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAct                   <  1:2997913/63‑1 (MQ=255)
                                                           cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTa           >  1:724880/1‑71 (MQ=255)
                                                           cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTa           >  1:524734/1‑71 (MQ=255)
                                                            tgTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTAt          >  1:4743265/2‑71 (MQ=255)
                                                             gTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATg         <  1:658033/71‑1 (MQ=255)
                                                             gTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATg         <  1:3003014/71‑1 (MQ=255)
                                                              tCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTc                                 >  1:1794025/1‑46 (MQ=255)
                                                              tCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCcatcatt                          <  1:4831082/53‑2 (MQ=255)
                                                              tCAGTTCAATCAAATAGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTa           <  1:3702415/68‑1 (MQ=255)
                                                                aGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACg      >  1:1631772/1‑71 (MQ=255)
                                                                aGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACg      >  1:2144389/1‑71 (MQ=255)
                                                                   tCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                       >  1:168352/1‑35 (MQ=255)
                                                                    cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                 <  1:290926/56‑1 (MQ=255)
                                                                    cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                 <  1:249314/56‑1 (MQ=255)
                                                                    cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                 <  1:3499983/56‑1 (MQ=255)
                                                                    cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTggg              >  1:4682924/1‑59 (MQ=255)
                                                                    cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCAATCATGGTGAACTCTggg              >  1:2521287/1‑59 (MQ=255)
                                                                     gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTAc   <  1:4813837/69‑1 (MQ=255)
                                                                     |                                                                     
GCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACG  >  minE/1933501‑1933639

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