Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,395,087 T→C 10.4% intergenic (‑160/‑208) fruB ← / → setB fused fructose‑specific PTS enzyme IIA component and HPr component/lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,395,0870TC10.4% 202.3 / 12.9 77intergenic (‑160/‑208)fruB/setBfused fructose‑specific PTS enzyme IIA component and HPr component/lactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base T (30/39);  new base C (3/5);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

AACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATCCGTTATCAGCAATTTC  >  minE/1395019‑1395156
                                                                    |                                                                     
aaCGTTTCAAGAAAGCATAATACTTCTTTTTCAGAACTCAAGGAAAGCCGTAAAGTTATTAGCAAAAGt                                                                       <  1:1382931/69‑1 (MQ=255)
aaCGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGttt                                                                     <  1:1419358/71‑1 (MQ=255)
aaCGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGttt                                                                     <  1:711577/71‑1 (MQ=255)
aaCGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGt                                                                       >  1:3289577/1‑69 (MQ=255)
aaCGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGctt                                                                     >  1:2428467/1‑71 (MQ=255)
  cGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGa                                                                   >  1:1607243/1‑71 (MQ=255)
    tttCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGAtg                                                                 <  1:2750746/71‑1 (MQ=255)
    tttCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGCTTGAtg                                                                 <  1:1232199/71‑1 (MQ=255)
      tCAAGAAAGGATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGAtgt                                                                <  1:3794849/70‑1 (MQ=255)
      tCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGAtgt                                                                <  1:122931/70‑1 (MQ=255)
       cAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGAtgtg                                                               <  1:2710319/70‑1 (MQ=255)
       cAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGAtgtg                                                               <  1:173493/70‑1 (MQ=255)
       cAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGAtgtg                                                               <  1:1957565/70‑1 (MQ=255)
           aaaGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGca                                                          >  1:571165/1‑71 (MQ=255)
           aaaGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGca                                                          >  1:1746735/1‑71 (MQ=255)
            aaGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGca                                                          >  1:2213229/1‑70 (MQ=255)
             aGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGcaca                                                        <  1:370921/71‑1 (MQ=255)
              gCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGtgctg                                                            <  1:1746903/66‑1 (MQ=255)
              gCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGtgctg                                                            <  1:466081/66‑1 (MQ=255)
              gCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACaa                                                       <  1:2745000/71‑1 (MQ=255)
                ataataCTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGttt                                                                     <  1:1525139/55‑1 (MQ=255)
                 taataCTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAAttt                                                    <  1:909375/71‑1 (MQ=255)
                   ataCTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAAttttt                                                  <  1:3029451/71‑1 (MQ=255)
                   ataCTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAAtttt                                                   <  1:516555/70‑1 (MQ=255)
                   ataCTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAAtttt                                                   <  1:2944210/70‑1 (MQ=255)
                   ataCTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAAtttt                                                   <  1:2163676/70‑1 (MQ=255)
                         cTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAt                                            >  1:3686559/1‑71 (MQ=255)
                         cTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAt                                            >  1:2053216/1‑71 (MQ=255)
                         cTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAt                                            >  1:2852826/1‑71 (MQ=255)
                            tttCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGtt                                                                      <  1:2507216/42‑1 (MQ=255)
                            tttCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAtttt                                         <  1:388545/71‑1 (MQ=255)
                            tttCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAtttt                                         <  1:2253915/71‑1 (MQ=255)
                            tttCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAttt                                          <  1:864593/70‑1 (MQ=255)
                            tttCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAttt                                          <  1:3459800/70‑1 (MQ=255)
                            tttCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAt                                            <  1:2546294/68‑1 (MQ=255)
                            tttCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGCTTGATGTGCTGCACAATTTTTTGCCAtttt                                         <  1:207532/71‑1 (MQ=255)
                              tCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGt                                                                       >  1:1940913/1‑39 (MQ=255)
                              tCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGt                                                                       >  1:2169696/1‑39 (MQ=255)
                              tCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAttttt                                        >  1:2261095/1‑70 (MQ=255)
                                aGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAttttt                                        >  1:2200991/1‑68 (MQ=255)
                                aGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAttttt                                        >  1:3487482/1‑68 (MQ=255)
                                aGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCt                                     <  1:3313871/71‑1 (MQ=255)
                                aGCACGAAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTg                                                                    <  1:2690718/40‑1 (MQ=255)
                                   aCGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGAt                                                                  >  1:1284425/1‑39 (MQ=255)
                                   aCGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAAt                                  <  1:3370218/71‑1 (MQ=255)
                                     gCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTg                                <  1:2642057/71‑1 (MQ=255)
                                     gCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTg                                <  1:489728/71‑1 (MQ=255)
                                     gCAAGGAAAGCCGTAAAGTTATGAGCAAAAGCTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTa                                    <  1:1207018/67‑1 (MQ=255)
                                     gCAAGGAAAGCCGTAAAGTTATGAGCAAAAGCTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTa                                    <  1:1667474/67‑1 (MQ=255)
                                       aaGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCt                              >  1:355761/1‑71 (MQ=255)
                                       aaGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCt                              >  1:3619652/1‑71 (MQ=255)
                                       aaGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCt                              >  1:1083969/1‑71 (MQ=255)
                                         ggAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGa                            <  1:2359837/71‑1 (MQ=255)
                                         ggAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGa                            <  1:1263380/71‑1 (MQ=255)
                                         ggAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGa                            <  1:3896336/71‑1 (MQ=255)
                                         ggAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGa                            <  1:955257/71‑1 (MQ=255)
                                         ggAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGa                            <  1:3261951/71‑1 (MQ=255)
                                          gAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAttttt                                        >  1:3408676/1‑58 (MQ=255)
                                             aGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATgg                         >  1:290561/1‑70 (MQ=255)
                                              gCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGa                       <  1:1767352/71‑1 (MQ=255)
                                                cGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAAttt                                                    <  1:2278554/40‑1 (MQ=255)
                                                  tAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTg                             >  1:2076586/1‑61 (MQ=255)
                                                  tAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTg                             >  1:3629012/1‑61 (MQ=255)
                                                  tAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAAt                   <  1:3557707/71‑1 (MQ=255)
                                                   aaaGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATc                  >  1:1788772/1‑71 (MQ=255)
                                                   aaaGTTATGAGCAAAAGCTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATc                  >  1:2155804/1‑71 (MQ=255)
                                                    aaGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAtttt                                         >  1:3223986/1‑47 (MQ=255)
                                                    aaGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAtttt                                         >  1:3219307/1‑47 (MQ=255)
                                                    aaGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATgg                         >  1:2672995/1‑63 (MQ=255)
                                                    aaGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATcc                 >  1:714642/1‑71 (MQ=255)
                                                    aaGTTATGAGCAAAAGCTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATcc                 >  1:2966481/1‑71 (MQ=255)
                                                    aaGTTATGAGCAAAAGCTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATcc                 >  1:3118790/1‑71 (MQ=255)
                                                       ttATGAGCTAAAGCTTGATGTGCTGCACAATTTTTTGCCAttttt                                        <  1:188187/45‑1 (MQ=38)
                                                         aTGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATcc                 >  1:2533812/1‑66 (MQ=255)
                                                           gAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTa                                    >  1:3638579/1‑45 (MQ=255)
                                                            aGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATCCGTTATCa          <  1:1029737/70‑1 (MQ=255)
                                                             gCAAAAGTTTGATGTGCTGCACAATTTTTTGCCAttttt                                        >  1:2668680/1‑39 (MQ=255)
                                                               aaaaGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAAtt                                 >  1:1254254/1‑44 (MQ=255)
                                                               aaaaGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGc                               <  1:1522185/46‑1 (MQ=255)
                                                                aaaGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATCCGTTATCAGCAat     >  1:2885633/1‑71 (MQ=255)
                                                                aaaGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATCCGTTATCAGCAat     >  1:2124387/1‑71 (MQ=255)
                                                                aaaGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATCCGTTATCAGCAat     <  1:2395465/71‑1 (MQ=255)
                                                                 aaGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATCCGTTATCAGCAatt    >  1:1474229/1‑71 (MQ=255)
                                                                 aaGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATCCGTTATCAGCAatt    <  1:2512450/71‑1 (MQ=255)
                                                                  aGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATCCGTTATCAGCAattt   <  1:988017/71‑1 (MQ=255)
                                                                   gTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATCCGTTATCAGCAatttc  <  1:2516100/71‑1 (MQ=255)
                                                                    |                                                                     
AACGTTTCAAGAAAGCATAATACTTCTGTTTCAGCACGCAAGGAAAGCCGTAAAGTTATGAGCAAAAGTTTGATGTGCTGCACAATTTTTTGCCATTTTTCCTAATTGCTGATGGGAAAATCCGTTATCAGCAATTTC  >  minE/1395019‑1395156

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: