Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,564,887 C→T 63.9% intergenic (‑99/+193) 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,8870CT63.9% 50.2 / 55.5 61intergenic (‑99/+193)yfdE/glkpredicted CoA‑transferase, NAD(P)‑binding/glucokinase
Reads supporting (aligned to +/- strand):  ref base C (14/8);  new base T (14/25);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 6.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGCATTTATGTAAAATCAATAACATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACTG  >  minE/1564831‑1564950
                                                        |                                                               
ggCATTTATGTAAAATCAATAACATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTa                                                   >  1:831183/1‑71 (MQ=255)
ggCATTTATGTAAAATCAATAACATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTa                                                   >  1:370087/1‑71 (MQ=255)
  cATTTATGTAAAATCAATAACATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAAt                                                     <  1:1946752/67‑1 (MQ=255)
      tATGTAAAATCAATAACATGATTAACGTTAATTCATCGTTCAGTGATAGTCGACGAGTGAAATTATGa                                                >  1:3969693/1‑68 (MQ=255)
       aTGTAAAATCAATAACATGATTAACGTTAATTCCTCGTTCAGTGATAGCCGACGAGTGAAATTATGAACg                                             <  1:5423354/70‑1 (MQ=255)
        tGTAAAATCAATAACATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACg                                             <  1:5243747/69‑1 (MQ=255)
        tGTAAAATCAATAACATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACg                                             <  1:3051495/69‑1 (MQ=255)
                  aTAACATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACa                                    <  1:1394095/68‑1 (MQ=255)
                    aaCATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACaaa                                  >  1:1090873/1‑68 (MQ=255)
                     aCATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAAt                                                     >  1:3716286/1‑48 (MQ=255)
                      cATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGcc                             >  1:3983791/1‑71 (MQ=255)
                        tGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCg                               >  1:1351107/1‑67 (MQ=255)
                            tAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTAt                                                  >  1:5340105/1‑44 (MQ=255)
                            tAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTAt                                                  >  1:4620577/1‑44 (MQ=255)
                            tAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTAt                                                  >  1:3882854/1‑44 (MQ=255)
                             atcGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTAc                        >  1:13374/3‑69 (MQ=255)
                             atcGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTAc                        >  1:5665730/3‑69 (MQ=255)
                               cGTTAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGaa                                               >  1:3674909/1‑44 (MQ=255)
                                gTTAAATCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAAt                                 >  1:989525/1‑57 (MQ=255)
                                 ttAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTa                                                   <  1:89866/38‑1 (MQ=255)
                                 ttAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTa                                                   <  1:2742161/38‑1 (MQ=255)
                                 ttAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTa                                                   <  1:2246871/38‑1 (MQ=255)
                                 ttAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTa                                                   <  1:75740/38‑1 (MQ=255)
                                 ttAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTa                                                   <  1:1281568/38‑1 (MQ=255)
                                 ttAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGc                                          <  1:1555840/47‑1 (MQ=255)
                                 ttAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGc                                          <  1:1549224/47‑1 (MQ=255)
                                 ttAATTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTa                         <  1:4799263/64‑1 (MQ=255)
                                 ttAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAAc                                              <  1:3607175/43‑1 (MQ=255)
                                 ttAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGc                              >  1:2518740/1‑59 (MQ=255)
                                   aaTTCCTCGTTCAGTGATAGTTGACGAGTGAAAtt                                                    >  1:5277583/1‑35 (MQ=255)
                                   aaTTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACaaa                >  1:745230/1‑71 (MQ=255)
                                   aaTTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACaaa                >  1:1866222/1‑71 (MQ=255)
                                    aTTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaat               <  1:620533/71‑1 (MQ=255)
                                    aTTCCTCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAaat               <  1:3780505/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)
                                         tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCaa                                        <  1:1514477/41‑2 (MQ=255)
                                         tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                      <  1:1181267/59‑1 (MQ=255)
                                         tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                      <  1:453067/59‑1 (MQ=255)
                                         tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                      <  1:1866768/59‑1 (MQ=255)
                                         tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                      <  1:3649770/59‑1 (MQ=255)
                                         tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCa                      <  1:3466279/59‑1 (MQ=255)
                                         tCGTTCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATg           >  1:1232916/1‑70 (MQ=255)
                                            ttCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCAtt       <  1:5621617/71‑1 (MQ=255)
                                            ttCAGTGATAGTTGACGAGTGAAATTATGAACGTGCAAAACAAATCGCCGTTACCAGACAAATAATGCAtt       <  1:2494953/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:2050402/39‑1 (MQ=38)
                                                 tGATAGTTGACGAGTGAAATTATGAACGTGCAAAACaaa                                  <  1:2188321/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)
                                                        |                                                               
GGCATTTATGTAAAATCAATAACATGATTAACGTTAATTCCTCGTTCAGTGATAGTCGACGAGTGAAATTATGAACGTGCAGAACAAATCGCCGTTACCAGACAAATAATGCATTTACTG  >  minE/1564831‑1564950

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: