Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,492 T→C 50.9% intergenic (‑117/+100) 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,4920TC50.9% 23.0 / 77.4 59intergenic (‑117/+100)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base T (17/12);  new base C (14/16);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCT  >  minE/1614422‑1614559
                                                                      |                                                                   
aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                                                     >  1:3856953/1‑71 (MQ=38)
aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                                                     >  1:1797315/1‑71 (MQ=38)
aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                                                     <  1:5201586/71‑1 (MQ=38)
aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                                                     >  1:2797065/1‑71 (MQ=38)
   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)
       tATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                                                     >  1:5187757/1‑64 (MQ=39)
       tATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCAc                                                                >  1:5126278/1‑69 (MQ=255)
        aTCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACcgc                                                             >  1:2750132/1‑71 (MQ=14)
        aTCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCAcc                                                               <  1:502711/69‑1 (MQ=255)
             cccTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCAc                                                        <  1:4959017/71‑1 (MQ=255)
                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)
                        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: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)
                                                         ttgCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTa                      <  1:2801232/59‑1 (MQ=255)
                                                           gCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAAtt           <  1:4596660/70‑2 (MQ=21)
                                                           gCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATcc          >  1:2122357/1‑71 (MQ=255)
                                                             aTCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCg        <  1:2698853/71‑1 (MQ=255)
                                                              tCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCg        >  1:4818242/1‑70 (MQ=18)
                                                              tCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCg        >  1:1423884/1‑70 (MQ=18)
                                                              tCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTAt                              >  1:762628/1‑48 (MQ=255)
                                                              tCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGt       <  1:1751831/71‑1 (MQ=21)
                                                                cGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTa      <  1:2153166/70‑1 (MQ=18)
                                                                cGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTa      <  1:1830344/70‑1 (MQ=255)
                                                                cGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTa      <  1:1247214/70‑1 (MQ=255)
                                                                cGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAg     <  1:3673076/71‑1 (MQ=18)
                                                                  gCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGc   >  1:3780197/1‑71 (MQ=255)
                                                                   cAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCt  >  1:4676485/1‑71 (MQ=255)
                                                                      |                                                                   
AACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCT  >  minE/1614422‑1614559

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: