Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,453 G→C 31.3% 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,4530GC31.3% 114.2 / 66.3 83intergenic (‑78/+139)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base C (12/14);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

CTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTG  >  minE/1614384‑1614521
                                                                     |                                                                    
cTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCg                                                                      <  1:1808626/70‑1 (MQ=255)
   ccTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAgg                                                                  <  1:865819/71‑1 (MQ=255)
       ttAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTg                                                               <  1:3206098/70‑1 (MQ=255)
              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:4221789/71‑1 (MQ=39)
                             tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCa                                        <  1:3641343/71‑1 (MQ=39)
                             tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCa                                        <  1:2416348/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:658721/1‑70 (MQ=11)
                                      aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                               >  1:671011/1‑70 (MQ=11)
                                      aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                               >  1:3353167/1‑70 (MQ=11)
                                      aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                               >  1:3491071/1‑70 (MQ=11)
                                      aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                               >  1:1797315/1‑71 (MQ=38)
                                      aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                               >  1:2797065/1‑71 (MQ=38)
                                      aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                               <  1:5201586/71‑1 (MQ=38)
                                      aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                               >  1:3856953/1‑71 (MQ=38)
                                       aCGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCatct                              >  1:1726505/1‑68 (MQ=11)
                                         gCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTg                             >  1:38098/1‑70 (MQ=11)
                                         gCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTg                             >  1:236485/1‑70 (MQ=11)
                                         gCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTg                             >  1:5672081/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:4944970/60‑1 (MQ=11)
                                              aTCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGc                                  <  1:2264733/60‑1 (MQ=11)
                                              aTCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGc                                  <  1:793690/60‑1 (MQ=11)
                                              aTCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCAcc                         <  1:502711/69‑1 (MQ=255)
                                                   cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACaa                                                <  1:4477590/41‑1 (MQ=9)
                                                   cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACaa                                                <  1:4348493/41‑1 (MQ=9)
                                                   cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACaa                                                <  1:1971297/41‑1 (MQ=9)
                                                   cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACaa                                                <  1:401474/41‑1 (MQ=9)
                                                   cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACaa                                                <  1:3316219/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:5323426/1‑71 (MQ=255)
                                                      tACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTgg               >  1:3925711/1‑71 (MQ=255)
                                                        cATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCg             <  1:3817910/71‑1 (MQ=14)
                                                        cATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCg             <  1:807848/71‑1 (MQ=14)
                                                        cATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCg             <  1:1671830/71‑1 (MQ=14)
                                                        cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGc                                              >  1:587354/1‑38 (MQ=40)
                                                        cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCg                                             >  1:5282379/1‑39 (MQ=40)
                                                        cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCg                                             >  1:3073076/1‑39 (MQ=40)
                                                        cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCg                                             >  1:3850333/1‑39 (MQ=40)
                                                        cATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCg             <  1:4556734/71‑1 (MQ=255)
                                                         aTGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCg             >  1:4747650/1‑70 (MQ=255)
                                                           gtgCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGcgt                                >  1:2478874/1‑47 (MQ=11)
                                                           gtgCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGcgt                                >  1:2978124/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:26239/1‑58 (MQ=255)
                                                                   ccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTgg               >  1:1764484/1‑58 (MQ=255)
                                                                   ccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGt   >  1:1779074/1‑70 (MQ=255)
                                                                    cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTg  >  1:4811753/1‑70 (MQ=11)
                                                                     |                                                                    
CTCCCTGTTAGTTGTTATTTATTGGCGGATGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTG  >  minE/1614384‑1614521

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: