Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,042,434 G→A 7.3% intergenic (+25/+18) yqiI → / ← glgS conserved hypothetical protein/predicted glycogen synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,042,4340GA 7.3% 168.0 / 2.1 54intergenic (+25/+18)yqiI/glgSconserved hypothetical protein/predicted glycogen synthesis protein
Reads supporting (aligned to +/- strand):  ref base G (21/29);  new base A (2/2);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAG  >  minE/2042365‑2042503
                                                                     |                                                                     
ggACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGc                                                                      <  1:1978340/71‑1 (MQ=255)
ggACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGc                                                                      <  1:805324/71‑1 (MQ=255)
ggACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGc                                                                      <  1:3187351/71‑1 (MQ=255)
 gACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGc                                                                      >  1:1945862/1‑70 (MQ=255)
  aCCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGcc                                                                     <  1:1469672/70‑1 (MQ=255)
  aCCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGcc                                                                     <  1:190154/70‑1 (MQ=255)
   ccGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCAt                                                                   <  1:3511188/71‑1 (MQ=255)
    cGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATc                                                                  <  1:2921990/71‑1 (MQ=255)
       ttAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAgg                                                               >  1:1070606/1‑71 (MQ=255)
       ttAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAgg                                                               >  1:701341/1‑71 (MQ=255)
           tAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCg                                                           <  1:3273886/71‑1 (MQ=255)
           tAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCg                                                           <  1:2901664/71‑1 (MQ=255)
             aaGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAg                                                                >  1:2359344/1‑64 (MQ=255)
             aaGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAg                                                                >  1:352614/1‑64 (MQ=255)
             aaGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGt                                                          >  1:1361337/1‑70 (MQ=255)
               gACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTtat                                                       <  1:1190105/71‑1 (MQ=255)
                aCGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTtata                                                      >  1:844686/1‑71 (MQ=255)
                        cTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGc                                                                      <  1:1295658/47‑1 (MQ=255)
                        cTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGc                                               <  1:520781/70‑1 (MQ=255)
                         tGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGctc                                             >  1:3181041/1‑71 (MQ=255)
                               tCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTc                                       <  1:747529/71‑1 (MQ=255)
                                 aCCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGt                                                          <  1:2732033/50‑1 (MQ=255)
                                 aCCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTcc                                      <  1:3476240/70‑1 (MQ=255)
                                 aCCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTc                                       >  1:3643146/1‑69 (MQ=255)
                                     ataACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAg                                    <  1:1649378/68‑1 (MQ=255)
                                     ataACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAg                                    <  1:1428593/68‑1 (MQ=255)
                                     ataACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAg                                    <  1:3389900/68‑1 (MQ=255)
                                     ataACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctc                                 >  1:1642660/1‑71 (MQ=255)
                                     ataACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctc                                 >  1:1579670/1‑71 (MQ=255)
                                        aCTAGTATCTAATACAAACACTAAAACGGGCCATCAGGc                                                              >  1:2520414/1‑39 (MQ=255)
                                        aCTAGTATCTAATACAAACACTAAAACGGGCCATCAGGc                                                              >  1:2586466/1‑39 (MQ=255)
                                        aCTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctctct                              <  1:2650749/71‑1 (MQ=255)
                                          tAGTATCTAATACAAACACTAAAACGGGCCATCAGGc                                                              >  1:681119/1‑37 (MQ=255)
                                                 tAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAAc                                         <  1:2311405/51‑1 (MQ=255)
                                                 tAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAAc                                         <  1:2118050/51‑1 (MQ=255)
                                                   aTACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcatca                   <  1:683524/71‑1 (MQ=255)
                                                   aTACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGc                       <  1:2408770/67‑1 (MQ=255)
                                                    tACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcatc                    <  1:1293595/69‑1 (MQ=255)
                                                    tACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcatc                    <  1:3333605/69‑1 (MQ=255)
                                                      cAAACACTAAAACGGGCCATCAGGCCCGTTATATCa                                                   >  1:602686/1‑36 (MQ=255)
                                                      cAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                <  1:2035602/71‑1 (MQ=255)
                                                      cAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                <  1:2933160/71‑1 (MQ=255)
                                                      cAAACACTAAAACGGACCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                >  1:3157438/1‑71 (MQ=255)
                                                      cAAACACTAAAACGGACCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                >  1:3608714/1‑71 (MQ=255)
                                                       aaaCACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTg               >  1:2294168/1‑71 (MQ=255)
                                                        aaCACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCGTCATCTGt              <  1:1829182/71‑1 (MQ=255)
                                                           acTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGa           >  1:823842/1‑71 (MQ=255)
                                                             tAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGAc          >  1:3121205/1‑70 (MQ=255)
                                                             tAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGAc          >  1:2653032/1‑70 (MQ=255)
                                                               aaaCGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCt                         >  1:3278011/1‑53 (MQ=255)
                                                               aaaCGGACCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACaa        <  1:2958320/70‑1 (MQ=255)
                                                               aaaCGGACCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACaa        <  1:2459091/70‑1 (MQ=255)
                                                                 aCGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAc     >  1:1475204/1‑71 (MQ=255)
                                                                  cGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAcc    <  1:1743299/71‑1 (MQ=255)
                                                                    ggCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAg  <  1:2524075/71‑1 (MQ=255)
                                                                     gCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTg               <  1:2518098/57‑1 (MQ=255)
                                                                     gCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTg               <  1:3149551/57‑1 (MQ=255)
                                                                     |                                                                     
GGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAG  >  minE/2042365‑2042503

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: