Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 428,765 T→G 38.1% G507G (GGT→GGG nagE → fused N‑acetyl glucosamine specific PTS enzyme IICBA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE428,7650TG38.1% 42.4 / 75.1 84G507G (GGT→GGGnagEfused N‑acetyl glucosamine specific PTS enzyme IICBA components
Reads supporting (aligned to +/- strand):  ref base T (20/32);  new base G (8/24);  total (28/56)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCTGTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGG  >  minE/428697‑428827
                                                                    |                                                              
ccTGTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATTGCGGAGCTGGTATCGCCGATTACCGGt                                                                <  1:1852224/69‑1 (MQ=255)
ccTGTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGt                                                                <  1:194438/69‑1 (MQ=255)
ccTGTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGt                                                                <  1:1924517/69‑1 (MQ=255)
ccTGTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGt                                                                <  1:1273939/69‑1 (MQ=255)
ccTGTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGt                                                                <  1:288435/69‑1 (MQ=255)
ccTGTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGt                                                                <  1:2325075/69‑1 (MQ=255)
ccTGTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGt                                                                <  1:3375814/69‑1 (MQ=255)
   gTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGAt                                                             <  1:479925/69‑1 (MQ=255)
   gTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGAt                                                             <  1:1818669/69‑1 (MQ=255)
    tAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGAt                                                             <  1:683585/68‑1 (MQ=255)
      gCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGt                                                           <  1:1540313/68‑1 (MQ=255)
        aaaaCCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGt                                                        >  1:3411861/1‑69 (MQ=255)
          aaCCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTgg                                                      <  1:1195414/69‑1 (MQ=255)
           aCCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGc                                                     <  1:1406108/69‑1 (MQ=255)
           aCCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGc                                                     <  1:1639662/69‑1 (MQ=255)
            ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                    >  1:1844493/1‑69 (MQ=255)
            ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                    >  1:2028288/1‑69 (MQ=255)
            ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                    >  1:69461/1‑69 (MQ=255)
            ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                    >  1:2573490/1‑69 (MQ=255)
            ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                    >  1:1472407/1‑69 (MQ=255)
            ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                    >  1:1424699/1‑69 (MQ=255)
            ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                    >  1:3076864/1‑69 (MQ=255)
            ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                    >  1:3090227/1‑69 (MQ=255)
                                 gTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGaa                               <  1:1651306/69‑1 (MQ=255)
                                 gTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGaa                               <  1:2623519/69‑1 (MQ=255)
                                  tatctatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAg                              <  1:1743406/69‑1 (MQ=255)
                                  tatctatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAg                              <  1:3428908/69‑1 (MQ=255)
                                   atctatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGc                             >  1:1822081/1‑69 (MQ=255)
                                   atctatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGc                             >  1:80604/1‑69 (MQ=255)
                                    tctatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCa                            >  1:3039154/1‑69 (MQ=255)
                                      tatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCAtt                          >  1:2534061/1‑69 (MQ=255)
                                      tatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCAtt                          >  1:566748/1‑69 (MQ=255)
                                      tatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCAtt                          >  1:1999442/1‑69 (MQ=255)
                                      tatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCAtt                          >  1:1117814/1‑69 (MQ=255)
                                               gCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCaa                 >  1:2878167/1‑69 (MQ=255)
                                               gCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCaa                 >  1:3230006/1‑69 (MQ=255)
                                               gCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCaa                 >  1:1511876/1‑69 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:768865/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:1260127/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:1261503/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:3038573/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:3025982/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:2962683/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:838004/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:2755058/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:2706012/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:3373512/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:1468954/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:148809/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:1542653/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:600832/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:2476656/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:1838843/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:1780394/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:1788817/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:235328/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:638441/69‑1 (MQ=255)
                                                  ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:2264477/69‑1 (MQ=255)
                                                   gTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:2004967/68‑1 (MQ=255)
                                                    tATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc              <  1:2722614/67‑1 (MQ=255)
                                                     aTCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggt           >  1:508937/1‑69 (MQ=255)
                                                     aTCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggt           >  1:3158763/1‑69 (MQ=255)
                                                     aTCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggt           >  1:327067/1‑69 (MQ=255)
                                                     aTCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggt           >  1:889947/1‑69 (MQ=255)
                                                     aTCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggt           >  1:2073616/1‑69 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGTATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:2471867/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtt          <  1:2218578/69‑2 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:3329779/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:421912/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:1821697/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:254666/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:2653793/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:1437467/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:1389707/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:2800521/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:3216397/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:3204419/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:3163833/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          >  1:1011402/1‑69 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:3097304/69‑1 (MQ=255)
                                                      tCGCCGATTACCGGTGATGTCGTGGCACGGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:3122168/69‑1 (MQ=255)
                                                        gCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAg                              <  1:1773882/47‑1 (MQ=255)
                                                        gCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAg                              <  1:1496322/47‑1 (MQ=255)
                                                        gCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAg                              <  1:3023591/47‑1 (MQ=255)
                                                        gCCGATTACCGGTGATGTCGTGGCACTGGATCAAGTTCCTGACGAAg                              <  1:2360906/47‑1 (MQ=255)
                                                        gCCGATTACCGGGGATGTCGTGGCACTGGATCAGGTTCCTGACGAAg                              <  1:469441/47‑1 (MQ=255)
                                                              tACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg          <  1:3429246/61‑1 (MQ=255)
                                                              tACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACgg  <  1:2296458/69‑1 (MQ=255)
                                                              tACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACgg  >  1:2683384/1‑69 (MQ=255)
                                                              tACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACgg  >  1:2761731/1‑69 (MQ=255)
                                                              tACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACgg  >  1:1277492/1‑69 (MQ=255)
                                                                    |                                                              
CCTGTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGG  >  minE/428697‑428827

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: