Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,466 G→A 42.9% 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,4660GA42.9% 91.4 / 101.7 91intergenic (‑91/+126)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base G (33/19);  new base A (15/24);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 2.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

TTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGC  >  minE/1614398‑1614535
                                                                    |                                                                     
ttatttatTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGAc                                                                     >  1:3972471/1‑71 (MQ=255)
  atttattGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGAcgc                                                                   >  1:1815634/1‑71 (MQ=255)
  atttattGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGAcg                                                                    <  1:5531776/70‑1 (MQ=255)
    ttattGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCgg                                                                 >  1:2696822/1‑71 (MQ=255)
       ttGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGc                                                                >  1:3541707/1‑69 (MQ=255)
       ctGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCa                                                               <  1:2234559/69‑1 (MQ=11)
        tGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAg                                                             >  1:1456548/1‑71 (MQ=39)
        tGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTAGATAAGACGCGGCAAg                                                             >  1:4947173/1‑71 (MQ=255)
          gCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCg                                                           >  1:3857576/1‑71 (MQ=14)
              aTGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGttg                                                        >  1:1982493/1‑68 (MQ=14)
               tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCa                                                      <  1:3901698/71‑1 (MQ=14)
               tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCa                                                      <  1:2416348/71‑1 (MQ=39)
               tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCa                                                      <  1:3641343/71‑1 (MQ=39)
               tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCa                                                      <  1:4221789/71‑1 (MQ=39)
                 cggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATc                                                    >  1:2534512/1‑71 (MQ=14)
                  ggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATcc                                                   <  1:3777624/71‑1 (MQ=39)
                  ggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATcc                                                   <  1:2929521/71‑1 (MQ=39)
                    cgTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCgg                                                 >  1:2523114/1‑71 (MQ=39)
                     gTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGc                                                <  1:4820010/71‑1 (MQ=14)
                      tAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGttg                                                        <  1:1639171/62‑3 (MQ=14)
                      tAAACGCCTTATCCGCACTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCa                                               <  1:4381879/71‑1 (MQ=11)
                       aaaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCaa                                              >  1:1555259/1‑71 (MQ=38)
                        gaCGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGcgt                                              <  1:5543437/69‑3 (MQ=11)
                        aaCGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTc                                                         >  1:1030981/1‑59 (MQ=11)
                        aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                                             >  1:2938918/1‑70 (MQ=11)
                        aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                                             >  1:3491071/1‑70 (MQ=11)
                        aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                                             >  1:671011/1‑70 (MQ=11)
                        aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                                             >  1:3353167/1‑70 (MQ=11)
                        aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                                             >  1:658721/1‑70 (MQ=11)
                        aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                             >  1:3856953/1‑71 (MQ=38)
                        aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                             >  1:2797065/1‑71 (MQ=38)
                        aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                             >  1:1797315/1‑71 (MQ=38)
                        aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                             <  1:5201586/71‑1 (MQ=38)
                         aCGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCatct                                            >  1:1726505/1‑68 (MQ=11)
                           gCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTg                                           >  1:38098/1‑70 (MQ=11)
                           gCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTg                                           >  1:5672081/1‑70 (MQ=11)
                           gCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTg                                           >  1:236485/1‑70 (MQ=11)
                            ccTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCa                                         <  1:1479959/71‑1 (MQ=14)
                            ccTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCa                                         <  1:3840454/71‑1 (MQ=14)
                            ccTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCa                                         >  1:1914965/1‑71 (MQ=255)
                             cTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGc                                                <  1:1635542/63‑1 (MQ=35)
                              ttATCCGCCCTACATGTGCAATTCCGTAACCTTGATAAAACGCGACAAGCg                                                           <  1:431/51‑1 (MQ=2)
                              ttATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCg                                                           <  1:378486/51‑1 (MQ=11)
                               tATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                             >  1:5187757/1‑64 (MQ=39)
                               tATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCAc                                        >  1:5126278/1‑69 (MQ=255)
                                aTCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACcgc                                     >  1:2750132/1‑71 (MQ=14)
                                aTCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGc                                                <  1:793690/60‑1 (MQ=11)
                                aTCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGc                                                <  1:4944970/60‑1 (MQ=11)
                                aTCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGc                                                <  1:2264733/60‑1 (MQ=11)
                                aTCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCAcc                                       <  1:502711/69‑1 (MQ=255)
                                     cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACaa                                                              <  1:1971297/41‑1 (MQ=9)
                                     cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACaa                                                              <  1:4477590/41‑1 (MQ=9)
                                     cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACaa                                                              <  1:4348493/41‑1 (MQ=9)
                                     cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACaa                                                              <  1:3316219/41‑1 (MQ=9)
                                     cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACaa                                                              <  1:401474/41‑1 (MQ=9)
                                     cccTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaa                                              >  1:364833/1‑57 (MQ=2)
                                     cccTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCAc                                <  1:4959017/71‑1 (MQ=255)
                                      ccTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAg                                                             <  1:2669289/41‑1 (MQ=9)
                                        tACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTc                                                         <  1:5638472/43‑1 (MQ=9)
                                        tACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGc                                                       <  1:502283/45‑1 (MQ=11)
                                        tACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTgg                             >  1:3925711/1‑71 (MQ=255)
                                        tACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTgg                             >  1:5323426/1‑71 (MQ=255)
                                          cATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCg                           <  1:1671830/71‑1 (MQ=14)
                                          cATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCg                           <  1:807848/71‑1 (MQ=14)
                                          cATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCg                           <  1:3817910/71‑1 (MQ=14)
                                          cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGc                                                            >  1:587354/1‑38 (MQ=40)
                                          cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCg                                                           >  1:3850333/1‑39 (MQ=40)
                                          cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCg                                                           >  1:5282379/1‑39 (MQ=40)
                                          cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCg                                                           >  1:3073076/1‑39 (MQ=40)
                                          cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCg                           <  1:4556734/71‑1 (MQ=255)
                                           aTGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCg                           >  1:4747650/1‑70 (MQ=255)
                                             gtgCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGcgt                                              >  1:2978124/1‑47 (MQ=11)
                                             gtgCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGcgt                                              >  1:2478874/1‑47 (MQ=11)
                                                cAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcg                     <  1:5055456/71‑1 (MQ=255)
                                                cAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgc                      <  1:2173726/70‑1 (MQ=255)
                                                 aaTCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcgg                    <  1:3088174/71‑1 (MQ=255)
                                                     ccGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTg                <  1:988845/71‑1 (MQ=11)
                                                     ccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTgg                             >  1:1764484/1‑58 (MQ=255)
                                                     ccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTgg                             >  1:26239/1‑58 (MQ=255)
                                                     ccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGt                 >  1:1779074/1‑70 (MQ=255)
                                                      cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTg                >  1:4811753/1‑70 (MQ=11)
                                                           gtcggATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTg                <  1:4506289/62‑1 (MQ=11)
                                                             cggATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCt         >  1:3433141/2‑70 (MQ=14)
                                                              ggATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTa       >  1:2150067/1‑71 (MQ=255)
                                                              ggATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCtt        <  1:705783/70‑1 (MQ=14)
                                                               gATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTg                              <  1:2918214/47‑1 (MQ=255)
                                                               gATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCg                           <  1:969245/50‑1 (MQ=255)
                                                               gATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTAt      >  1:4773168/1‑71 (MQ=255)
                                                               gATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTAt      >  1:1463460/1‑71 (MQ=255)
                                                                aTAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATc     >  1:2344458/1‑71 (MQ=18)
                                                                 tAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATcc    >  1:4978755/1‑71 (MQ=255)
                                                                    gACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGc  >  1:4877157/1‑70 (MQ=255)
                                                                    |                                                                     
TTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGC  >  minE/1614398‑1614535

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: