Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,466 G→A 35.2% intergenic (‑91/+126) 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,4660GA35.2% 76.3 / 47.4 54intergenic (‑91/+126)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base G (15/20);  new base A (12/7);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

ATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCT  >  minE/1614400‑1614528
                                                                  |                                                              
atttattGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGAcgc                                                            >  1:1700644/1‑71 (MQ=255)
     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:249682/1‑48 (MQ=255)
                                                 tCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCa                                  >  1:1231018/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:3131435/61‑1 (MQ=255)
                                                           tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcggcg           <  1:1222409/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)
                                                                  |                                                              
ATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCT  >  minE/1614400‑1614528

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: