Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,480 C→T 21.4% 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,4800CT21.4% 115.4 / 25.6 56intergenic (‑105/+112)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base C (21/23);  new base T (10/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 4.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGC  >  minE/1614413‑1614546
                                                                   |                                                                  
tGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCa                                                                 >  1:737851/1‑71 (MQ=14)
 tcggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCa                                                                 <  1:1101341/69‑1 (MQ=255)
 gcggcgTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCAt                                                                >  1:3368692/1‑71 (MQ=14)
 aaggcgTTCACGCCGCATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCAt                                                                >  1:1572756/3‑71 (MQ=2)
     cgTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCgg                                                            >  1:2173851/1‑71 (MQ=14)
      gTAAACGCCTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCgg                                                            >  1:2058987/1‑70 (MQ=14)
      gTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGc                                                           <  1:1384635/71‑1 (MQ=39)
      gTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTCGATAAGACGCGGCAAGCGTCGCATCCGGc                                                           <  1:3520332/71‑1 (MQ=255)
       tAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCa                                                          <  1:807156/71‑1 (MQ=38)
         aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCaa                                                         <  1:2802027/70‑1 (MQ=37)
         aaCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAt                                                        <  1:3107984/71‑1 (MQ=38)
          aCGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAAtt                                                       >  1:2060079/1‑71 (MQ=39)
              cTTATCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCAc                                                   >  1:2111889/1‑71 (MQ=14)
              cTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCAc                                                   >  1:3467131/1‑71 (MQ=255)
               ttATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGc                                                                  >  1:3689995/1‑55 (MQ=11)
                 aTCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCAc                                                   >  1:3483256/1‑68 (MQ=11)
                  tCCGCCCTACATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCa                                                                 >  1:2040593/1‑53 (MQ=2)
                    cGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGTCAAGCGTCGc                                                                  >  1:3795167/1‑50 (MQ=255)
                      cccTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCa                                            <  1:2301621/70‑1 (MQ=14)
                      cccTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCAc                                           >  1:1945120/1‑71 (MQ=14)
                      cccTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGc                                                           <  1:2341939/55‑1 (MQ=11)
                        cTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTg                                         >  1:2991800/1‑71 (MQ=255)
                        cTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGAAAGCGTCGCATCCGGCAATTGCACCGCGCCACTg                                         <  1:1665207/71‑1 (MQ=255)
                          aCATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGc                                       >  1:222925/1‑71 (MQ=255)
                           cATGTGCAATCCTCCTGGCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTgg                                        <  1:3574051/69‑1 (MQ=14)
                                 cAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcg                                <  1:3180184/71‑1 (MQ=255)
                                   aTCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTg                                         <  1:2503289/60‑1 (MQ=255)
                                    tCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCa                                                    >  1:249682/1‑48 (MQ=255)
                                    tCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCa                                                    >  1:1231018/1‑48 (MQ=255)
                                    tCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcggcg                             >  1:3905088/1‑71 (MQ=255)
                                      ccGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTg                           <  1:3388477/71‑1 (MQ=255)
                                       cGTAGGGTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGa                          <  1:2687519/71‑1 (MQ=255)
                                        gTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCAc                                           >  1:2109157/1‑53 (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)
                                              tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcggcg                             <  1:3131435/61‑1 (MQ=255)
                                              tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATgcggcg                             <  1:1222409/61‑1 (MQ=255)
                                              tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGc                      <  1:2296956/68‑1 (MQ=255)
                                              tGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCt                    >  1:2963760/1‑70 (MQ=255)
                                                gATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGGGCGTGAACGcc                     >  1:2853831/1‑66 (MQ=17)
                                                      aCGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGc                                       <  1:3119725/43‑1 (MQ=40)
                                                      aCGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGc                                       <  1:904734/43‑1 (MQ=40)
                                                      aCGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGccc           >  1:2719307/1‑71 (MQ=255)
                                                        gcgGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTa         >  1:661283/1‑71 (MQ=255)
                                                          gacAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACa       <  1:297707/69‑1 (MQ=18)
                                                          gacAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACa       >  1:854993/3‑71 (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)
                                                                   |                                                                  
TGCGGCGTAAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGC  >  minE/1614413‑1614546

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: