Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,480 C→T 34.2% intergenic (‑105/+112) 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,4800CT34.2% 94.0 / 63.6 76intergenic (‑105/+112)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base C (28/22);  new base T (13/13);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCA  >  minE/1614413‑1614547
                                                                   |                                                                   
tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCa                                                                  <  1:3901698/71‑1 (MQ=14)
tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCa                                                                  <  1:3641343/71‑1 (MQ=39)
tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCa                                                                  <  1:4221789/71‑1 (MQ=39)
tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCa                                                                  <  1:2416348/71‑1 (MQ=39)
  cggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATc                                                                >  1:2534512/1‑71 (MQ=14)
   ggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATcc                                                               <  1:2929521/71‑1 (MQ=39)
   ggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATcc                                                               <  1:3777624/71‑1 (MQ=39)
     cgTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCgg                                                             >  1:2523114/1‑71 (MQ=39)
      gTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGc                                                            <  1:4820010/71‑1 (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:3353167/1‑70 (MQ=11)
         aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                                                         >  1:3491071/1‑70 (MQ=11)
         aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                                                         >  1:2938918/1‑70 (MQ=11)
         aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                                                         >  1:658721/1‑70 (MQ=11)
         aaCGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaac                                                         >  1:671011/1‑70 (MQ=11)
         aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                                         >  1:2797065/1‑71 (MQ=38)
         aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                                         >  1:1797315/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)
                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)
                      cccTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCaa                                                          >  1:364833/1‑57 (MQ=2)
                      cccTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCAc                                            <  1:4959017/71‑1 (MQ=255)
                         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)
                           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)
                                            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:1463460/1‑71 (MQ=255)
                                                gATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTAt                  >  1:4773168/1‑71 (MQ=255)
                                                 aTAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATc                 >  1:2344458/1‑71 (MQ=18)
                                                  tAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATcc                >  1:4978755/1‑71 (MQ=255)
                                                   aaaaCGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCg               >  1:1016883/4‑71 (MQ=18)
                                                   aaaaCGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCg               >  1:5645550/4‑71 (MQ=18)
                                                     gACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGc              >  1:4877157/1‑70 (MQ=255)
                                                       cgcgACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGccc            >  1:2658491/1‑70 (MQ=17)
                                                         cgGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTAc         <  1:5262566/71‑1 (MQ=255)
                                                          gacAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACa        >  1:2468543/3‑71 (MQ=18)
                                                             aaGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACAtgt     <  1:3885686/71‑1 (MQ=255)
                                                             aaGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACAtgt     >  1:606767/1‑71 (MQ=21)
                                                               gCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGc   <  1:3671484/71‑1 (MQ=21)
                                                                 gTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCa  <  1:939301/70‑1 (MQ=255)
                                                                   |                                                                   
TGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCA  >  minE/1614413‑1614547

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: