Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,294,050 A→C 26.0% intergenic (+77/‑198) gldA → / → yijF glycerol dehydrogenase, NAD/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,294,0500AC26.0% 103.4 / 53.2 73intergenic (+77/‑198)gldA/yijFglycerol dehydrogenase, NAD/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/24);  new base C (14/5);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGG  >  minE/2293983‑2294115
                                                                   |                                                                 
aaCTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAAtc                                                                >  1:995316/1‑71 (MQ=255)
aaCTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAAtc                                                                >  1:3502041/1‑71 (MQ=255)
 aCTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAAtct                                                               <  1:3113736/71‑1 (MQ=255)
      cGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCa                                                          <  1:379824/71‑1 (MQ=255)
       ggCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCaa                                                         >  1:3294144/1‑71 (MQ=255)
        gCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCaa                                                         >  1:1973101/1‑70 (MQ=255)
            gTCCAGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAatata                                                     >  1:1330705/1‑70 (MQ=255)
              ccGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATAtt                                                   >  1:3253906/1‑70 (MQ=255)
                 ggAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAt                                               >  1:1401993/1‑71 (MQ=255)
                 ggAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAt                                               >  1:2611737/1‑71 (MQ=255)
                 ggAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAt                                               >  1:2128710/1‑71 (MQ=255)
                   aGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAttt                                             <  1:3367968/71‑1 (MQ=255)
                    gTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTg                                            <  1:1289934/71‑1 (MQ=255)
                     tttGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAttt                                             <  1:1598888/69‑1 (MQ=255)
                      ttGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGc                                           <  1:1254907/70‑1 (MQ=255)
                      ttGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCg                                          >  1:3073908/1‑71 (MQ=255)
                       tGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAAtctc                                                              <  1:2943384/50‑1 (MQ=255)
                       tGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGc                                           >  1:47669/1‑69 (MQ=255)
                         accGCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGaa                                                >  1:1897338/3‑62 (MQ=14)
                         accGCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGc                                       <  1:3647828/69‑1 (MQ=14)
                         aaCGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGc                                       >  1:657255/1‑71 (MQ=255)
                          ccGCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCt                                      >  1:1373522/2‑71 (MQ=255)
                          aCGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGc                                           >  1:3782285/1‑66 (MQ=255)
                          aCGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCt                                      <  1:352729/71‑1 (MQ=255)
                          aCGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCt                                      <  1:2571491/71‑1 (MQ=255)
                           cGCAAAAATGCCTGATGCGCTACGCTTATCAGGCGTACGCCATCTCTGCAATATATTGAAtt                                              >  1:1733740/1‑62 (MQ=14)
                            gCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATAtt                                                   <  1:2066470/56‑1 (MQ=255)
                              aaaaATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                  >  1:62027/1‑71 (MQ=255)
                              aaaaATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                  >  1:3496625/1‑71 (MQ=255)
                               aaaaTGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGt                                 <  1:799777/71‑1 (MQ=255)
                               aaaTTGCCTTATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                  >  1:2628688/1‑70 (MQ=255)
                               aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                  >  1:3815840/1‑70 (MQ=255)
                                aaatGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                  >  1:960635/4‑69 (MQ=255)
                                aaTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCttt                                    >  1:1960770/1‑67 (MQ=255)
                                aaTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                  >  1:311150/1‑69 (MQ=255)
                                aaTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                  >  1:3048736/1‑69 (MQ=255)
                                    gCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTa                                >  1:2573273/1‑67 (MQ=255)
                                    gCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTa                                >  1:3302748/1‑67 (MQ=255)
                                     ccTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCtt                                     <  1:2686658/61‑1 (MQ=255)
                                     ccTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGt                                 <  1:465331/65‑1 (MQ=255)
                                         aTGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAAt                                               <  1:1818995/47‑1 (MQ=255)
                                         aTGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGAt                        >  1:1089080/1‑70 (MQ=255)
                                         aTGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGAt                        >  1:2431973/1‑70 (MQ=255)
                                             gcTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTg                                  <  1:1548148/56‑1 (MQ=255)
                                               tACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAg                     >  1:1461516/1‑67 (MQ=255)
                                               tACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGATTTTGTAGGCAGGATAAg                     >  1:2979552/1‑67 (MQ=14)
                                                 cGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTc               <  1:1533679/71‑1 (MQ=255)
                                                 cGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGtt                >  1:1554278/1‑70 (MQ=255)
                                                 cGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGtt                >  1:293116/1‑70 (MQ=255)
                                                 cGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGt                 <  1:3173641/69‑1 (MQ=255)
                                                    ttATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCg                  <  1:321679/65‑1 (MQ=255)
                                                    ttATCAGGCCTAAGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACt            >  1:2043816/1‑71 (MQ=255)
                                                      aTCAGGCCTACGCAATCTCTGCAATATATTGAAtt                                              <  1:1490344/35‑1 (MQ=255)
                                                      aTCAGGCCTACGCAATCTCTGCAATATATTGAAtt                                              <  1:2318912/35‑1 (MQ=255)
                                                      aTCAGGCCTACGCAATCTCTGCAATATATTGAATTTg                                            <  1:1669030/37‑1 (MQ=255)
                                                      aTCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGc                                       >  1:529753/1‑42 (MQ=255)
                                                       tCAGGCCTACGCAATCTCTGCAATATATTGAATTTGc                                           >  1:3107459/1‑37 (MQ=255)
                                                       tCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTcc          <  1:1103155/70‑1 (MQ=255)
                                                       tCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCg         <  1:1062469/71‑1 (MQ=255)
                                                         aGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGc                   >  1:414829/1‑59 (MQ=17)
                                                         aGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCa       >  1:2111750/1‑71 (MQ=14)
                                                         aGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCa       >  1:2116405/1‑71 (MQ=14)
                                                         aGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCAc             >  1:1103871/1‑65 (MQ=255)
                                                         aGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCa       >  1:972255/1‑71 (MQ=255)
                                                          ggCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCAt      <  1:3878604/71‑1 (MQ=255)
                                                           gCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATc     <  1:2697775/71‑1 (MQ=255)
                                                            ccTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATcc    >  1:560551/1‑71 (MQ=11)
                                                              tACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTCGTa                                <  1:960924/41‑1 (MQ=255)
                                                               aCGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCgg  >  1:3306699/1‑70 (MQ=255)
                                                                cGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGt                                 >  1:3496042/1‑38 (MQ=255)
                                                                cGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCgg  <  1:3844082/69‑1 (MQ=255)
                                                                  cAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATa                       <  1:3845639/46‑1 (MQ=255)
                                                                  cAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATcc    >  1:2099575/1‑65 (MQ=255)
                                                                   aaTCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATa                       <  1:769919/45‑1 (MQ=255)
                                                                   aaTCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATa                       <  1:1773593/45‑1 (MQ=255)
                                                                   |                                                                 
AACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGG  >  minE/2293983‑2294115

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: