Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,459 T→C 31.5% intergenic (‑84/+133) 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,4590TC31.5% 86.1 / 42.5 54intergenic (‑84/+133)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base T (14/23);  new base C (12/5);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 3.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

TTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCT  >  minE/1614391‑1614528
                                                                    |                                                                     
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)
                            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)
                                                                   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)
                                                                    |                                                                     
TTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCT  >  minE/1614391‑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: