Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,594 A→C 9.0% A283A (GCT→GCG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,5940AC 9.0% 264.4 / 26.1 111A283A (GCT→GCGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (56/45);  new base C (7/3);  total (63/48)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 112 total aligned reads displayed.

GAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTA  >  minE/1933525‑1933663
                                                                     |                                                                     
gagaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAg                                                                      >  1:3251239/1‑71 (MQ=255)
gagaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAg                                                                      >  1:787454/1‑71 (MQ=255)
   aaTATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                                                                   >  1:3586161/1‑71 (MQ=255)
     tatCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                                                                 >  1:493825/1‑71 (MQ=255)
      atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                                                     >  1:803082/1‑66 (MQ=255)
      atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                                                                 >  1:710732/1‑70 (MQ=255)
      atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                                                                 >  1:1540856/1‑70 (MQ=255)
            tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                                          >  1:3873091/1‑71 (MQ=255)
            tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                                          >  1:2024923/1‑71 (MQ=255)
            tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                                          >  1:3752500/1‑71 (MQ=255)
             gCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                                               <  1:48026/65‑1 (MQ=255)
             gCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                                               <  1:743968/65‑1 (MQ=255)
             gCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTc                                                         >  1:808657/1‑71 (MQ=255)
              ccAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTc                                                         >  1:3769239/1‑70 (MQ=255)
               cAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTaga                                                             >  1:2596937/1‑65 (MQ=255)
                agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                                      <  1:405187/71‑1 (MQ=255)
                agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                                      <  1:2654301/71‑1 (MQ=255)
                 gagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTc                                                         >  1:2552437/1‑67 (MQ=255)
                   gtgcGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCcatcat                                                   <  1:4062946/68‑1 (MQ=255)
                    tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                                                     <  1:2990111/52‑1 (MQ=255)
                     aCGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg                                                 >  1:777398/1‑71 (MQ=255)
                      cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGcc                                                                    >  1:2928884/1‑51 (MQ=255)
                      cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                                                <  1:3369666/71‑1 (MQ=255)
                        gAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcatc                                                     >  1:4052373/1‑64 (MQ=255)
                          aCAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCcatc                                                     >  1:1070820/1‑62 (MQ=255)
                          aCAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAc                                            >  1:1482309/1‑71 (MQ=255)
                           cAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAct                                           >  1:3825783/1‑71 (MQ=255)
                             gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                                                            >  1:2408506/1‑52 (MQ=255)
                             gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                                         <  1:1261435/71‑1 (MQ=255)
                                aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                                               <  1:1948008/46‑1 (MQ=255)
                                 ttCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGt                                     >  1:2920777/1‑71 (MQ=255)
                                  tCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                                               >  1:3181064/1‑60 (MQ=255)
                                  tCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGt                                     >  1:831731/1‑70 (MQ=255)
                                   cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                                          >  1:2669447/1‑48 (MQ=255)
                                   cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg                                       >  1:909813/1‑67 (MQ=255)
                                     gTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATg                                 <  1:2209943/71‑1 (MQ=255)
                                     gTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATg                                 <  1:852328/71‑1 (MQ=255)
                                        aGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg                                       <  1:2262911/62‑1 (MQ=255)
                                         gCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGt                             <  1:89657/71‑1 (MQ=255)
                                         gCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGt                             <  1:1667434/71‑1 (MQ=255)
                                          cTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                                                            >  1:400692/1‑39 (MQ=255)
                                          cTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                                          >  1:406979/1‑41 (MQ=255)
                                          cTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTAt                                  <  1:1580681/65‑1 (MQ=255)
                                          cTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTa                            >  1:18510/1‑71 (MQ=255)
                                           tCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                                               <  1:1292786/51‑1 (MQ=255)
                                           tCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTAc                           >  1:3990410/1‑71 (MQ=255)
                                             gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTATGACGTACGGa                        <  1:1870443/71‑1 (MQ=255)
                                              aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACTTACGGa                        <  1:533436/71‑1 (MQ=255)
                                              aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                        <  1:3899779/71‑1 (MQ=255)
                                              aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                        <  1:1979392/71‑1 (MQ=255)
                                              aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                        <  1:3216173/71‑1 (MQ=255)
                                              aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                        <  1:968724/71‑1 (MQ=255)
                                              aTCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAAcag                         >  1:3822874/1‑68 (MQ=255)
                                              aTCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGa                        <  1:2550477/71‑1 (MQ=255)
                                               tCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTc                                                         >  1:1651916/1‑37 (MQ=12)
                                                 aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                      >  1:590360/1‑70 (MQ=255)
                                                 aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                        <  1:3124829/68‑1 (MQ=255)
                                                 aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAAt                     <  1:2070404/71‑1 (MQ=255)
                                                 aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAAt                     <  1:310532/71‑1 (MQ=255)
                                                 aGATCTTTGTAAACTGCGTAAGCCATGTAGAGTTCCATCATGGt                                                >  1:2915336/1‑44 (MQ=255)
                                                  gATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                        <  1:2337203/67‑1 (MQ=255)
                                                   aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcatcat                                                   <  1:2705904/39‑1 (MQ=255)
                                                   aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                                      <  1:1104750/36‑1 (MQ=255)
                                                   aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                                      <  1:2762994/36‑1 (MQ=255)
                                                   aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                                               <  1:2020246/43‑1 (MQ=255)
                                                   aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                        <  1:2238646/66‑1 (MQ=255)
                                                   aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAc                   >  1:3417673/1‑71 (MQ=255)
                                                   aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAc                   >  1:3863377/1‑71 (MQ=255)
                                                   aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAc                   >  1:3634999/1‑71 (MQ=255)
                                                   aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAc                   >  1:870379/1‑71 (MQ=255)
                                                   aTCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATAc                   >  1:1151704/1‑71 (MQ=255)
                                                   aTCTTTGTAATATGCGTAAGCCATGTAGAGTTCCATAATGGTGAACTCTGGGTTATg                                 <  1:1832478/57‑1 (MQ=255)
                                                   aTCATTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                                      <  1:1360294/36‑1 (MQ=255)
                                                     cTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGAc                               >  1:617691/1‑57 (MQ=255)
                                                     cTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                      >  1:356757/1‑66 (MQ=255)
                                                      gtggTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCtt                >  1:371941/4‑71 (MQ=255)
                                                      gtggTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCt                 >  1:485024/4‑70 (MQ=255)
                                                       ttGTAATCTGCGTAAGCCATGTAGAGTTCcatcat                                                   >  1:2554648/1‑35 (MQ=255)
                                                       ttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATg                                                  <  1:1406816/36‑1 (MQ=255)
                                                       ttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATg                                                  >  1:3621195/1‑36 (MQ=255)
                                                       tggTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTc               <  1:1457386/69‑1 (MQ=255)
                                                        tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                                         <  1:2881580/44‑1 (MQ=255)
                                                         gTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGt             >  1:2164326/1‑71 (MQ=255)
                                                           aaTCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTa           >  1:1180402/1‑71 (MQ=255)
                                                            aTCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg                                       <  1:1988331/42‑1 (MQ=255)
                                                            aTCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTa           <  1:3319082/70‑1 (MQ=255)
                                                            aTCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTAc          >  1:3393670/1‑71 (MQ=255)
                                                            aTCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTAc          >  1:1745258/1‑71 (MQ=255)
                                                             tCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                      <  1:213005/58‑1 (MQ=255)
                                                              cTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACg         >  1:1379981/1‑70 (MQ=255)
                                                               tGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACgg        >  1:886015/1‑70 (MQ=255)
                                                               tGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACgg        >  1:717129/1‑70 (MQ=255)
                                                               tGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACgg        >  1:2049827/1‑70 (MQ=255)
                                                               tGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACgg        >  1:510700/1‑70 (MQ=255)
                                                                 cGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACg         >  1:363158/1‑67 (MQ=255)
                                                                 cGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAg     >  1:2787459/1‑71 (MQ=255)
                                                                   tAAGCCATGTTGAGTTCCATCATGGTGAACTCTTGGTTATGACGTACGGAAATACCTTCGt             <  1:232386/61‑1 (MQ=255)
                                                                   tAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATg                                 <  1:1330016/41‑1 (MQ=255)
                                                                   tAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGt    >  1:4001823/1‑70 (MQ=255)
                                                                     aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTa  <  1:3801058/70‑1 (MQ=255)
                                                                     |                                                                     
GAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTA  >  minE/1933525‑1933663

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: