Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,391,735 T→C 24.7% K370E (AAA→GAA)  fruA ← fused fructose‑specific PTS enzyme IIB'BC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,391,7350TC24.7% 128.0 / 46.4 73K370E (AAA→GAA) fruAfused fructose‑specific PTS enzyme IIB'BC components
Reads supporting (aligned to +/- strand):  ref base T (26/29);  new base C (7/11);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGATCGGTTT  >  minE/1391666‑1391804
                                                                     |                                                                     
gAACCTCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTa                                                                      <  1:3080924/71‑1 (MQ=255)
   ccGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTAcc                                                                    >  1:3699283/1‑70 (MQ=255)
   ccGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCg                                                                   <  1:384437/71‑1 (MQ=255)
   ccGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTCACCg                                                                   <  1:443106/71‑1 (MQ=255)
    cGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTCACCGa                                                                  <  1:1428275/71‑1 (MQ=255)
        ttCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTCTCCGATCAg                                                              <  1:1052977/71‑1 (MQ=255)
         tCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAgg                                                             >  1:3176837/1‑71 (MQ=255)
          cGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGAt                                                                 >  1:167157/1‑66 (MQ=255)
          cGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAgg                                                             >  1:1370503/1‑70 (MQ=255)
           gCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTa                                                           <  1:1146796/71‑1 (MQ=255)
            cAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTCACCGATCAGGTa                                                           <  1:1939825/70‑1 (MQ=255)
             aGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTCACCGATCAGGTAGa                                                         <  1:204521/71‑1 (MQ=255)
              gTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGat                                                        <  1:164556/71‑1 (MQ=255)
               tCCCCATGGGCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGatc                                                       <  1:926662/71‑1 (MQ=255)
                 cccATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCg                                                                   <  1:598093/57‑1 (MQ=255)
                 cccATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGatcat                                                     <  1:205375/71‑1 (MQ=255)
                    aTGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGc                                                  >  1:3647111/1‑71 (MQ=255)
                       gTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATAGc                                                  >  1:496278/1‑68 (MQ=255)
                       gTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTCACCGATCAGGTAGATCATCGc                                                  >  1:2494588/1‑68 (MQ=255)
                        tCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAga                                              >  1:2779581/1‑71 (MQ=255)
                        tCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAga                                              >  1:1117237/1‑71 (MQ=255)
                         cTGCGGCCAGTGAGTCAGCCCTTCGAGAATGCCCGCAACTGGTTCACCGATCAGGTAGATCATCGCCAgac                                             >  1:3817218/1‑71 (MQ=255)
                          tGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTCACCGATCAGGTAGatc                                                       <  1:2101642/60‑1 (MQ=255)
                            cAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTa                                                                      >  1:3670439/1‑43 (MQ=255)
                            cAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGt                                                            <  1:3427440/53‑1 (MQ=255)
                             aGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAgaccgac                                         <  1:4002423/71‑1 (MQ=255)
                               ccAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGacca                                       >  1:3290093/1‑71 (MQ=255)
                               ccAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGacca                                       >  1:3285513/1‑71 (MQ=255)
                               ccAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGacca                                       >  1:2967843/1‑71 (MQ=255)
                                  gTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGAccacca                                    <  1:1398400/71‑1 (MQ=255)
                                  gTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGAccacca                                    <  1:619701/71‑1 (MQ=255)
                                  gTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGAccacca                                    <  1:3187719/71‑1 (MQ=255)
                                   tGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGaccac                                      >  1:2595212/1‑68 (MQ=255)
                                    gAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTCACCGATCAGGTAGat                                                        >  1:1838699/1‑49 (MQ=255)
                                    gAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTCACCGATCAGGTAGat                                                        >  1:3189556/1‑49 (MQ=255)
                                    gAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTCACCGATCAGGTAGa                                                         >  1:3169146/1‑48 (MQ=255)
                                         aGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGa                             <  1:2900723/71‑1 (MQ=255)
                                         aGCCCTTCGAGAATGCCAGCAACTGGTTCACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTgg                              >  1:2625078/1‑70 (MQ=255)
                                         aGCCCTTCGAGAATGCCAGCAACTGGTTCACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGa                             <  1:2373791/71‑1 (MQ=255)
                                          gCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGaa                            <  1:1637135/71‑1 (MQ=255)
                                             cTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTa                                                           <  1:3597245/37‑1 (MQ=255)
                                             cTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAg                                                          <  1:2617439/38‑1 (MQ=255)
                                             cTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAgaccgac                                         <  1:1841050/55‑1 (MQ=255)
                                             cTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGaa                            <  1:308102/68‑1 (MQ=255)
                                             cTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAAtt                         <  1:1933753/71‑1 (MQ=255)
                                               tCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGa                                                         <  1:3265583/37‑1 (MQ=255)
                                               tCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAg                       >  1:3061802/1‑71 (MQ=255)
                                                cGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGaaa                           >  1:1826567/1‑66 (MQ=255)
                                                cGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAg                       >  1:2138020/1‑70 (MQ=255)
                                                cGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGc                      >  1:1561371/1‑71 (MQ=255)
                                                cGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGc                      >  1:2984885/1‑71 (MQ=255)
                                                cGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGc                      >  1:1553748/1‑71 (MQ=255)
                                                cGAGAATGCCAGCAACTGGTTCACCGATCAGGTAGatcat                                                     >  1:598690/1‑40 (MQ=255)
                                                cGAGAATGCCAGCAACTGGTTCACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAg                       <  1:337157/70‑1 (MQ=255)
                                                  agaATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCgg                    <  1:1055288/71‑1 (MQ=255)
                                                  agaATGCCAGCAACTGGTTCACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCgg                    <  1:4143893/71‑1 (MQ=255)
                                                   gaATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAgaccgac                                         >  1:2940641/1‑49 (MQ=255)
                                                   gaATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGaccac                                      <  1:2139990/52‑1 (MQ=255)
                                                     aTGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAAt                 >  1:475509/1‑71 (MQ=255)
                                                      tGCCAGCAACTGGTTTACCTATCAGGTAGATCATCGCCAgaccgac                                         <  1:3737732/46‑1 (MQ=255)
                                                      tGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGc                                                  <  1:2013622/37‑1 (MQ=255)
                                                      tGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGaccac                                      <  1:2990845/49‑1 (MQ=255)
                                                      tGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGaccac                                      <  1:467238/49‑1 (MQ=255)
                                                      tGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTgg                              >  1:781664/1‑57 (MQ=255)
                                                       gCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCa                                                >  1:1691233/1‑38 (MQ=255)
                                                        ccAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGAt              >  1:2943981/1‑71 (MQ=255)
                                                        ccAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGAt              >  1:1391918/1‑71 (MQ=255)
                                                         cAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATc             <  1:2576145/71‑1 (MQ=255)
                                                               cTGGTTTACCGATCAGGTAGATCATCGCCAgaccg                                           <  1:348640/35‑1 (MQ=255)
                                                                tGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAgg          <  1:2025912/67‑1 (MQ=255)
                                                                  gTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGATCGGt    >  1:3851760/1‑71 (MQ=255)
                                                                  gTTCACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGATCg      <  1:4065136/69‑1 (MQ=255)
                                                                  gTTCACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGATCg      <  1:2181933/69‑1 (MQ=255)
                                                                     tACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGATCGGttt  <  1:1499428/70‑1 (MQ=255)
                                                                     |                                                                     
GAACCGCATTCGCAGTCCCCATGGTCTGCAGCCAGTGAGTCAGCCCTTCGAGAATGCCAGCAACTGGTTTACCGATCAGGTAGATCATCGCCAGACCGACCACCAGACTGGAAATTAGCGGAATGATCAGGATCGGTTT  >  minE/1391666‑1391804

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: