Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,522 A→G 14.3% intergenic (‑147/+70) 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,5220AG14.3% 91.9 / 9.3 42intergenic (‑147/+70)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base A (19/17);  new base G (4/2);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

CGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGATGCGGCAAGCGTCGCATCCGG  >  minE/1614452‑1614588
                                                                      |                                                                  
cGTAGGGTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGa                                                                    <  1:2687519/71‑1 (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)
      ccggATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCt                                                              >  1:3689378/3‑71 (MQ=17)
      cTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCt                                                              <  1:1488824/71‑1 (MQ=255)
       tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGc                                                                <  1:2296956/68‑1 (MQ=255)
       tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCt                                                              >  1:2963760/1‑70 (MQ=255)
         gATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGGGCGTGAACGcc                                                               >  1:2853831/1‑66 (MQ=17)
               aCGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGccc                                                     >  1:2719307/1‑71 (MQ=255)
                 gcgGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTa                                                   >  1:661283/1‑71 (MQ=255)
                   gacAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACa                                                 >  1:854993/3‑71 (MQ=18)
                   gacAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACa                                                 <  1:297707/69‑1 (MQ=18)
                      aaGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACAtgt                                              >  1:165621/1‑71 (MQ=255)
                        gCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACAtgtg                                             >  1:1576231/1‑70 (MQ=255)
                        gCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGc                                            <  1:2718191/71‑1 (MQ=255)
                        gCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGc                                            <  1:2490576/71‑1 (MQ=255)
                          gTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTa                                                            >  1:3626005/1‑53 (MQ=255)
                            cGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACAt                                                >  1:3702230/1‑63 (MQ=255)
                             gCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAAtt                                        <  1:15493/70‑2 (MQ=21)
                               aTCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCg                                     <  1:1921948/71‑1 (MQ=255)
                               aTCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCg                                     >  1:3449703/1‑71 (MQ=255)
                                tCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCg                                     >  1:2230178/1‑70 (MQ=18)
                                  cGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTaa                                  >  1:3194727/1‑70 (MQ=18)
                                    gCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGc                                                       >  1:3748466/1‑48 (MQ=255)
                                      aaTTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTg                              >  1:443419/1‑71 (MQ=255)
                                       aTTGCTCCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTgg                             <  1:3856386/71‑1 (MQ=255)
                                          gCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGAt                           <  1:2509347/70‑1 (MQ=255)
                                           cACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTa                                                            >  1:1603914/1‑36 (MQ=255)
                                           cACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGAt                           <  1:2057435/69‑1 (MQ=17)
                                            aCCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAg                        <  1:2203437/71‑1 (MQ=17)
                                              cgcgCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATTTGTGTTCCCGTAGGTCGGATAAGac                      <  1:203905/71‑2 (MQ=14)
                                              cgcgCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGAt                      <  1:1091102/71‑1 (MQ=255)
                                                cgcCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATaa                         <  1:2515505/66‑1 (MQ=255)
                                                     cTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACaa               <  1:899545/71‑1 (MQ=11)
                                                     cTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGt                                    <  1:2281915/50‑1 (MQ=40)
                                                     cTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGATGCGGc                 >  1:1167919/1‑69 (MQ=255)
                                                        gCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGATGCGGCAAGc             <  1:2220183/70‑1 (MQ=40)
                                                          ggATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGATGCGGCAAGCGTc          >  1:2382209/1‑71 (MQ=39)
                                                            aTGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGATg                     >  1:3285937/1‑58 (MQ=40)
                                                                ggcgTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATaa                         <  1:1779004/50‑1 (MQ=11)
                                                                  agTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGATGCGGCAAGCGTCGCATCCgg  >  1:2324430/2‑71 (MQ=39)
                                                                   gTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGATGCGGCAAGCGTCGCATCCgg  >  1:1865315/1‑70 (MQ=40)
                                                                      |                                                                  
CGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGATGCGGCAAGCGTCGCATCCGG  >  minE/1614452‑1614588

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: