Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,453 G→C 19.2% intergenic (‑78/+139) 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,4530GC19.2% 112.0 / 20.6 52intergenic (‑78/+139)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base G (17/25);  new base C (7/3);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAAC  >  minE/1614383‑1614524
                                                                      |                                                                       
aCTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCg                                                                         >  1:3132824/1‑71 (MQ=255)
    ccTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAgg                                                                     <  1:1286241/71‑1 (MQ=255)
    ccTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAgg                                                                     <  1:2418148/71‑1 (MQ=255)
    ccTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAgg                                                                     <  1:3537392/71‑1 (MQ=255)
        ttAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTgg                                                                 <  1:2373627/71‑1 (MQ=255)
         tAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGa                                                                >  1:1885955/1‑71 (MQ=255)
           gttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATa                                                              <  1:3582369/71‑1 (MQ=255)
            ttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGAt                                                               >  1:2110406/1‑69 (MQ=255)
            ttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATaa                                                             <  1:2063967/71‑1 (MQ=255)
            ttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATaa                                                             <  1:3809581/71‑1 (MQ=255)
                 atttattGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGAcgc                                                        >  1:1700644/1‑71 (MQ=255)
                    tattGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCg                                                                         >  1:882715/1‑51 (MQ=255)
                      ctGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGTCCGGATAAAACGCGGCaa                                                   <  1:3373962/70‑1 (MQ=11)
                         gCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCg                                                <  1:3663009/71‑1 (MQ=39)
                         gCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCg                                                <  1:489257/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)
                                    gTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAgg                                                                     >  1:2210462/1‑39 (MQ=40)
                                    gTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAgg                                                                     >  1:2925568/1‑39 (MQ=40)
                                    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)
                                                                      |                                                                       
ACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAAC  >  minE/1614383‑1614524

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: