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

AACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTG‑CCGGAGGTCACAATCACCGCCACCGGCTGCTTGCT‑GACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTC  >  NC_000913/2378821‑2379097
                                                                                                                                                       |                                                                                                                                    
aaCGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTTAGA‑GAAT‑CA‑GTTTTTTT‑cc                                                                                                                                  >  1:225767/1‑151 (MQ=255)
  cGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAGG‑GGAA‑CAGGTTTTTTT‑CCg                                                                                                                                 >  2:485669/1‑151 (MQ=255)
   gTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGt                                                                                                                               <  2:291630/151‑1 (MQ=255)
     gCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTa                                                                                                                             <  1:534680/151‑1 (MQ=255)
       ggggATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GGTTTTTT‑CCgg                                                                                                                                >  2:268206/1‑146 (MQ=255)
                  gggTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CC‑GGTTTTCT‑CCCGTTAACCCCGcct                                                                                                                    >  2:185971/1‑145 (MQ=255)
                               cAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCCGTTAACCCGGCTTCAATCAATGcc                                                                                                         >  1:260685/1‑145 (MQ=255)
                                aaTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTTAGAGGAAT‑CG‑GTTTTTTT‑CCGGTTACCCCGCGTTCa                                                                                                                  >  2:429202/1‑136 (MQ=255)
                                aaTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAGA‑GAAA‑TA‑CGTTTTTT‑CCCGGTAACCCCGCCTCAATTAGGGcc                                                                                                         >  2:722029/1‑144 (MQ=255)
                                 aTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTTAAGG‑GAAT‑CG‑GTTTTTTT‑CTGG‑TAACCCGGGTTTAATCAGGGGCgg                                                                                                       >  2:63665/1‑145 (MQ=255)
                                   gAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTGCCGGATAgaga                                                                                                <  1:465905/150‑1 (MQ=255)
                                            tgCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTTT‑CCGGGTAAACCGGCTTCAATCAGTGGCGGATAGAGATTTTCCAcc                                                                                       >  1:293249/1‑150 (MQ=255)
                                                    gTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTTT‑CCCGGTAAACCGGCTTCAATCAGTGCCGGATAGAGATTTTCCAccccccggcc                                                                               >  1:140871/1‑145 (MQ=255)
                                                            cGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GGAT‑CA‑GGTTTTCT‑CCCGTTTAACCGGCGTCAAACAGGGGCGGGGAGAGATTTGTCCCCCCCGTG‑CCGgggg                                                                         >  2:493131/1‑147 (MQ=255)
                                                             gAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAAA‑GAAT‑CA‑GTTTTTTT‑CCCGTTAATCCGGGCTCCATCAATGCCGGATAGAGATTTGCCCCCCCCCTG‑CCCGGGGTcac                                                                     >  1:812091/1‑150 (MQ=255)
                                                             gAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCCGGTAACCCGGGCTCAAAAAATGCCGGATAGAGATTTGTCACCCCCGCG‑cccgggg                                                                         >  1:251522/1‑141 (MQ=255)
                                                                 aTTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTG‑CCGGAGGTCACAATc                                                                 <  1:451254/150‑1 (MQ=255)
                                                                    cccGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGGGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTG‑CCGGAGGTCACAATCACCg                                                             <  2:689184/151‑1 (MQ=255)
                                                                      cGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAGA‑GGAT‑TA‑GGTTTTCT‑CCGGGTAACCCGGGCTCAATCCAGGGCGGATAGAGATTTTCCACCCCCGGG‑CCGGAGGGCACCAACACcccc                                                           >  2:274821/1‑151 (MQ=255)
                                                                         cTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGGTAACCCGGCCTCAATCAGTGCCGGATAGAGATTTGCCACCCCCCTG‑CCGGAGGGCACAAACACCGCcccc                                                        >  1:414793/1‑151 (MQ=255)
                                                                          tGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAGA‑GAAT‑AA‑GTTTTTCTCCCGTTTAAACCGGC‑TCAATCAGTGGCGGATAGAGAATTTGCGCCCCCCTG‑GCGGAGGAGACAATCAcacaccac                                                        >  1:271346/1‑143 (MQ=255)
                                                                           ggggatttGCTTGTTCCGCGCCGCATTCAATTAGCTCGGCCGGGCAATCGGC‑GGTTAAGA‑GAAT‑AA‑GTTTTTCT‑CCGGTTACCCCGCTTTCAATCAGGGCGGGAAAGAGATGTCCCACCGCCGTG‑CCGGAGGTCACAATCACCGCCACCg                                                       <  2:107930/144‑1 (MQ=255)
                                                                             cGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAA‑CC‑GTTTTTCT‑CCCGTTAACCCGGCTTCCATCAGTGCCGGATAAAGATTTTCCACCGCCCCG‑CCGGAGGTCACAATCACCGCCACCGgc                                                     >  2:75723/1‑150 (MQ=255)
                                                                              gAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGA‑GGAT‑CG‑GGTTTTTT‑CCGGGGAAACCGGGCTTTATCCGTGCCGGACAGAGATTTTGCACCGCCGTG‑CCGGGAGGCACAATCAACGcccaacgc                                                     >  2:166632/1‑143 (MQ=255)
                                                                              gAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGA‑GAAT‑CA‑GTTTTTTT‑CCCGGTAAACCCGCGTCCATCAGGGCCGGATAAAGATTTTCCACCCCCCGG‑CCGGGGGTCACAATCACCCCCCCCGgctg                                                   >  1:925333/1‑151 (MQ=255)
                                                                               aaTTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGA‑GAAT‑CA‑GGTTTTTT‑TCGGGTTACCCGGCCTCAATCAGTGCCGGATAAAGATTTTCCACCCCCGTG‑CCGGGGGGGACCATCACCGCCCCCCgctg                                                   >  2:136473/1‑150 (MQ=255)
                                                                                  tGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTTT‑TCCGTTAAACCCGCTTCAATCAGTGCCGGATAGAGATTTGCCCCCCCCCGG‑CCGGAGGGCACAATCACCGCCCCCCGCtgctt                                                >  2:409808/1‑150 (MQ=255)
                                                                                     cTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑AA‑GTTTTTCT‑CCGGGTTAACCGGCTTCAATCAGTGCCGGGTAAAGAATTTCCCCCGCCGTG‑CCGGAGGTCCCAATCACCGCCCCCGGCtgctt                                                >  1:392683/1‑147 (MQ=255)
                                                                                      tGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGA‑GAAT‑CA‑GTTTTTTT‑CCCGTTAAACCGGCTTCAATCAGTGGCCGAAAAAGATTTTTCACCCCCGTG‑CCGGAGGGCACAATCACCCCCCCCcgccgcct                                                >  1:326302/1‑144 (MQ=255)
                                                                                             cgccgcAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGA‑GAAT‑CC‑GTTTTTTT‑CCCGGTAACCCGGCCTCCATCAGTGCCGGGTAGAGATTTTCCACCGCCGTG‑CCGGGGGTCACAATCAACGCCCCCGGCCGCTTTCC‑GACttt                                      >  2:205290/1‑148 (MQ=255)
                                                                                             cgccgcAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGCTTAGA‑GGATCCA‑GTTTTT‑T‑CCGTGTAACACGGCCTCCATCCGTGACGGCTAGAGATTTGTCACCGCCGTG‑CCCGAGGTCACAAACAACGCCCACCGCtgcttggtcgcttttg                                     >  2:356124/1‑140 (MQ=255)
                                                                                               ccgcAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTACCCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTG‑CCGGAGGTCACAATCACCGCCACCGGCTGCTTGCT‑GACTTTCGCCa                                 <  1:30121/151‑1 (MQ=255)
                                                                                                 gcAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAA‑CA‑GTTTTTTT‑CCGGGTAACCCGGCTTCCATCAGTGCCGGATAGAGATTTGCCACCGCCCTG‑CCGGGGGGCACAATCACCGCCACCCGCCGCTTGCT‑GACTTTTGcc                                  >  1:369918/1‑148 (MQ=255)
                                                                                                    gTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTTAGA‑GGAT‑AC‑GTTTTTTTCCCGGTTACCCCGCCTTC‑ATAAAGGCCGGGTAGAGATTTTCCCCCCCCGTGCCCGGGGGTC‑CAAAAACCCCCACCCGCCGCTTTCT‑GGACTTCGCCAGcccc                            >  2:473279/1‑151 (MQ=255)
                                                                                                       aaTTAGCTCCGGCGGGCGATCGGC‑GGTTAAAA‑GAAT‑CA‑GTTTTTTC‑CCGGCTAACCCCGCTTCAATCAGTGCCGGATAGAGATCTGCCACCCTCCCG‑CCGGAGGGCACAAAAAACGCCAACGGGTTCTTGTT‑GGCTTTCTCCAGccc                             >  1:588770/1‑147 (MQ=255)
                                                                                                       aaTTAGCTCCGGCGGGCGATCGGC‑GGTTAAAA‑GAAT‑CA‑GATTTTTT‑CCCGTTAACCCCGCCTCCATCATTGCCCGATAAAAAATTTCCACCCCCGTG‑CCGGGGGGACCAATCACCCCCACCGGCCTCTTGCTCGACttcgcccccccccgc                          >  2:870304/1‑138 (MQ=255)
                                                                                                       aaTTAGCTCCGGCGGGCGATCGGC‑GGGTAAAA‑GAAT‑CC‑GCTTTTCT‑CCCGGTAAACCGGGTTCAATCAGTGGCGGGTAGAGATTTGCCCCCGCCCTG‑CCGGAGGGCACAATCCCCCCCCCCCGCTGCCTGCC‑GACTTTTGCCCGCCCCCgc                         >  2:908604/1‑151 (MQ=255)
                                                                                                       aaTTAGCCCCGGCGGGCGATCGGC‑GGTAAAGA‑GAAT‑CA‑GTTTTTCT‑CCGCTTACCCCGGCTTCAATCAGGGCCGGATAGAGATTTGCCACCGCCGTG‑CCGGAGGTCACAATCACCGCCACCGGCTGCTTGCT‑GACTTTCGCCAGCCCCAg                          <  2:142720/150‑1 (MQ=255)
                                                                                                        aTTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTTCCACCGCCGTG‑CCGGAGGGCACAATCACCGCCACCGGCTGCTTGCT‑GACTTTTGCCCGCCCCAgc                         >  1:491600/1‑150 (MQ=255)
                                                                                                        aTTAGCTCCGGCGGGCGATCGGC‑GGTTAAAA‑GAAT‑CA‑GTTTTTTT‑CCCGTTAACCCCGCCTCCATCAGTGCCGGATAGAGATTTTCCACCCCCGTG‑CCCGAGGTCACAAACACCGCCACCCGCCGCTTGCT‑GACTTTCCCCAGCCCCCgc                         >  2:924880/1‑150 (MQ=255)
                                                                                                                      gcgATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTTT‑CCCGTTTACCCGGCCTCAATCAGTGCCGGGTAGAGATTTTCCACCGCCCTG‑CCGGGGGGCCCCAACACCGCCACCGGCTGCTTGCT‑GGCTTTTGCCCGGCCCCGCGCCAGATGcccc             >  2:937753/1‑148 (MQ=255)
                                                                                                                        gATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTG‑CCGGAGGTCACAATCACCGCCACCGGCTGCTTGCT‑GACTTTCGCCAGCCCCCGCGCCAGATGCCCCAAccc        >  1:445526/1‑151 (MQ=255)
                                                                                                                         aTCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTG‑CCGGAGGTCACAATCACCGCCACCGGCTGCTTGCT‑GACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCa       <  1:345775/151‑1 (MQ=255)
                                                                                                                           cggc‑ggTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCGCCGCCGTG‑CCGGAGGTCACAATCACCGCCACCGGCTGCTTGCT‑GACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCAc      <  2:421181/150‑1 (MQ=255)
                                                                                                                              c‑ggTTAAGA‑GAAT‑CA‑TTTTTTCT‑CCGGTTACCCCGGCTGCAATCAGTGCCGGATAGAGATTTCCCACCGCCGTG‑CCGGAGGTCACAATCACCGCCACCGGCTGCTTGCT‑GACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTc  <  2:225703/151‑1 (MQ=255)
                                                                                                                              c‑ggTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTG‑CCGGAGGTCACAATCACCGCCACCGGCTGCTTGCT‑GACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTc  >  2:350946/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                    
AACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGA‑GAAT‑CA‑GTTTTTCT‑CCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTG‑CCGGAGGTCACAATCACCGCCACCGGCTGCTTGCT‑GACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTC  >  NC_000913/2378821‑2379097

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