Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,906,647 T→C 23.0% A306A (GCA→GCG treB ← fused trehalose(maltose)‑specific enzyme IIBC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,906,6470TC23.0% 97.7 / 24.3 78A306A (GCA→GCGtreBfused trehalose(maltose)‑specific enzyme IIBC component of PTS
Reads supporting (aligned to +/- strand):  ref base T (36/24);  new base C (4/14);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 6.78e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGCC  >  minE/2906579‑2906709
                                                                    |                                                              
tCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAAtgct                                                                >  1:3177116/1‑69 (MQ=255)
tCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAAtgct                                                                >  1:1789382/1‑69 (MQ=255)
tCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAAtgct                                                                >  1:890681/1‑69 (MQ=255)
          gCGTGGTCTGGTGTACACCGGTGATCACCAGCGTGGCGTACATGAAGCCAAACAATGCTGCGCCAATCg                                                      <  1:953284/69‑1 (MQ=255)
          gCGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCg                                                      <  1:966654/69‑1 (MQ=255)
          gCGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCg                                                      <  1:2634488/69‑1 (MQ=255)
           cGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCgg                                                     >  1:1711542/1‑69 (MQ=255)
           cGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCgg                                                     >  1:603500/1‑69 (MQ=255)
           cGTGGTCTGGTGTACACCGGTGATAACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCgg                                                     >  1:417143/1‑69 (MQ=255)
               gTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAg                                                   >  1:3348724/1‑67 (MQ=255)
               gTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAg                                                   >  1:506675/1‑67 (MQ=255)
               gTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAg                                                   >  1:432859/1‑67 (MQ=255)
               gTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAg                                                   >  1:1937784/1‑67 (MQ=255)
               gTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAg                                                   >  1:1008097/1‑67 (MQ=255)
               gTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAg                                                   >  1:2805770/1‑67 (MQ=255)
               gTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAg                                                   >  1:3113945/1‑67 (MQ=255)
               gTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAg                                                   >  1:308921/1‑67 (MQ=255)
                      gtACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCt                                            <  1:3376163/67‑1 (MQ=255)
                             ggTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGtc                                     >  1:264242/1‑67 (MQ=255)
                             ggTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGtc                                     >  1:994829/1‑67 (MQ=255)
                             ggTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGtc                                     >  1:590032/1‑67 (MQ=255)
                             ggTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGtc                                     >  1:599409/1‑67 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGt                              <  1:1490016/68‑1 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGt                              <  1:1762737/68‑1 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGt                              <  1:2026467/68‑1 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGt                              <  1:2072943/68‑1 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGt                              <  1:3407885/68‑1 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGt                              <  1:249509/68‑1 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTg                             >  1:3198778/1‑69 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTg                             >  1:1398832/1‑69 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTg                             >  1:2600802/1‑69 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTg                             >  1:3343049/1‑69 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTg                             >  1:3215767/1‑69 (MQ=255)
                                   cACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTg                             >  1:1506518/1‑69 (MQ=255)
                                    aCCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGt                              <  1:1440602/67‑1 (MQ=255)
                                    aCCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGa                            <  1:1645359/69‑1 (MQ=255)
                                    aCCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGa                            <  1:3079214/69‑1 (MQ=255)
                                    aCCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGa                            <  1:63899/69‑1 (MQ=255)
                                    aCCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGa                            <  1:2189647/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:2915679/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:893299/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:959150/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:564453/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:1054637/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:530609/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:3066002/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:3039010/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:2907837/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:1209650/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:1214494/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:256711/69‑1 (MQ=255)
                                             gCGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCAGTCATCAGGTGACGTACCGCa                   <  1:808697/69‑1 (MQ=255)
                                              cGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:3342616/68‑1 (MQ=255)
                                              cGTACAGGAAGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:3050690/68‑1 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCaaa                                               >  1:734632/1‑36 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGtcatca                                 >  1:1466364/1‑50 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCaa                  >  1:2320379/1‑65 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:84123/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:1475674/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:1473963/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:758974/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:165020/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:1452554/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:1219825/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:1153382/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:3056889/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              <  1:1720239/69‑1 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:2391740/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:2394054/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:2474775/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:3294594/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:3282145/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:2841645/1‑69 (MQ=255)
                                                  cAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc              >  1:3073238/1‑69 (MQ=255)
                                                       aGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGc         <  1:1459041/69‑1 (MQ=255)
                                                       aGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGc         <  1:2718517/69‑1 (MQ=255)
                                                       aGCCAAACAATGCCGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                   <  1:2856040/59‑1 (MQ=255)
                                                        gCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGcc        <  1:817457/69‑1 (MQ=255)
                                                              cAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGcc  <  1:1542620/69‑1 (MQ=255)
                                                              cAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGcc  <  1:1464047/69‑1 (MQ=255)
                                                              cAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGcc  <  1:1313146/69‑1 (MQ=255)
                                                                    |                                                              
TCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGCC  >  minE/2906579‑2906709

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: