Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,378,9680CT47.5% 16.8 / 19.7 40E98K (GAA→AAA) menD2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase; SEPHCHC synthase
Reads supporting (aligned to +/- strand):  ref base C (14/7);  new base T (19/0);  total (33/7)
Fisher's exact test for biased strand distribution p-value = 8.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGT  >  NC_000913/2378821‑2379114
                                                                                                                                                   |                                                                                                                                                   
aaCGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTc‑c                                                                                                                                                 >  1:90306/1‑151 (MQ=255)
      cGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTC‑CGGTTaa                                                                                                                                           <  2:648948/151‑1 (MQ=255)
                     tCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGCTAAGAGAAGCAGTGTTTTTCTCGGTTAACCGGCCttaatta                                                                                                                              >  2:539141/1‑145 (MQ=255)
                      cGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAAACAGTTTTTTTC‑CCGGTAACCCGGCCTccaac                                                                                                                              >  2:494581/1‑144 (MQ=255)
                      cGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTC‑CGGTTAAACCGGCTTCAATCAGt                                                                                                                           >  1:582568/1‑151 (MQ=255)
                             ggCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATACGTTTTTTTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGAt                                                                                                                    >  2:659595/1‑151 (MQ=255)
                                aaTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTTTC‑CGGTTAACCCGGCTTCAATCAGTGGCggg                                                                                                                     >  1:505290/1‑146 (MQ=255)
                                 aTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTTTC‑CGGGTTAACCGGCTTCAATCAGTGCCGGATa                                                                                                                   >  1:221537/1‑148 (MQ=255)
                                        aCTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAACAAAAAAAGTTTTTTTCTGGGTTACACCGGGTT‑AAACAAAGCCGGAGAGAGATTTTcc                                                                                                         >  2:448689/1‑151 (MQ=255)
                                          tgtgCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAAACAGGTTTTTTC‑CGGGTAACCCGGCTTGAATCCGTGCCGGATAAAGAttttg                                                                                                          >  2:597312/1‑146 (MQ=255)
                                           gtgCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCAcc                                                                                                      <  1:279556/151‑1 (MQ=255)
                                               gTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCAccgccg                                                                                                  <  2:545597/151‑1 (MQ=255)
                                                 ggggTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTAAGAGAATCAGTTTTTTTC‑CGGGTAACCCGGCTTCAATCAAGGCCGGATAGAGATTTGCCACCGCCGTg                                                                                                >  1:52133/1‑151 (MQ=255)
                                                     tGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGGTTTTTTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCgg                                                                                            >  1:153532/1‑151 (MQ=255)
                                                            cGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATCAGTTTTTCTC‑CGGGTTACCCGGCTTCAATCAGTGCCGGATAGAGAATTGCCACCCCCCTGCCGGAgg                                                                                         >  1:362423/1‑147 (MQ=255)
                                                            cGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCCGGTTTTTTT‑CGGGTAACCCCGCGTCCATCCGTGCCGGAGAGAGCTTTTTCACCCCCCCGCCGGAGGTcac                                                                                     >  2:232843/1‑151 (MQ=255)
                                                                  ttCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAAAGCATCAGTTTTTCTC‑CCGTTAACCCGGCCTCAATCAGCGACGGATAAAGATTTGCCACCCCCCTGCCCGAGGGCACAATCa                                                                                >  1:461227/1‑150 (MQ=255)
                                                                    cccGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGGTTTTCTC‑CGGTTCAACCCGCTTCAATCAGTGGCGGGGAGAGAATTGCCACCGCCGCGGCGGGGGGCACaa                                                                                   >  2:46062/1‑145 (MQ=255)
                                                                     ccGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTTTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGGTAAAGATTTGCCACCCCCGTGCCGGGGGGCACAAACACCGc                                                                            >  1:286351/1‑151 (MQ=255)
                                                                     ccGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGCTTTTTTC‑CCGGTAACCCGGCCTCAATCACTTCCCGGCAGAGATTTTCCACCGCCGTGCCCGGGGGGCCAATCAcc                                                                              >  2:628237/1‑149 (MQ=255)
                                                                       ggCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTC‑CGGGTAACCCGGGCTCCATCAGTGCCGGGTAGAGAGTTTCCCCCGCCGGGCCGGGGGGCAAAAACACCGcc                                                                           >  2:480027/1‑150 (MQ=255)
                                                                        gCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTTTC‑CGGGTAAACCCGCCTCCATCAGGGGCGGGTAGAGATTTTCCCCCGCCCTGCCGGAGGGCACAATCACACCCAc                                                                         >  2:444149/1‑151 (MQ=255)
                                                                          tGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGGCACAATCACCGCCAcc                                                                        >  1:348811/1‑150 (MQ=255)
                                                                          tGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGGTTTTTTC‑CCGTTAACCCCGCTTCAATCAGTGCCGGATAGAGATTTTCCCCCCCCGTGCCGGGGGGCACAAACACCGCCAc                                                                         >  2:278514/1‑149 (MQ=255)
                                                                             cGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATCAGTTTTTCTC‑CGGTTAACCCCGCCTCAATCAGTGCCGGGTAGAGATTTTCCACCCCCGTGCCGGGGGGCACAATCACCGCCCCCgg                                                                      >  2:324839/1‑149 (MQ=255)
                                                                             cGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTC‑CCGGTTACCCGGCTTCAATCAGTGCCGGATAGAGATTTTCCCCCGCCCTGCCGGGGGGCACaaacccccccccccg                                                                      >  2:212255/1‑136 (MQ=255)
                                                                                 ttGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGTAGCAGGTTTTTTT‑CCGTTTACCCCGCTTCAAACAGTGCCGGATAGAGATTTTCCACCCCCCTGCCGGGGGTCACAATCAcccccccc                                                                        >  2:558891/1‑143 (MQ=255)
                                                                                         ttCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTATGAGAATCAGTTTTTTTC‑CGGTTAAACCCGCTTCCATCAGTGCCGGGTAGAGATTTTCCACCCCCCGGCCGGGGGGGACAAACACCCCCCCCCGCcgctttcttgct                                                         >  1:624178/1‑146 (MQ=255)
                                                                                               ccgcAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTTTC‑CGGGTAACCCGGCTTCAATCAGTGCCGGGTAGAGATTTTCCACCGCCCTGCCGGGGGGCACAATCACCGCCCCCCGCTGCTTGCTGACTTTcccc                                                   >  1:215511/1‑150 (MQ=255)
                                                                                               ccgcAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCa                                                  <  1:271949/151‑1 (MQ=255)
                                                                                                cgcAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTTTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTTCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCCGCTGCTTGCTGACTTTTGcc                                                   >  1:502882/1‑149 (MQ=255)
                                                                                                  cAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAACCATTTTTTCTC‑CGGTTAACCCGGCTTCAATAAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGc                                                <  2:406500/150‑1 (MQ=255)
                                                                                                       aaTTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATAAGTTTTTTTC‑CGGGTAACCCGGCTTCCATTAGTGCCGGATAGAAATTTTCCACCCCCCTGCCGGTGGGCACAAACACCCCCTCCCGCTCCCTGCTGACCTTCCCCacaccccac                                          >  2:93415/1‑143 (MQ=255)
                                                                                                        aTTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTTCCACCCCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCCGCCCCAgc                                          >  1:228084/1‑150 (MQ=255)
                                                                                                                            ggcgtTTAAGAGAATCAGTTTTTCCC‑CGGTGACCCCGGCTTCAATCAGTGCGGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGc                     <  2:461051/151‑1 (MQ=255)
                                                                                                                                   aaGAGAATCAGTTTTTCTC‑CGGTTACCCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATc                <  2:90289/149‑1 (MQ=255)
                                                                                                                                    agagAATCAGTTTTTCTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGa              >  2:558130/1‑150 (MQ=255)
                                                                                                                                        caTCAGTTTTTCTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATgg         >  2:218628/2‑151 (MQ=255)
                                                                                                                                             gTTTTTCTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGtgtg    >  1:104886/1‑151 (MQ=255)
                                                                                                                                               ttttCTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTgtggt  >  1:425847/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCAGTTTTTCTC‑CGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGT  >  NC_000913/2378821‑2379114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑