Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,581 T→G 16.4% K288Q (AAA→CAA) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,5810TG16.4% 215.3 / 46.2 112K288Q (AAA→CAA) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (45/48);  new base G (8/10);  total (53/59)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 129 total aligned reads displayed.

ACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCG  >  minE/1933511‑1933651
                                                                      |                                                                      
aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                                                        >  1:2426161/1‑71 (MQ=255)
aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                                                        >  1:1174200/1‑71 (MQ=255)
aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                                                        <  1:5454919/71‑1 (MQ=255)
    ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                                                    >  1:4152859/1‑71 (MQ=255)
    ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                                     <  1:4825934/70‑1 (MQ=255)
    ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                                     <  1:3235233/70‑1 (MQ=255)
     cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                                                   <  1:1125447/71‑1 (MQ=255)
     cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                                                   <  1:5629324/71‑1 (MQ=255)
     cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                                                   <  1:4032619/71‑1 (MQ=255)
       tCTTACCGAGAATATCCTGTGACAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                                                        <  1:1831827/64‑1 (MQ=255)
         ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGt                                                                      >  1:363106/1‑64 (MQ=255)
         ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGc                                                               >  1:2824540/1‑71 (MQ=255)
           aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                                                       <  1:2674842/61‑1 (MQ=255)
               agaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                                         >  1:251775/1‑71 (MQ=255)
                 aaTATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                                                       <  1:3865103/71‑1 (MQ=255)
                    atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                                                     >  1:4320873/1‑70 (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:4458211/70‑1 (MQ=255)
                       cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                                                  <  1:2918867/70‑1 (MQ=255)
                       cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                                                  <  1:1164551/70‑1 (MQ=255)
                        tgtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTagag                                                >  1:489102/1‑71 (MQ=255)
                        tgtgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                                                   >  1:3499498/1‑52 (MQ=255)
                         gtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGt                                               >  1:2505897/1‑71 (MQ=255)
                         gtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGt                                               >  1:4180999/1‑71 (MQ=255)
                         gtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGt                                               >  1:2978081/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:5620026/38‑1 (MQ=255)
                                  tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                                                       <  1:4564623/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: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)
                                      gAACAGTGACTCTGTCAGTTCAATCAAATCGTAGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg                                   <  1:207200/70‑1 (MQ=255)
                                      gAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                                   >  1:4119688/1‑70 (MQ=255)
                                       aaCAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGaa                                 >  1:4333181/1‑71 (MQ=255)
                                         cAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAAct                               <  1:4462152/71‑1 (MQ=255)
                                         cAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAct                               >  1:573568/1‑71 (MQ=255)
                                          aGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActc                              >  1:3809577/1‑71 (MQ=255)
                                           gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                             >  1:4744690/1‑71 (MQ=255)
                                           gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                             >  1:2719605/1‑71 (MQ=255)
                                           gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                             >  1:3805927/1‑71 (MQ=255)
                                           gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                             >  1:5368644/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:2781024/1‑49 (MQ=255)
                                              aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                                                >  1:2863214/1‑49 (MQ=255)
                                              aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                                                  <  1:1526961/47‑1 (MQ=255)
                                                 cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATTGTGAAct                               <  1:4907958/63‑1 (MQ=255)
                                                 cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAct                               <  1:2997913/63‑1 (MQ=255)
                                                 cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTa                       >  1:524734/1‑71 (MQ=255)
                                                 cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTa                       >  1:724880/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)
                                                      aGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACg                  >  1:1337695/1‑71 (MQ=255)
                                                         tCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                                   >  1:168352/1‑35 (MQ=255)
                                                          cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                             <  1:3499983/56‑1 (MQ=255)
                                                          cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                             <  1:290926/56‑1 (MQ=255)
                                                          cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTggg                          >  1:4682924/1‑59 (MQ=255)
                                                          cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACg              >  1:5580556/1‑71 (MQ=255)
                                                          cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCAATCATGGTGAACTCTggg                          >  1:2521287/1‑59 (MQ=255)
                                                           gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTAc               <  1:4813837/69‑1 (MQ=255)
                                                           gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTAc               <  1:2618209/69‑1 (MQ=255)
                                                            aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg                                     <  1:3918940/46‑1 (MQ=255)
                                                            aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg                                     <  1:3696514/46‑1 (MQ=255)
                                                             taaGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa            <  1:2749264/68‑1 (MQ=255)
                                                               aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa          >  1:4235832/1‑70 (MQ=255)
                                                                 aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                          <  1:5137049/36‑1 (MQ=255)
                                                                 aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGaa                                 <  1:1727524/45‑1 (MQ=255)
                                                                   cTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACg              >  1:3107546/1‑62 (MQ=255)
                                                                    tttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACg                  <  1:5237807/57‑1 (MQ=255)
                                                                     ttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTa                       >  1:110065/1‑51 (MQ=255)
                                                                      tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                                   <  1:1784391/38‑1 (MQ=255)
                                                                      tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                                   <  1:70909/38‑1 (MQ=255)
                                                                      tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTa                >  1:4402170/1‑57 (MQ=255)
                                                                      tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTc   <  1:281655/70‑1 (MQ=255)
                                                                      tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCg  >  1:4408375/1‑71 (MQ=255)
                                                                      |                                                                      
ACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCG  >  minE/1933511‑1933651

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: