Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,247 T→A 32.2% S637S (TCT→TCA thrA → fused aspartokinase I and homoserine dehydrogenase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,2470TA32.2% 63.8 / 43.3 62S637S (TCT→TCAthrAfused aspartokinase I and homoserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base T (10/32);  new base A (3/17);  total (13/49)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

GAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCT  >  minE/2185‑2312
                                                              |                                                                 
gagaACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCg                                                             >  1:2032932/1‑69 (MQ=255)
gagaACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCg                                                             >  1:2773219/1‑69 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:1690082/69‑1 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:2496709/69‑1 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:1456974/69‑1 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:1482304/69‑1 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:240329/69‑1 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:2679190/69‑1 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:501088/69‑1 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:2832265/69‑1 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:2290945/69‑1 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:1637966/69‑1 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:412499/69‑1 (MQ=255)
    aCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                         <  1:1042771/69‑1 (MQ=255)
     ccTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCttt                                                         >  1:2811814/1‑68 (MQ=255)
     ccTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCttt                                                         >  1:412942/1‑68 (MQ=255)
      cTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCt                                                       <  1:2446871/69‑1 (MQ=255)
      cTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCt                                                       <  1:2598539/69‑1 (MQ=255)
                       tGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGt                                       <  1:2262425/68‑1 (MQ=255)
                        gCAGGTGATGAATTGATGAAGTTCTCCGTCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTa                                     <  1:3058158/69‑1 (MQ=255)
                        gCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTa                                     <  1:1702408/69‑1 (MQ=255)
                        gCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTa                                     <  1:2055293/69‑1 (MQ=255)
                        gCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTa                                     <  1:1284694/69‑1 (MQ=255)
                        gCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTa                                     <  1:1255670/69‑1 (MQ=255)
                        gCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTa                                     <  1:2286734/69‑1 (MQ=255)
                        gCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTa                                     <  1:1063607/69‑1 (MQ=255)
                        gCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTa                                     <  1:2059303/69‑1 (MQ=255)
                              gatgaATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGaa                               >  1:2080863/1‑69 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAgg                             <  1:2013351/68‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAgg                             <  1:1858158/68‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAgg                             <  1:2091881/68‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAgg                             <  1:1562029/68‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAgg                             <  1:2659940/68‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAgg                             <  1:2989196/68‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAgg                             <  1:3343197/68‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTCATATCTTCGGCAAGTTAGACGAAgg                             <  1:586148/68‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTCATATCTTCGGCAAGTTAGACGAAgg                             <  1:590453/68‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGc                            <  1:736109/69‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGc                            <  1:3188161/69‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGc                            <  1:1096211/69‑1 (MQ=255)
                                 gaATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGc                            <  1:1667338/69‑1 (MQ=255)
                                  aaTTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAgg                             <  1:3044812/67‑1 (MQ=255)
                                    ttGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                                 <  1:3071360/61‑1 (MQ=255)
                                    ttGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                                 <  1:808126/61‑1 (MQ=255)
                                    ttGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                                 <  1:1363753/61‑1 (MQ=255)
                                    ttGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                                 <  1:1368229/61‑1 (MQ=255)
                                    ttGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                                 <  1:2552125/61‑1 (MQ=255)
                                    ttGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                                 <  1:2235754/61‑1 (MQ=255)
                                    ttGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                                 <  1:756540/61‑1 (MQ=255)
                                    ttGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                                 <  1:1051219/61‑1 (MQ=255)
                                    ttGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                                 <  1:1940671/61‑1 (MQ=255)
                                     tgatgaAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAgg                             <  1:2174529/64‑1 (MQ=255)
                                     tgatgaAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGa                        >  1:1600111/1‑69 (MQ=255)
                                        tgaAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGc                            <  1:75923/62‑1 (MQ=255)
                                         gaAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGttt                    >  1:2105109/1‑69 (MQ=255)
                                         gaAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGttt                    >  1:319780/1‑69 (MQ=255)
                                         gaAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGttt                    >  1:1462296/1‑69 (MQ=255)
                                               ctcCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGa              >  1:310385/1‑69 (MQ=255)
                                                tcCGACATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                                 >  1:1938184/1‑49 (MQ=38)
                                                tcCGACATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGa                                >  1:3115963/1‑50 (MQ=39)
                                                tcCGACATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGa                                >  1:924245/1‑50 (MQ=39)
                                                   ggCATTCTTTCTTGTTCGCTTTCTTATATCTTCGGCAAGTTa                                     <  1:1600797/42‑1 (MQ=255)
                                                       ttctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCa      <  1:892198/69‑1 (MQ=255)
                                                       ttctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCa      <  1:1587165/69‑1 (MQ=255)
                                                       ttctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCa      <  1:1006653/69‑1 (MQ=255)
                                                       ttctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCa      <  1:160987/69‑1 (MQ=255)
                                                       ttctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCa      <  1:1916099/69‑1 (MQ=255)
                                                        tctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCAc     <  1:1249233/69‑1 (MQ=255)
                                                        tctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCAc     <  1:1121264/69‑1 (MQ=255)
                                                        tctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCAc     <  1:2842823/69‑1 (MQ=255)
                                                        tctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCAc     <  1:2265596/69‑1 (MQ=255)
                                                        tctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCAc     <  1:2139047/69‑1 (MQ=255)
                                                        tctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCAc     <  1:107676/69‑1 (MQ=255)
                                                        tctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCAc     <  1:1729126/69‑1 (MQ=255)
                                                        tctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCAc     <  1:1860132/69‑1 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:426810/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:782212/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:51395/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:740060/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:988001/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:1821418/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:1015225/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:1102572/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:1209758/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:1616357/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:1671633/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:3374133/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:1895360/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:1976042/1‑69 (MQ=255)
                                                         ctttctGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACg    >  1:2732983/1‑69 (MQ=255)
                                                              tGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCt  <  1:555341/66‑1 (MQ=255)
                                                              tGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCt  <  1:1690574/66‑1 (MQ=255)
                                                              tGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCt  <  1:2180050/66‑1 (MQ=255)
                                                              tGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCt  <  1:1627679/66‑1 (MQ=255)
                                                              tGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCt  <  1:2394940/66‑1 (MQ=255)
                                                              tGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCt  <  1:1267405/66‑1 (MQ=255)
                                                              tGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCt  <  1:3152492/66‑1 (MQ=255)
                                                              tGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCt  <  1:3291847/66‑1 (MQ=255)
                                                              |                                                                 
GAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCT  >  minE/2185‑2312

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: