Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,294,089 C→A 43.6% intergenic (+116/‑159) 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,0890CA43.6% 74.1 / 70.5 71intergenic (+116/‑159)gldA/yijFglycerol dehydrogenase, NAD/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/21);  new base A (18/13);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

TGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATT  >  minE/2294022‑2294159
                                                                   |                                                                      
tGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGa                                                                     >  1:1324775/1‑71 (MQ=255)
tGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGa                                                                     >  1:2831294/1‑71 (MQ=255)
 gATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGa                                                                     >  1:4901081/1‑70 (MQ=255)
  aTGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGAt                                                                    <  1:4455300/70‑1 (MQ=255)
  aTGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATa                                                                   <  1:973445/71‑1 (MQ=255)
  aTGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGAt                                                                    >  1:1621580/1‑70 (MQ=255)
   tGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATaa                                                                  <  1:4767790/71‑1 (MQ=255)
     cgcTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAg                                                                 <  1:2057525/70‑1 (MQ=255)
       cTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCg                                                              <  1:1017138/71‑1 (MQ=255)
       cTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCg                                                              <  1:581835/71‑1 (MQ=255)
          cGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGtt                                                            <  1:4114046/70‑1 (MQ=255)
            cTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCAc                                                         >  1:3613057/1‑71 (MQ=255)
              tATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCg                                                              <  1:54883/64‑1 (MQ=255)
              tATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCg                                                              <  1:2644355/64‑1 (MQ=255)
               aTCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTcc                                                      >  1:1376865/1‑71 (MQ=255)
                tCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGAt                                                                    <  1:2993334/56‑1 (MQ=17)
                tCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGtt                                                            <  1:1285377/64‑1 (MQ=255)
                 cAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCg                                                     >  1:4720792/1‑70 (MQ=255)
                  aGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATaa                                                                  >  1:1977081/1‑56 (MQ=17)
                  aGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATaa                                                                  >  1:592788/1‑56 (MQ=17)
                  aGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCa                                                   <  1:3097563/71‑1 (MQ=255)
                    gCCTACGCCATCTCTGCAATATATTGAATTTGCTTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATc                                                 <  1:5565589/71‑1 (MQ=14)
                    gCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATc                                                 <  1:350415/71‑1 (MQ=14)
                      cTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATaa                                                                  >  1:4650485/1‑52 (MQ=14)
                      cTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCg                                               >  1:3390699/1‑71 (MQ=11)
                       tACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGtt                                                            <  1:2706362/57‑1 (MQ=255)
                       tACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGtt                                                            <  1:5361768/57‑1 (MQ=255)
                        aCGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGAt                                                                    <  1:776566/48‑1 (MQ=255)
                         cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATaa                                                                  >  1:827485/1‑49 (MQ=14)
                         cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATaa                                                                  >  1:4195964/1‑49 (MQ=14)
                         cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATaa                                                                  >  1:2887121/1‑49 (MQ=14)
                         cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTc                                                           >  1:4768019/1‑56 (MQ=14)
                         cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCgg                                              <  1:1258126/69‑1 (MQ=11)
                         cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCgg                                              <  1:2120277/69‑1 (MQ=11)
                         cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCa                                            >  1:3472742/1‑71 (MQ=11)
                         cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCa                                            >  1:2113870/1‑71 (MQ=11)
                         cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCa                                            >  1:4381388/1‑71 (MQ=11)
                         cGCAATCTCTGCAATATAATGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCTTCCgg                                              <  1:4109549/69‑1 (MQ=255)
                              tctcTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGAt                                                                    <  1:1518826/42‑1 (MQ=255)
                              tctcTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATa                                                                   <  1:4591468/43‑1 (MQ=255)
                               ctctGTAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCg                                                              >  1:749650/1‑47 (MQ=255)
                               ctctGCAATATATTGAATTTGCGTGCTTTTGTAAGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCtgt                                      >  1:2420119/1‑71 (MQ=255)
                                tctGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCa                                            >  1:4102370/1‑64 (MQ=40)
                                tctGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATaa                                                                  >  1:4010011/1‑42 (MQ=255)
                                 ctGCAATATATTGAATTTGCGTGCTTTTTTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGc                                    >  1:5678167/1‑71 (MQ=255)
                                 ctGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGtt                                                            <  1:5568961/47‑1 (MQ=38)
                                 ctGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGtt                                                            <  1:3646079/47‑1 (MQ=38)
                                  tGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGc                                    >  1:4100520/1‑70 (MQ=255)
                                    cAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCa                                                          <  1:2137637/46‑1 (MQ=38)
                                    cAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCa                                                          <  1:5045231/46‑1 (MQ=38)
                                    cAATATATTGAATTTGCGTGCTTTTGTAGGCAGGAt                                                                    <  1:50652/36‑1 (MQ=255)
                                     aaTATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCActct                                        >  1:3623469/1‑63 (MQ=255)
                                      atatatTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGt                                >  1:1044988/1‑70 (MQ=255)
                                      atatatTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGt                                <  1:1771631/70‑1 (MQ=255)
                                      atatatTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTa                               <  1:3342011/71‑1 (MQ=255)
                                      atatatTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTa                               <  1:2666049/71‑1 (MQ=255)
                                         tatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCa                                                          >  1:4161596/1‑41 (MQ=255)
                                         tatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCa                                                          >  1:1820872/1‑41 (MQ=255)
                                          atTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCa                                                          >  1:704235/1‑40 (MQ=40)
                                           ttGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCt                           >  1:540453/1‑70 (MQ=255)
                                            tGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTc                                                           >  1:5021569/1‑37 (MQ=255)
                                             gAATTTGCGTGCTTTTGTAGGCCTGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGa                         <  1:770221/70‑1 (MQ=255)
                                             gAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGa                         <  1:2242911/70‑1 (MQ=255)
                                               aTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCt                           >  1:3459319/1‑66 (MQ=255)
                                               aTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCt                           >  1:3457347/1‑66 (MQ=255)
                                                tttGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTc                                                       <  1:150415/37‑1 (MQ=37)
                                                 ttGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAAt                             >  1:2507610/1‑62 (MQ=255)
                                                 ttGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAAt                             >  1:815692/1‑62 (MQ=255)
                                                  tGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATcc                                                >  1:3887064/1‑42 (MQ=38)
                                                            tGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACgttgtt          >  1:543986/1‑70 (MQ=255)
                                                             gTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCActc                                         <  1:1144228/38‑1 (MQ=40)
                                                                  ccGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTc         >  1:1958157/1‑65 (MQ=255)
                                                                   cGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAAtt  >  1:5140115/1‑71 (MQ=255)
                                                                   |                                                                      
TGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATT  >  minE/2294022‑2294159

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