Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,292,617 A→G 51.7% intergenic (+168/‑253) fsaB → / → gldA fructose‑6‑phosphate aldolase 2/glycerol dehydrogenase, NAD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,292,6170AG51.7% 4.6 / 76.0 58intergenic (+168/‑253)fsaB/gldAfructose‑6‑phosphate aldolase 2/glycerol dehydrogenase, NAD
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base G (14/16);  total (27/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

TATTAGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGTGCTGGC  >  minE/2292549‑2292687
                                                                    |                                                                      
tattaGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGAc                                                                      >  1:1806757/1‑71 (MQ=255)
tattaGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGAc                                                                      >  1:1717600/1‑71 (MQ=255)
  ttaGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACgg                                                                    <  1:2191904/71‑1 (MQ=255)
   taGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATACCCAGGGCGGAAACgg                                                                    >  1:3649662/1‑70 (MQ=255)
     gCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCAt                                                                 >  1:2736979/1‑71 (MQ=255)
     gCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCAt                                                                 >  1:1264819/1‑71 (MQ=255)
     gCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCAt                                                                 >  1:3872384/1‑71 (MQ=255)
       cGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTc                                                               >  1:1099856/1‑71 (MQ=255)
       cGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTc                                                               >  1:2359933/1‑71 (MQ=255)
        gCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCg                                                              <  1:2436150/71‑1 (MQ=255)
        gCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCg                                                              >  1:1162525/1‑71 (MQ=255)
           cTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTAc                                                           >  1:1889125/1‑71 (MQ=255)
                aGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGtt                                                       >  1:1522448/1‑70 (MQ=255)
                     cAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAg                                                 <  1:2490910/71‑1 (MQ=255)
                     cAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAg                                                 >  1:2686493/1‑71 (MQ=255)
                       aaaTACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGc                                                <  1:3519674/70‑1 (MQ=255)
                            cGTTGCTCCGTATGTTAACCGCGTAGATGACCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCaaa                                          >  1:804589/1‑71 (MQ=255)
                             gTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCaaa                                          >  1:1622517/1‑70 (MQ=255)
                                gCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCt                                      >  1:3718908/1‑71 (MQ=255)
                                 cTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCAtt                                                                <  1:3338192/44‑1 (MQ=255)
                                 cTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCAtt                                                                <  1:2547203/44‑1 (MQ=255)
                                 cTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCAtt                                                                <  1:3877418/44‑1 (MQ=255)
                                 cTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCt                                      <  1:1339555/70‑1 (MQ=255)
                                 cTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCa                                            <  1:416351/64‑1 (MQ=255)
                                  tCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAgg                                                   <  1:131046/56‑1 (MQ=255)
                                  tCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAg                                                 <  1:1059568/58‑1 (MQ=255)
                                    cGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGt                                    >  1:681528/1‑69 (MQ=255)
                                    cGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGc                                             >  1:3431525/1‑60 (MQ=255)
                                    cGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGt                                    <  1:3690559/69‑1 (MQ=255)
                                     gTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGtt                                   <  1:1984820/69‑1 (MQ=255)
                                       aTGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTc                                                      <  1:1725823/48‑1 (MQ=255)
                                       aTGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGa                                >  1:3084988/1‑70 (MQ=255)
                                       aTGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGa                                >  1:922830/1‑70 (MQ=255)
                                        tGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGaa                               <  1:94970/70‑1 (MQ=255)
                                        tGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGaa                               <  1:3300024/70‑1 (MQ=255)
                                        tGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGaa                               <  1:1368391/70‑1 (MQ=255)
                                          ttAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAAt                             <  1:3882769/70‑1 (MQ=255)
                                           tAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCg                                                              >  1:1487228/1‑36 (MQ=255)
                                            aaCCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGc                           <  1:1163760/70‑1 (MQ=255)
                                             aCCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCAc                         <  1:2731668/71‑1 (MQ=255)
                                             aCCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCAc                         >  1:1969961/1‑71 (MQ=255)
                                             aCCGCCTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACCCTGTTAGAAATGCAc                         <  1:501305/71‑1 (MQ=255)
                                                 cgTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGTAGCTGCAAACGCTGTTAGAAATGCACgcgc                     <  1:178876/71‑1 (MQ=255)
                                                  gTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCAcgcg                      >  1:2832897/1‑69 (MQ=255)
                                                   tAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCt                                               >  1:2700980/1‑43 (MQ=255)
                                                   tAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCa                                                     <  1:962730/37‑1 (MQ=255)
                                                   tAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCa                                                     <  1:3085428/37‑1 (MQ=255)
                                                       tGCCCAGGGCGGAAACGTCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGaaa               <  1:499470/71‑1 (MQ=255)
                                                        gCCCAGGGCGGAGACGGCATTCGTATGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAg              >  1:383819/1‑71 (MQ=255)
                                                         cccAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAg              <  1:2537589/70‑1 (MQ=255)
                                                         cccAGGGCGGAGACGACATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAg              <  1:1091303/70‑1 (MQ=255)
                                                          ccAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCa            >  1:2009602/1‑71 (MQ=255)
                                                            aGGGCGGAGACGACATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATg          >  1:1373837/1‑71 (MQ=255)
                                                            aGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATg          >  1:3626611/1‑71 (MQ=255)
                                                              ggCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGt        <  1:705746/71‑1 (MQ=255)
                                                              ggCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGt        <  1:3157331/71‑1 (MQ=255)
                                                               gCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGTg       <  1:3375811/71‑1 (MQ=255)
                                                                cGGAAAGGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATTGTGc      <  1:1687632/71‑1 (MQ=255)
                                                                   aaaCGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGTGCTgg   >  1:2364351/1‑71 (MQ=255)
                                                                    aaCGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGTGCTGGc  <  1:2197832/71‑1 (MQ=255)
                                                                    |                                                                      
TATTAGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGTGCTGGC  >  minE/2292549‑2292687

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: