Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 428,771 C→T 31.5% V509V (GTC→GTT 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,7710CT31.5% 90.6 / 48.0 98V509V (GTC→GTTnagEfused N‑acetyl glucosamine specific PTS enzyme IICBA components
Reads supporting (aligned to +/- strand):  ref base C (34/33);  new base T (8/23);  total (42/56)
Fisher's exact test for biased strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGA  >  minE/428705‑428838
                                                                  |                                                                   
aaaaCCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGt                                                                   >  1:3411861/1‑69 (MQ=255)
  aaCCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTgg                                                                 <  1:1195414/69‑1 (MQ=255)
   aCCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGc                                                                <  1:1639662/69‑1 (MQ=255)
   aCCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGc                                                                <  1:1406108/69‑1 (MQ=255)
    ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                               >  1:69461/1‑69 (MQ=255)
    ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                               >  1:2028288/1‑69 (MQ=255)
    ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                               >  1:1844493/1‑69 (MQ=255)
    ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                               >  1:3090227/1‑69 (MQ=255)
    ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                               >  1:2573490/1‑69 (MQ=255)
    ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                               >  1:3076864/1‑69 (MQ=255)
    ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                               >  1:1472407/1‑69 (MQ=255)
    ccGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTTGTGGCa                                                               >  1:1424699/1‑69 (MQ=255)
                         gTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGaa                                          <  1:1651306/69‑1 (MQ=255)
                         gTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGaa                                          <  1:2623519/69‑1 (MQ=255)
                          tatctatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAg                                         <  1:3428908/69‑1 (MQ=255)
                          tatctatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAg                                         <  1:1743406/69‑1 (MQ=255)
                           atctatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGc                                        >  1:80604/1‑69 (MQ=255)
                           atctatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGc                                        >  1:1822081/1‑69 (MQ=255)
                            tctatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCa                                       >  1:3039154/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)
                              tatcGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCAtt                                     >  1:2534061/1‑69 (MQ=255)
                                       gCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCaa                            >  1:3230006/1‑69 (MQ=255)
                                       gCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCaa                            >  1:2878167/1‑69 (MQ=255)
                                       gCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCaa                            >  1:1511876/1‑69 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:600832/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:3038573/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:638441/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:2476656/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:2706012/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:2755058/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:235328/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:2962683/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:2264477/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:3025982/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:3373512/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:768865/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:1468954/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:148809/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:1261503/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:1542653/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:838004/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:1838843/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:1260127/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:1788817/69‑1 (MQ=255)
                                          ggTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:1780394/69‑1 (MQ=255)
                                           gTATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:2004967/68‑1 (MQ=255)
                                            tATCGCCGATTACCGGGGATGTTGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGc                         <  1:2722614/67‑1 (MQ=255)
                                             aTCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggt                      >  1:889947/1‑69 (MQ=255)
                                             aTCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggt                      >  1:508937/1‑69 (MQ=255)
                                             aTCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggt                      >  1:327067/1‑69 (MQ=255)
                                             aTCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggt                      >  1:3158763/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:421912/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:3204419/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     >  1:1011402/1‑69 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:3329779/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:3216397/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:3097304/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:3163833/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:2800521/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:1389707/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:2653793/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:254666/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:1437467/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:1821697/69‑1 (MQ=255)
                                              tCGCCGATTACCGGTGATGTCGTGGCACGGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:3122168/69‑1 (MQ=255)
                                                gCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAg                                         <  1:3023591/47‑1 (MQ=255)
                                                gCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAg                                         <  1:1496322/47‑1 (MQ=255)
                                                gCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAg                                         <  1:1773882/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:1277492/1‑69 (MQ=255)
                                                      tACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACgg             >  1:2761731/1‑69 (MQ=255)
                                                              aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtg                     <  1:2185440/53‑1 (MQ=255)
                                                              aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCg      >  1:874842/1‑68 (MQ=255)
                                                              aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCg      >  1:702165/1‑68 (MQ=255)
                                                              aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCg      >  1:984030/1‑68 (MQ=255)
                                                              aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCg      >  1:495362/1‑68 (MQ=255)
                                                              aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCg      >  1:308349/1‑68 (MQ=255)
                                                              aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCg      >  1:2204574/1‑68 (MQ=255)
                                                              aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCg      >  1:2108090/1‑68 (MQ=255)
                                                              aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCg      >  1:1366213/1‑68 (MQ=255)
                                                                 tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGa  >  1:52428/1‑69 (MQ=255)
                                                                 tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGa  >  1:2894148/1‑69 (MQ=255)
                                                                 tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGa  <  1:2798970/69‑1 (MQ=255)
                                                                 tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGa  <  1:251283/69‑1 (MQ=255)
                                                                 tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGa  >  1:2125819/1‑69 (MQ=255)
                                                                 tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGa  >  1:1971772/1‑69 (MQ=255)
                                                                 tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGa  >  1:1952872/1‑69 (MQ=255)
                                                                 tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGa  <  1:1804227/69‑1 (MQ=255)
                                                                 tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGa  >  1:1095392/1‑69 (MQ=255)
                                                                  |                                                                   
AAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGA  >  minE/428705‑428838

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: