Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,454 T→C 18.2% intergenic (‑79/+138) 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,4540TC18.2% 119.5 / 20.5 55intergenic (‑79/+138)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base T (19/26);  new base C (7/3);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAAC  >  minE/1614387‑1614524
                                                                   |                                                                      
ccTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAgg                                                                     <  1:2418148/71‑1 (MQ=255)
ccTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAgg                                                                     <  1:1286241/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:3809581/71‑1 (MQ=255)
        ttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATaa                                                             <  1:2063967/71‑1 (MQ=255)
             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)
                                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: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)
                                                                   |                                                                      
CCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAAC  >  minE/1614387‑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: