Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,452 C→T 19.2% intergenic (‑77/+140) 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,4520CT19.2% 111.9 / 19.5 52intergenic (‑77/+140)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base C (19/23);  new base T (7/3);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

ACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGA  >  minE/1614383‑1614522
                                                                     |                                                                      
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:3809581/71‑1 (MQ=255)
            ttgttATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATaa                                                           <  1:2063967/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)
                                                                     |                                                                      
ACTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGA  >  minE/1614383‑1614522

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: