Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,316,241 T→C 23.1% A40A (GCT→GCC 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,2410TC23.1% 147.3 / 37.0 112A40A (GCT→GCChisFimidazole glycerol phosphate synthase, catalytic subunit with HisF
Reads supporting (aligned to +/- strand):  ref base T (61/25);  new base C (17/9);  total (78/34)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

Only 100 of 120 total aligned reads displayed.

GAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGT  >  minE/1316175‑1316300
                                                                  |                                                           
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:2327301/68‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:3259459/68‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:1621760/68‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:3220290/68‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:3137102/68‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:343066/68‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:484941/68‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:2988166/68‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:735749/68‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:2495345/68‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                           <  1:2771786/69‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                           <  1:3319325/69‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                           <  1:1081371/69‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                           <  1:1928300/69‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                           <  1:1014874/69‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                           <  1:1580153/69‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                           <  1:1454834/69‑1 (MQ=255)
gAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCga                                                           <  1:151073/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)
 aaaGGCGTACAGTATCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:1873204/67‑1 (MQ=255)
    ggCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTgaagaa                                                       >  1:3161976/1‑69 (MQ=255)
    ggCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTgaagaa                                                       >  1:1046165/1‑69 (MQ=255)
    ggCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTgaagaa                                                       >  1:397026/1‑69 (MQ=255)
    ggCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTgaaga                                                        >  1:1112606/1‑68 (MQ=255)
          cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                                 >  1:335646/1‑69 (MQ=255)
          cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                                 >  1:1254514/1‑69 (MQ=255)
          cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                                 >  1:1280318/1‑69 (MQ=255)
          cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                                 >  1:3423386/1‑69 (MQ=255)
          cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                                 >  1:710778/1‑69 (MQ=255)
          cAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCt                                                 >  1:1636646/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGa                                            >  1:721281/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGa                                            >  1:2866380/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGa                                            >  1:3387736/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGatg                                             >  1:3191027/1‑66 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:3397259/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:2626506/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:2574373/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:926920/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:2365800/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:2560506/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:150320/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:1432017/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:1973648/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:1953138/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:1994539/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:1945068/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:1804803/1‑69 (MQ=255)
               tCGCAACCATGAAATCATTGGCGATATCGTGCCACTGGCAAAACGCTACGCCGAAGAAGGTGCAGATGa                                            >  1:2235109/1‑69 (MQ=255)
                  cAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTg                                                            <  1:2136395/50‑1 (MQ=255)
                      cATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgt                                      >  1:1887689/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:3187635/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:1621837/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:1451122/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:1851002/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:1439842/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:1435297/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:3226056/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:1956467/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:1989036/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:1384436/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:406950/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:1242615/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtgt                                    >  1:915027/1‑68 (MQ=255)
                        tGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGgtg                                     >  1:1657747/1‑67 (MQ=255)
                              cATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGa                             <  1:256602/69‑1 (MQ=255)
                                        atCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGatat                          >  1:2747796/1‑62 (MQ=255)
                                        atCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGatat                          >  1:362123/1‑62 (MQ=255)
                                        atCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGatat                          >  1:546969/1‑62 (MQ=255)
                                        atCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGatat                          >  1:1308425/1‑62 (MQ=255)
                                        atCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGatat                          >  1:1616861/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)
                                         tCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCt                   >  1:2432980/1‑68 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTGCGATATCACCGCTTcc                >  1:991504/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:1083667/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:2182301/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:742862/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:2578107/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:2022199/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:2022418/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:1338678/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:54602/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:2081216/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:563055/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:635759/1‑69 (MQ=255)
                                           gTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTcc                >  1:2424740/1‑69 (MQ=255)
                                              ccGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGc             <  1:3276121/69‑1 (MQ=255)
                                              ccGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGc             <  1:3117532/69‑1 (MQ=255)
                                              ccGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGc             <  1:258959/69‑1 (MQ=255)
                                                         acgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCgtgt  <  1:1918070/69‑1 (MQ=255)
                                                         acgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCgtgt  <  1:1319812/69‑1 (MQ=255)
                                                         acgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCgtgt  <  1:2821204/69‑1 (MQ=255)
                                                         acgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCgtgt  <  1:2865931/69‑1 (MQ=255)
                                                         acgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCgtgt  <  1:930903/69‑1 (MQ=255)
                                                         acgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCgtgt  <  1:969262/69‑1 (MQ=255)
                                                         acgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCgtgt  <  1:2951343/69‑1 (MQ=255)
                                                                  |                                                           
GAAAGGCGTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGT  >  minE/1316175‑1316300

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: