Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,564,912 G→A 76.3% intergenic (‑124/+168) yfdE ← / ← glk predicted CoA‑transferase, NAD(P)‑binding/glucokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,564,9120GA76.3% 130.0 / 39.2 76intergenic (‑124/+168)yfdE/glkpredicted CoA‑transferase, NAD(P)‑binding/glucokinase
Reads supporting (aligned to +/- strand):  ref base G (11/7);  new base A (29/29);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAACATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGT  >  minE/1564849‑1564978
                                                               |                                                                  
aTAACATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACa                                                                <  1:1394095/68‑1 (MQ=255)
  aaCATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACaaa                                                              >  1:1090873/1‑68 (MQ=255)
    cATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGcc                                                         >  1:3983791/1‑71 (MQ=255)
      tGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCg                                                           >  1:1351107/1‑67 (MQ=255)
           atcGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTAc                                                    >  1:5665730/3‑69 (MQ=255)
           atcGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTAc                                                    >  1:13374/3‑69 (MQ=255)
              gTTAAATCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAAt                                                             >  1:989525/1‑57 (MQ=255)
               ttAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTa                                                     <  1:4799263/64‑1 (MQ=255)
               ttAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGc                                                          >  1:2518740/1‑59 (MQ=255)
                 aaTTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACaaa                                            >  1:745230/1‑71 (MQ=255)
                 aaTTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACaaa                                            >  1:1866222/1‑71 (MQ=255)
                  aTTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaat                                           <  1:3780505/71‑1 (MQ=255)
                  aTTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaat                                           <  1:620533/71‑1 (MQ=255)
                  aTTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaat                                           <  1:4329885/71‑1 (MQ=255)
                  aTTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaat                                           <  1:3513230/71‑1 (MQ=255)
                     ccTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTa                                                     >  1:3833993/1‑58 (MQ=255)
                     ccTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACaa                                             <  1:5404833/66‑1 (MQ=255)
                     ccTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaataat                                        >  1:3860437/1‑71 (MQ=255)
                     ccTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaataat                                        >  1:1453044/1‑71 (MQ=255)
                     ccTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaataa                                         >  1:3845153/1‑70 (MQ=255)
                     ccTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaataa                                         <  1:4794704/70‑1 (MQ=255)
                     ccTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACAaataat                                        >  1:5530491/1‑71 (MQ=255)
                       tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                                                  <  1:453067/59‑1 (MQ=255)
                       tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                                                  <  1:3466279/59‑1 (MQ=255)
                       tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                                                  <  1:3649770/59‑1 (MQ=255)
                       tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                                                  <  1:1866768/59‑1 (MQ=255)
                       tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                                                  <  1:1181267/59‑1 (MQ=255)
                       tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATg                                       >  1:1232916/1‑70 (MQ=255)
                          ttCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCAtt                                   <  1:2494953/71‑1 (MQ=255)
                          ttCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCAtt                                   <  1:5621617/71‑1 (MQ=255)
                          ttCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCAtt                                   <  1:433230/71‑1 (MQ=255)
                             aGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTa                                 <  1:2541359/70‑1 (MQ=255)
                              gTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACt                               <  1:4275528/71‑1 (MQ=255)
                               tGATAGTTGACGAGTGAAATTATGAACGTGCAAAACaaa                                                              <  1:2188321/39‑1 (MQ=38)
                               tGATAGTTGACGAGTGAAATTATGAACGTGCAAAACaaa                                                              <  1:2050402/39‑1 (MQ=38)
                                gATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACTg                              >  1:4525058/1‑70 (MQ=255)
                                 aTAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGt                                                       >  1:4576650/1‑44 (MQ=38)
                                 aTAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaataa                                         >  1:1927597/1‑58 (MQ=255)
                                 aTAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTg                              >  1:1972136/1‑69 (MQ=255)
                                    gttgacgaGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGgaga                          <  1:210699/67‑1 (MQ=255)
                                    gttgacgaGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGgaga                          <  1:2779091/67‑1 (MQ=255)
                                      tgacgaGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAAt                        >  1:3404159/2‑70 (MQ=255)
                                       gacgaGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCAttt                                  <  1:5086357/59‑1 (MQ=255)
                                       gacgaGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAAtt                       >  1:2494595/1‑70 (MQ=255)
                                       gacgaGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTa                      <  1:3564703/71‑1 (MQ=255)
                                       gacgaGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAAtt                       >  1:4805922/1‑70 (MQ=255)
                                        acgaGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTaa                     >  1:2720629/1‑71 (MQ=255)
                                        acgaGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTaa                     <  1:2357501/71‑1 (MQ=255)
                                        acgaGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTaa                     <  1:934415/71‑1 (MQ=255)
                                         ctaGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGgaga                          <  1:2761392/63‑1 (MQ=255)
                                             tGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaata                                          >  1:5284013/1‑45 (MQ=255)
                                             tGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCg                >  1:5608278/1‑71 (MQ=255)
                                              gAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGc               >  1:3449051/1‑71 (MQ=255)
                                               aaaTTATGAACGTGCAGAACAAATCGCCGTTACCAGACAaat                                           >  1:238978/1‑42 (MQ=255)
                                               aaaTTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCt              >  1:1245752/1‑71 (MQ=255)
                                               aaaTTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaataa                                         >  1:5503957/1‑44 (MQ=255)
                                               aaaTTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTa                                 >  1:1244948/1‑52 (MQ=255)
                                                 aTTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGgaga                          <  1:4067804/57‑1 (MQ=255)
                                                 aTTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaata                                          >  1:2916374/1‑41 (MQ=255)
                                                 aTTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaata                                          >  1:3826324/1‑41 (MQ=255)
                                                 aTTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTg             >  1:4397406/1‑70 (MQ=255)
                                                   tATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTAc                                <  1:4301292/49‑1 (MQ=255)
                                                   tATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTAc                                <  1:5383225/49‑1 (MQ=255)
                                                   tATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTAc                                <  1:4165229/49‑1 (MQ=255)
                                                    aTGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTa                                 >  1:3071688/1‑47 (MQ=255)
                                                    aTGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTa                                 >  1:3426928/1‑47 (MQ=255)
                                                    aTGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCt              >  1:906641/1‑66 (MQ=255)
                                                    aTGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCt              >  1:2168664/1‑66 (MQ=255)
                                                    aTGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCt              >  1:3120076/1‑66 (MQ=255)
                                                     tGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTg                              >  1:4545645/1‑49 (MQ=255)
                                                     tGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACgg        <  1:5235651/71‑1 (MQ=255)
                                                      gAACGTGCAAAACAAATCGCCGTTACCAGACAaataa                                         <  1:337291/37‑1 (MQ=255)
                                                        aCGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAg     >  1:5247916/1‑71 (MQ=255)
                                                        aCGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAg     <  1:678046/71‑1 (MQ=255)
                                                            gCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAgtgt  <  1:4726409/70‑1 (MQ=255)
                                                            gCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAgtgt  <  1:3085498/70‑1 (MQ=255)
                                                               |                                                                  
ATAACATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGT  >  minE/1564849‑1564978

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: