Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,579 T→G 16.9% K288N (AAA→AAC) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,5790TG16.9% 211.0 / 47.7 112K288N (AAA→AAC) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (45/48);  new base G (8/11);  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 = 1.00e+00

Only 100 of 123 total aligned reads displayed.

CACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAA  >  minE/1933510‑1933643
                                                                     |                                                                
cACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCtt                                                                 >  1:2054902/1‑71 (MQ=255)
 aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                                                >  1:1174200/1‑71 (MQ=255)
     ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                                            >  1:4152859/1‑71 (MQ=255)
     ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                                            >  1:3082879/1‑71 (MQ=255)
     ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                                            >  1:2175117/1‑71 (MQ=255)
     ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                             <  1:4825934/70‑1 (MQ=255)
     ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                             <  1:4282210/70‑1 (MQ=255)
      cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                                            <  1:525005/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)
        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:2674842/61‑1 (MQ=255)
                agaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                                 >  1:251775/1‑71 (MQ=255)
                  aaTATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                                               <  1:3865103/71‑1 (MQ=255)
                     atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                                         <  1:1886964/58‑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:1164551/70‑1 (MQ=255)
                         tgtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTagag                                        >  1:489102/1‑71 (MQ=255)
                         tgtgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                             <  1:4200930/50‑1 (MQ=255)
                          gtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGt                                       >  1:4180999/1‑71 (MQ=255)
                          gtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGt                                       >  1:2505897/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)
                               agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                                               >  1:4795297/1‑58 (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:3404465/38‑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:4610910/71‑1 (MQ=255)
                                     cGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGt                            <  1:5607190/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)
                                       gAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg                           <  1:3540865/70‑1 (MQ=255)
                                        aaCAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGaa                         >  1:4333181/1‑71 (MQ=255)
                                        aaCAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGaa                         <  1:655701/71‑1 (MQ=255)
                                         aCAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                            <  1:4098631/67‑1 (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: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: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:249314/56‑1 (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:2618209/69‑1 (MQ=255)
                                                            gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTAc       <  1:4813837/69‑1 (MQ=255)
                                                             aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg                             <  1:3696514/46‑1 (MQ=255)
                                                             aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg                             <  1:3918940/46‑1 (MQ=255)
                                                              taaGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa    <  1:2749264/68‑1 (MQ=255)
                                                              tCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGa    >  1:5617105/1‑70 (MQ=255)
                                                                  aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                  <  1:5137049/36‑1 (MQ=255)
                                                                  aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                  <  1:4623723/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)
                                                                     tttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACg          <  1:1190534/57‑1 (MQ=255)
                                                                     |                                                                
CACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAA  >  minE/1933510‑1933643

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: