Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,614,492 T→C 32.6% 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,4920TC32.6% 73.6 / 40.6 49intergenic (‑117/+100)eutH/eutGpredicted inner membrane protein/predicted alcohol dehydrogenase in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base T (17/16);  new base C (11/5);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

AACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGG  >  minE/1614422‑1614561
                                                                      |                                                                     
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)
        aTCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCAc                                                                  >  1:3483256/1‑68 (MQ=11)
             cccTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCa                                                           <  1:2301621/70‑1 (MQ=14)
             cccTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCAc                                                          >  1:1945120/1‑71 (MQ=14)
               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:1231018/1‑48 (MQ=255)
                           tCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCa                                                                   >  1:249682/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:2490576/71‑1 (MQ=255)
                                                      gCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGc                 <  1:2718191/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)
                                                                      |                                                                     
AACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGGATAAGACGCGGCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGCAATCCCGTAGGCTGG  >  minE/1614422‑1614561

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: