Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,790,053 C→G 22.5% N366K (AAC→AAG pgi → glucosephosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,790,0530CG22.5% 226.1 / 60.3 129N366K (AAC→AAGpgiglucosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (70/30);  new base G (19/10);  total (89/40)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

Only 100 of 131 total aligned reads displayed.

CGGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGCACG  >  minE/2789992‑2790114
                                                             |                                                             
cggcgTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         <  1:2429319/68‑1 (MQ=255)
cggcgTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         <  1:2843537/68‑1 (MQ=255)
cggcgTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         <  1:1736965/68‑1 (MQ=255)
cggcgTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         <  1:2874315/68‑1 (MQ=255)
cggcgTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTCACCGTAACGGTAACGTtgtg                                                         <  1:3127556/68‑1 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:437766/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:2338337/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:3026920/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:2546622/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:2212610/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:203365/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:2604518/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:2931778/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:3432650/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:556847/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:924304/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:629055/1‑61 (MQ=255)
       cTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTtgtg                                                         >  1:278278/1‑61 (MQ=255)
         tCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGa                                               >  1:2773103/1‑69 (MQ=255)
         tCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGa                                               >  1:295646/1‑69 (MQ=255)
         tCCAGCAGGGCAATATGGAATCCAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGa                                               <  1:2858469/69‑1 (MQ=255)
         tCCAGCAGGGCAATATGGAATCCAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGa                                               <  1:3050684/69‑1 (MQ=255)
         tCCAGCAGGGCAATATGGAATCCAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGa                                               <  1:746418/69‑1 (MQ=255)
         tCCAGCAGGGCAATATGGAATCCAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGa                                               <  1:1299835/69‑1 (MQ=255)
         tCCAGCAGGGCAATATGGAATCCAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGa                                               <  1:3395479/69‑1 (MQ=255)
         tCCAGCAGGGCAATATGGAATCCAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGa                                               <  1:1450207/69‑1 (MQ=255)
         tCCAGCAGGGCAATATGGAATCCAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGa                                               <  1:460468/69‑1 (MQ=255)
         tCCAGCAGGGCAATATGGAATCCAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGa                                               <  1:160873/69‑1 (MQ=255)
          ccAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGAc                                              >  1:3239239/1‑69 (MQ=255)
          ccAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGAc                                              >  1:2993700/1‑69 (MQ=255)
                  gCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCg                                       >  1:605388/1‑68 (MQ=255)
                          gAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTg                               >  1:3226862/1‑68 (MQ=255)
                          gAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTg                               >  1:545282/1‑68 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTAcc                                                  >  1:2749525/1‑45 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGt                           >  1:2832786/1‑68 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTg                          >  1:2722164/1‑69 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTg                          >  1:2552915/1‑69 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTg                          >  1:2826808/1‑69 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTg                          >  1:2099699/1‑69 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTg                          >  1:238584/1‑69 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTg                          >  1:1380209/1‑69 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTg                          >  1:1509914/1‑69 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTg                          >  1:1667127/1‑69 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTg                          >  1:3420831/1‑69 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTg                          >  1:1941365/1‑69 (MQ=255)
                              ccAACGGTAAGTATGTTGACCGTAACGGTAAGGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTg                          >  1:3392516/1‑69 (MQ=255)
                               cAACGGTAAGTATGTTGCCCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGa                         >  1:590031/1‑69 (MQ=255)
                               cAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTTGGGTGa                         >  1:795344/1‑69 (MQ=255)
                               cAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGa                         >  1:3084181/1‑69 (MQ=255)
                               cAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGa                         >  1:683415/1‑69 (MQ=255)
                               cAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGa                         >  1:2764593/1‑69 (MQ=255)
                               cAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGa                         >  1:2626813/1‑69 (MQ=255)
                               cAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGa                         >  1:3103408/1‑69 (MQ=255)
                               cAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGa                         >  1:2433982/1‑69 (MQ=255)
                               cAACGGTAAGTATGTTGACAGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGa                         >  1:2459934/1‑69 (MQ=255)
                                        gTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCa                 <  1:2495244/68‑1 (MQ=255)
                                         tATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACt               <  1:2488687/69‑1 (MQ=255)
                                         tATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACt               <  1:3007124/69‑1 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACGAGGCACTAAc            >  1:2068488/1‑68 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAAc            >  1:3299690/1‑68 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:970041/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:650338/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:946787/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:1706219/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:1948382/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:2238386/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:2251945/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:225977/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:1606042/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:1347860/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:2739865/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:1282278/1‑69 (MQ=255)
                                             ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:2772095/1‑69 (MQ=255)
                                             atGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:3002473/2‑69 (MQ=255)
                                              tGAGCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGACCCAGGCACTaa             >  1:1360756/1‑66 (MQ=255)
                                              tGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGCACCAGGCACTAACg           >  1:2073815/1‑68 (MQ=255)
                                              tGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:1196060/1‑68 (MQ=255)
                                              tGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:1717774/1‑68 (MQ=255)
                                              tGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:1820218/1‑68 (MQ=255)
                                              tGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:1919210/1‑68 (MQ=255)
                                              tGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:2925224/1‑68 (MQ=255)
                                              tGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:2384021/1‑68 (MQ=255)
                                              tGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACg           >  1:260076/1‑68 (MQ=255)
                                                  cGTTACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:1510654/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:3337175/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:2856511/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:1223724/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:1243701/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:1327862/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:2572766/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:1603997/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:3020504/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:242969/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:1847831/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:3145479/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:2050695/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:202821/69‑1 (MQ=255)
                                                  cGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGGGAACCAGGCACTAACGGTCAg      <  1:121297/69‑1 (MQ=255)
                                                   gtaacggtaacgTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAg      <  1:2482252/68‑1 (MQ=255)
                                                      acggtaacgTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGCAcg  >  1:1572724/1‑69 (MQ=255)
                                                             |                                                             
CGGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGCACG  >  minE/2789992‑2790114

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: