Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,571,036	0	T	G	45.5%	76.4 / 17.0	55	N239T (AAC→ACC)	eutE	aldehyde oxidoreductase, ethanolamine utilization protein
Reads supporting (aligned to +/- strand):  ref base T (6/24);  new base G (24/1);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 7.93e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTA  >  NC_000913/2570895‑2571182                                                                                                                                               |                                                                                                                                                    aTCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGcc                                                                                                                                                 >  1:969420/1‑146 (MQ=255)       tCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCcgcg                                                                                                                                       >  2:306461/1‑150 (MQ=255)         agcGGCTACGCTACCAACACCAACCGGTACCTTTACGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGGCTGAGCTGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGGGAGCGCCTGCGGc                                                                                                                                     <  2:361420/148‑1 (MQ=255)         aTCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCa                                                                                                                                    <  2:525711/151‑1 (MQ=255)         aTCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCCGCGCCCGCgg                                                                                                                                      >  1:104270/1‑149 (MQ=255)            ggCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATc                                                                                                                                 <  2:448932/151‑1 (MQ=255)                aCGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGa                                                                                                                              <  1:358202/150‑1 (MQ=255)                 cGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGGCCGCGCCCGCGGcca                                                                                                                                   >  2:859787/1‑142 (MQ=255)                 cGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCTGCGGCAATCAGa                                                                                                                              >  1:931773/1‑149 (MQ=255)                  gCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCTGCGGCAATCAGa                                                                                                                              >  1:716498/1‑148 (MQ=255)                     aTCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGtt                                                                                                                          <  1:749994/149‑1 (MQ=255)                     aTCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCCGCCGCAAACAGACGttt                                                                                                                         >  2:556142/1‑150 (MQ=255)                      tCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGCTTCATCCACCACTACCGGC‑GGGGTGCCGGCGCCGGCCGccaataaccagtttt                                                                                                                       >  2:933518/1‑137 (MQ=255)                       caacaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTAt                                                                                                                       <  1:26222/150‑1 (MQ=255)                       caacaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCCGCGCCTGCGGGAATCAGaaggttttt                                                                                                                      >  2:207331/1‑143 (MQ=255)                        aacaacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCCGCGGGCATCACACGtttttt                                                                                                                      >  2:565470/1‑150 (MQ=255)                          caacaaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGCGGGGTTGCC‑CCGCCCCCGGCAAtttgtagttttttag                                                                                                                    >  2:402678/1‑136 (MQ=255)                             caaTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGgtgtg                                                                                                                <  1:672178/151‑1 (MQ=255)                               aTCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGtt                                                                                                              <  2:627487/151‑1 (MQ=255)                                 cAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGGCCGCGCCTGCGGCAAACCGACGGTTATTTGTGTGtt                                                                                                              >  2:536379/1‑149 (MQ=255)                                 cAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCTGCGGGCATCCGACGGTTTTTGGTGGGttt                                                                                                             >  1:991334/1‑150 (MQ=255)                                  aGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTAc                                                                                                           <  2:259637/151‑1 (MQ=255)                                   gTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTAcg                                                                                                          <  2:139102/151‑1 (MQ=255)                                     aCCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTAcgc                                                                                                         <  1:112701/150‑1 (MQ=255)                                            cgtcgGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTc                                                                                                 <  2:626221/151‑1 (MQ=255)                                             gtcgGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCCGCGCCTGCGGCCATCCGACGGTTTTTTGTGTGTTTTCCCGCCGCtt                                                                                                  >  2:94367/1‑149 (MQ=255)                                               cgGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTGcta                                                                                               <  2:428968/150‑1 (MQ=255)                                                 gCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTTCCAGCGCCTGCGGCCATCAGACGTTTTTTTGTGTGTTTTAGcgccgc                                                                                                    >  2:21756/1‑143 (MQ=255)                                                  cacaAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCCGCGGCAATCAGACGTTTTTTGGGGTGTTTACGCGCCCCTTct                                                                                                >  2:814382/1‑146 (MQ=255)                                                  cacaAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCCGCGGGGTTGCCGGCCCCT‑CGGCGATCACGCGGTTCTTTGGGGTTTTACGCCCCCCtttttc                                                                                              >  2:949709/1‑144 (MQ=255)                                                   acaAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCTGCGGCAATCAGAAGGTTATTTGGGTGTTTTCGCGCCGCCTCTACTAcc                                                                                          >  1:617288/1‑151 (MQ=255)                                                   acaAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCCGCGGGAATCAGACGTTTTTTGGGGTGTTTACGCGCCGCCTCTACTAcc                                                                                          >  1:680204/1‑151 (MQ=255)                                                    caAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGGCAGGGGCTGTGGGGATTAGAAGGTTTTTTGTGGGTTTACGCGCCGCTTtttctaac                                                                                          >  2:92395/1‑148 (MQ=255)                                                     aaaTGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGc                                                                                        <  1:515647/151‑1 (MQ=255)                                                              tgttATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCCCTTCTACTACCGCCTCgccgcc                                                                               >  1:86078/1‑151 (MQ=255)                                                              tgttATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCTGCGGCAATCAGACGGTTATTTGTGGGTTTTCGCGCCGCTTCTACTACCGCTTCgccgcc                                                                               >  1:635092/1‑151 (MQ=255)                                                                           agCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACcagc                                                                  <  2:769721/151‑1 (MQ=255)                                                                                 tttGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCTGCGGCAATCAGACGGTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCCTCGCCGCCGGGTACCACCCGCAGGcc                                                             >  1:182652/1‑150 (MQ=255)                                                                                    gACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCACTCACCAATAACGGCGGGGTTCCCGGCG‑CGGCGGCGATTAGGCGGTTATTTGAGGGTTTTCGCGCCCGCTCCACTTCCGCCTCGGCGCCCGGTTCCCCCCGCAgggccgt                                                          >  2:393112/1‑145 (MQ=255)                                                                                     aCGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGGCTGCGGCAATCAGAAGGTTATTTGGGGGTTTACGCGCCGCTTTTACTACCGCTTCGCCGCCGGGTACCCCCAGCCGGGCGGTa                                                         >  1:919240/1‑150 (MQ=255)                                                                                      cGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCCGCGCCTGCGGCAATCAGACGTTTTTTTGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCCGGTACCACCAGCAGGCCGATAcc                                                       >  2:135679/1‑151 (MQ=255)                                                                                        aTGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑g                                                    <  2:960996/151‑1 (MQ=255)                                                                                            aCTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGaaa                                                <  1:31901/151‑1 (MQ=255)                                                                                                    gCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTGAACa                                        <  1:874556/151‑1 (MQ=255)                                                                                                    gCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTTCCAGCGCCTGCGGCAATCAGAAGGTTATTTGGGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTTaa                                          >  1:742626/1‑149 (MQ=255)                                                                                                      aCGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGGTGCCAGCGCCCGCGGCAATCAGACGTTTTTTTGGGGGTTTACGCGCCGCTTCTCATAACGCTTTGCCGCCCGGTATCACCAGCATGGCGGTACCCGGGAAACtga                                            >  2:792890/1‑144 (MQ=255)                                                                                                          gCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTGAACAAGCGtt                                  <  2:104127/151‑1 (MQ=255)                                                                                                              ggTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTGAACAAACGTTgcg                               >  1:887310/1‑150 (MQ=255)                                                                                                                  ggcggTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTGAACAAGCGTTGCGCGGtt                          <  1:508329/151‑1 (MQ=255)                                                                                                                   gcggTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTGAACAAGCGTTGCGCGGttt                         <  2:387183/151‑1 (MQ=255)                                                                                                                             ccaccaCTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCg                <  2:832476/150‑1 (MQ=255)                                                                                                                              caccacTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGa              <  2:440940/151‑1 (MQ=255)                                                                                                                                    tACCGGC‑GGGGTGCCAGCGCCTGCGGCCATCAGACGCTTATTTGTGTGTTTACGCGGCGCTTCTTCTACCGCTTCCCCCCCGGTTACCACCCGCAGGCCGATACCC‑GGAAACTTGAACAAACGTTGCGCGGGTTCGATATCCGGGTTTGcc        >  2:631285/1‑151 (MQ=255)                                                                                                                                         gc‑ggGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGt   <  1:584291/151‑1 (MQ=255)                                                                                                                                          c‑ggGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTa  <  1:449623/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    ATCAGTTCATCGGCTACGCTATCAACAACAATCAGTACCTTTTCGTCGGCACAAATGATGTTGTTATCGAAAGAAGCGCCTTTGACGATGGACTGAGCGGCACGGGCGAGGTCGGCGGTTTCATCCACCACTACCGGC‑GGGTTGCCAGCGCCTGCGGCAATCAGACGTTTATTGGTGTGTTTACGCGCCGCTTCTACTACCGCTTCGCCGCCGGTTACCACCAGCAGGCCGATACCC‑GGAAACTTGAACAAGCGTTGCGCGGTTTCGATATCCGGATTTGCCACAGTA  >  NC_000913/2570895‑2571182

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