Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,316,223 G→A 24.0% P34P (CCG→CCA hisF → imidazole glycerol phosphate synthase, catalytic subunit with HisF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,316,2230GA24.0% 157.4 / 40.4 120P34P (CCG→CCAhisFimidazole glycerol phosphate synthase, catalytic subunit with HisF
Reads supporting (aligned to +/- strand):  ref base G (66/25);  new base A (17/12);  total (83/37)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 120 total aligned reads displayed.

TCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGC  >  minE/1316157‑1316289
                                                                  |                                                                  
tCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACt                                                                  <  1:2691097/69‑1 (MQ=255)
             gtggtgAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAacgctac                                                     >  1:2008796/1‑69 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:343066/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:2988166/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:3137102/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:3220290/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:3259459/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:3309695/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:1621760/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:423933/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:2495345/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:2327301/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:484941/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:619549/68‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                <  1:1014874/69‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                <  1:2771786/69‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                <  1:1580153/69‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                <  1:1928300/69‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                <  1:1454834/69‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                <  1:3319325/69‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                <  1:1081371/69‑1 (MQ=255)
                  gAAAGGCGTACAGTTTCGCAACCATGAAATCATAGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:3396918/68‑1 (MQ=255)
                   aaaGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGTCAAAACGCTACGCTg                                                 <  1:463980/67‑1 (MQ=255)
                   aaaGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:2798359/67‑1 (MQ=255)
                   aaaGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:2719347/67‑1 (MQ=255)
                   aaaGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:1240359/67‑1 (MQ=255)
                   aaaGGCGTACAGTATCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:1873204/67‑1 (MQ=255)
                      ggCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCt                                                                  >  1:1475204/1‑47 (MQ=255)
                      ggCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTgaagaa                                            >  1:3161976/1‑69 (MQ=255)
                      ggCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTgaagaa                                            >  1:397026/1‑69 (MQ=255)
                      ggCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTgaagaa                                            >  1:1046165/1‑69 (MQ=255)
                      ggCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTgaaga                                             >  1:1112606/1‑68 (MQ=255)
                      ggCGAACACTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAacgct                                                       >  1:3140386/1‑58 (MQ=255)
                            cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                      >  1:710778/1‑69 (MQ=255)
                            cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                      >  1:1280318/1‑69 (MQ=255)
                            cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                      >  1:3423386/1‑69 (MQ=255)
                            cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                      >  1:1254514/1‑69 (MQ=255)
                            cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                      >  1:335646/1‑69 (MQ=255)
                            cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                      >  1:1636646/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGa                                 >  1:721281/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGa                                 >  1:3196223/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGa                                 >  1:3387736/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGa                                 >  1:2866380/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:2626506/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:150320/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:2574373/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:2560506/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:2558889/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:1945068/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:2365800/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:1804803/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:1432017/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:1953138/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:2235109/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:787311/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:1973648/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:926920/1‑69 (MQ=255)
                                 tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                 >  1:1994539/1‑69 (MQ=255)
                                    cAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                 <  1:2136395/50‑1 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:1451122/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:406950/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:1439842/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:3187635/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:1384436/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:1435297/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:3226056/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:1621837/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:1641473/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:261559/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:1956467/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:2290407/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:1242615/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:915027/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:1915228/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:1851002/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                         >  1:1989036/1‑68 (MQ=255)
                                          tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtg                          >  1:1657747/1‑67 (MQ=255)
                                                cATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGa                  <  1:256602/69‑1 (MQ=255)
                                                          atCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGatat               >  1:1616861/1‑62 (MQ=255)
                                                          atCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGatat               >  1:362123/1‑62 (MQ=255)
                                                          atCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGatat               >  1:1308425/1‑62 (MQ=255)
                                                           tCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCt        >  1:1359087/1‑68 (MQ=255)
                                                           tCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCt        >  1:2275041/1‑68 (MQ=255)
                                                           tCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCt        >  1:3231121/1‑68 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGa                                    >  1:3402483/1‑38 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTGCGATATCACCGCTTcc     >  1:991504/1‑69 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCAcc           >  1:2814960/1‑63 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc     >  1:1338678/1‑69 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc     >  1:54602/1‑69 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc     >  1:563055/1‑69 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc     >  1:1083667/1‑69 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc     >  1:635759/1‑69 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc     >  1:2424740/1‑69 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc     >  1:742862/1‑69 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc     >  1:2182301/1‑69 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc     >  1:2081216/1‑69 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc     >  1:2022199/1‑69 (MQ=255)
                                                             gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTc      >  1:2056421/1‑68 (MQ=255)
                                                                ccGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGc  <  1:3117532/69‑1 (MQ=255)
                                                                ccGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGc  <  1:258959/69‑1 (MQ=255)
                                                                  |                                                                  
TCGTGATGGTCAGGTGGTGAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGC  >  minE/1316157‑1316289

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: