Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,575 G→A 6.7% L290L (CTG→TTG)  lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,5750GA 6.7% 288.9 / 7.6 89L290L (CTG→TTG) lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (50/33);  new base A (4/2);  total (54/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

TAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAAT  >  minE/1933505‑1933644
                                                                      |                                                                     
tAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAg                                                                       <  1:1616318/71‑1 (MQ=255)
    tCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCt                                                                   <  1:1739048/71‑1 (MQ=255)
         tCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGa                                                                      >  1:2497037/1‑63 (MQ=255)
         tCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                                                >  1:959493/1‑69 (MQ=255)
          ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                                             >  1:3578698/1‑71 (MQ=255)
           cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATc                                                                    <  1:2226686/63‑1 (MQ=255)
            gTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCtt                                                                  >  1:2053609/1‑64 (MQ=255)
              cTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTg                                                         >  1:2775208/1‑71 (MQ=255)
               ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGc                                                        >  1:2941054/1‑71 (MQ=255)
               ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGAGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGc                                                        >  1:1184279/1‑71 (MQ=255)
                 aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGc                                                        <  1:436360/69‑1 (MQ=255)
                 aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCg                                                       >  1:1048979/1‑70 (MQ=255)
                 aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCg                                                       >  1:3987714/1‑70 (MQ=255)
                  ccGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                              >  1:998915/1‑62 (MQ=255)
                  ccGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                              >  1:3244948/1‑62 (MQ=255)
                  ccGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGt                                                      >  1:3504792/1‑70 (MQ=255)
                   cGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                              <  1:3081496/61‑1 (MQ=255)
                    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:1540856/1‑70 (MQ=255)
                          atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                                              >  1:710732/1‑70 (MQ=255)
                             cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGAt                                                                     >  1:3978467/1‑44 (MQ=255)
                                tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                              <  1:2126579/48‑1 (MQ=255)
                                tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                       >  1:2024923/1‑71 (MQ=255)
                                tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                       >  1:3752500/1‑71 (MQ=255)
                                tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                       >  1:3873091/1‑71 (MQ=255)
                                 gCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                                  >  1:311818/1‑59 (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)
                                       gTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                                       <  1:4030477/64‑1 (MQ=255)
                                        tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                                              <  1:3224893/40‑1 (MQ=255)
                                        tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                                  <  1:2990111/52‑1 (MQ=255)
                                         gcGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCcatcattg                              >  1:805300/2‑69 (MQ=255)
                                         aCGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg                              >  1:777398/1‑71 (MQ=255)
                                          cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATc                                                           >  1:1096470/1‑41 (MQ=255)
                                          cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATc                                                           >  1:1953407/1‑41 (MQ=255)
                                          cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                                          >  1:533369/1‑42 (MQ=255)
                                          cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGcc                                                 >  1:2928884/1‑51 (MQ=255)
                                          cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                             >  1:582045/1‑71 (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)
                                               cAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcatc                                  <  1:272688/61‑1 (MQ=255)
                                               cAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAct                        >  1:3825783/1‑71 (MQ=255)
                                                 gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                                         >  1:2408506/1‑52 (MQ=255)
                                                 gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                      <  1:1261435/71‑1 (MQ=255)
                                                   gATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg                    >  1:3818034/1‑71 (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:852328/71‑1 (MQ=255)
                                                         gTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATg              <  1:2209943/71‑1 (MQ=255)
                                                            aGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg                    <  1:2262911/62‑1 (MQ=255)
                                                             gCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACg           <  1:3572921/70‑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:3216173/71‑1 (MQ=255)
                                                                  aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa     <  1:968724/71‑1 (MQ=255)
                                                                  aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa     <  1:3899779/71‑1 (MQ=255)
                                                                  aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa     <  1:1979392/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)
                                                                      |                                                                     
TAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAAT  >  minE/1933505‑1933644

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: