Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,472 G→A 14.1% intergenic (‑97/+120) eutH ← / ← eutG predicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,614,4720GA14.1% 130.0 / 13.5 57intergenic (‑97/+120)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base A (4/4);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTA  >  minE/1614405‑1614539
                                                                   |                                                                   
ctGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGTCCGGATAAAACGCGGCaa                                                                  <  1:3373962/70‑1 (MQ=11)
   gCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCg                                                               <  1:489257/71‑1 (MQ=39)
   gCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCg                                                               <  1:3663009/71‑1 (MQ=39)
    cGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGt                                                              <  1:2718424/71‑1 (MQ=39)
       aTGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGttg                                                            >  1:3091395/1‑68 (MQ=14)
        tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCa                                                          >  1:737851/1‑71 (MQ=14)
         tcggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCa                                                          <  1:1101341/69‑1 (MQ=255)
         gcggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCAt                                                         >  1:3368692/1‑71 (MQ=14)
         aaggcgTTCACGCCGCATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCAt                                                         >  1:1572756/3‑71 (MQ=2)
             cgTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCgg                                                     >  1:2173851/1‑71 (MQ=14)
              gTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCgg                                                     >  1:2058987/1‑70 (MQ=14)
              gTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGc                                                    <  1:1384635/71‑1 (MQ=39)
              gTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTCGATAAGACGCGGCAAGCGTCGCATCCGGc                                                    <  1:3520332/71‑1 (MQ=255)
               tAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCa                                                   <  1:807156/71‑1 (MQ=38)
                 aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCaa                                                  <  1:2802027/70‑1 (MQ=37)
                 aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                                 <  1:3107984/71‑1 (MQ=38)
                  aCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAtt                                                >  1:2060079/1‑71 (MQ=39)
                      cTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCAc                                            >  1:2111889/1‑71 (MQ=14)
                      cTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCAc                                            >  1:3467131/1‑71 (MQ=255)
                       ttATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGc                                                           >  1:3689995/1‑55 (MQ=11)
                         aTCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCAc                                            >  1:3483256/1‑68 (MQ=11)
                          tCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCa                                                          >  1:2040593/1‑53 (MQ=2)
                            cGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCg                                                               <  1:2638082/46‑1 (MQ=11)
                            cGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGTCAAGCGTCGc                                                           >  1:3795167/1‑50 (MQ=255)
                              cccTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCa                                     <  1:2301621/70‑1 (MQ=14)
                              cccTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCAc                                    >  1:1945120/1‑71 (MQ=14)
                              cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGc                                                    <  1:2341939/55‑1 (MQ=11)
                              cccTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGc                                                                <  1:3373022/43‑1 (MQ=9)
                                cTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTg                                  >  1:2991800/1‑71 (MQ=255)
                                cTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGAAAGCGTCGCATCCGGCAATTGCACCGCGCCACTg                                  <  1:1665207/71‑1 (MQ=255)
                                 tACATGTGCAATTCCGTAACCTGGATAAAACGCGACa                                                                   <  1:600062/37‑1 (MQ=9)
                                  aCATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGc                                >  1:222925/1‑71 (MQ=255)
                                   cATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTgg                                 <  1:3574051/69‑1 (MQ=14)
                                         cAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcg                         <  1:3180184/71‑1 (MQ=255)
                                           aTCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTg                                  <  1:2503289/60‑1 (MQ=255)
                                            tCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCa                                             >  1:1231018/1‑48 (MQ=255)
                                            tCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCa                                             >  1:249682/1‑48 (MQ=255)
                                            tCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcggcg                      >  1:3905088/1‑71 (MQ=255)
                                              ccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTg                    <  1:3388477/71‑1 (MQ=255)
                                               cGTAGGGTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGa                   <  1:2687519/71‑1 (MQ=255)
                                                gTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCAc                                    >  1:2109157/1‑53 (MQ=255)
                                                gTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGaa                  <  1:2877095/71‑1 (MQ=255)
                                                gTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGAATTGCACCGCGCCACTGGCGGATGCGGCGTGAAc                 >  1:3734222/1‑71 (MQ=255)
                                                cctggCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGa                   >  1:1652619/4‑70 (MQ=14)
                                                 ctggCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATTCGGCGTGAAc                 <  1:3133703/69‑1 (MQ=14)
                                                   ggCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGc               >  1:3817798/1‑71 (MQ=14)
                                                     ccggATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCt             >  1:3689378/3‑71 (MQ=17)
                                                     cTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCt             <  1:1488824/71‑1 (MQ=255)
                                                      tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcggcg                      <  1:1222409/61‑1 (MQ=255)
                                                      tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcggcg                      <  1:3131435/61‑1 (MQ=255)
                                                      tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGc               <  1:2296956/68‑1 (MQ=255)
                                                      tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCt             >  1:2963760/1‑70 (MQ=255)
                                                        gATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGGGCGTGAACGcc              >  1:2853831/1‑66 (MQ=17)
                                                              aCGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGc                                <  1:3119725/43‑1 (MQ=40)
                                                              aCGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGc                                <  1:904734/43‑1 (MQ=40)
                                                              aCGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGccc    >  1:2719307/1‑71 (MQ=255)
                                                                gcgGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTa  >  1:661283/1‑71 (MQ=255)
                                                                   |                                                                   
TTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTA  >  minE/1614405‑1614539

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: