Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,101,514	0	G	T	23.1%	147.2 / 14.5	52	intergenic (+84/‑176)	sodA/kdgT	superoxide dismutase (Mn)/2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (4/36);  new base T (12/0);  total (16/36)
Fisher's exact test for biased strand distribution p-value = 8.82e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTG  >  NC_000913/4101377‑4101643                                                                                                                                          |                                                                                                                                  tGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTtgtg                                                                                                                                 >  2:495880/1‑140 (MQ=255)        tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAAtt                                                                                                                          >  2:1317922/1‑140 (MQ=255)        tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAAtt                                                                                                                          >  1:1193292/1‑140 (MQ=255)        tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGTCAAtt                                                                                                                          >  2:718135/1‑140 (MQ=255)        tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTTATCAAtt                                                                                                                          >  2:328493/1‑140 (MQ=255)           tgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTCa                                                                                                                       >  1:1366934/1‑140 (MQ=255)            gAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCaa                                                                                                                      <  2:813089/140‑1 (MQ=255)                  ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaaa                                                                                                                <  1:328493/140‑1 (MQ=255)                  ggACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATaaa                                                                                                                >  1:767169/1‑140 (MQ=255)                    aCGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAAATTCAAAATaaaaa                                                                                                              >  2:257199/1‑140 (MQ=255)                       aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGTCAATTTTAAAATAAAAACaa                                                                                                           >  2:1023548/1‑140 (MQ=255)                       aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAATAAAAACaa                                                                                                           >  1:1134930/1‑140 (MQ=255)                        agcagcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAATAAAAACaaa                                                                                                          >  2:1198423/1‑139 (MQ=255)                        agcagcGGCACGTTTTGCGGCGAAAAAATAATAATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATaaaaaaaaa                                                                                                          >  1:1463927/1‑139 (MQ=255)                           agcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCGTCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGAt                                                                                                       <  2:1283962/140‑1 (MQ=255)                            gcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATTAATTTTAAAATAAAAAAAATGATc                                                                                                      >  2:955604/1‑140 (MQ=255)                             cGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGTCAATTTTAAAAAAAAAAAAAAGAAcc                                                                                                     >  2:1074175/1‑140 (MQ=255)                             cGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAAAAAAAAAAATGATcc                                                                                                     >  1:944434/1‑140 (MQ=255)                             cGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGTTCAATTTTAAAAAAAAAAAAATGATcc                                                                                                     >  2:1266474/1‑140 (MQ=255)                                cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                  <  1:483381/140‑1 (MQ=255)                                cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                  <  1:842800/140‑1 (MQ=255)                                cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat                                                                                                  <  2:767169/140‑1 (MQ=255)                                         cggcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAATAAAAACAAAGATCCGctaaaaataaa                                                                                         >  1:575758/1‑140 (MQ=255)                                             gAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAg                                                                                     <  1:750852/140‑1 (MQ=255)                                             gAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAg                                                                                     <  1:1074175/140‑1 (MQ=255)                                             gAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAg                                                                                     <  1:1266474/140‑1 (MQ=255)                                                      atCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAAt                                                                            <  2:533470/140‑1 (MQ=255)                                                        cATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTg                                                                          <  1:1165608/140‑1 (MQ=255)                                                           ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                       <  1:1168974/140‑1 (MQ=255)                                                           ttGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtg                                                                       <  1:317300/140‑1 (MQ=255)                                                              ccgccTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTgg                                                                    <  1:1568464/140‑1 (MQ=255)                                                                    gctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTTAATTGATGTGGTTTTGa                                                              <  2:1242590/140‑1 (MQ=255)                                                                     ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt                                                             <  1:718135/140‑1 (MQ=255)                                                                       gcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCa                                                           <  1:495880/140‑1 (MQ=255)                                                                       gcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCa                                                           <  1:144005/140‑1 (MQ=255)                                                                       gcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCa                                                           <  2:1183263/140‑1 (MQ=255)                                                                             aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa                                                     <  1:955604/140‑1 (MQ=255)                                                                               gCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTAtt                                                   <  2:1366934/140‑1 (MQ=255)                                                                               gCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATAATCCGATAAAAAAAAAAAAACGGTTTAAATTATGTGGGTTTTATCACTtttttt                                                   >  1:763549/1‑140 (MQ=255)                                                                                cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAATAAAAAAAATGATCCGATAAAAAAAAAAAAGCGTTTCAATTGATGTGGTTTTGATCACTttttttt                                                  >  1:997884/1‑139 (MQ=255)                                                                                       gCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATTAATTTCAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGCGGTTTAATTGATGTGGTTTTGATCACTTTTTTTGATTAaat                                           >  1:230904/1‑140 (MQ=255)                                                                                        ccGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatg                                          <  2:182646/140‑1 (MQ=255)                                                                                         cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatga                                         <  2:463945/140‑1 (MQ=255)                                                                                         cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAaatga                                         <  2:770965/140‑1 (MQ=255)                                                                                             tatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGt                                     <  1:315667/140‑1 (MQ=255)                                                                                               tCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCt                                   <  2:645980/140‑1 (MQ=255)                                                                                                     tAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGATAAAAATAAAAAAGCGTTTTAATTTATGTGGTTTTTATCACTTTTTTTTATTAAATTAATGTCTAttttt                             >  2:1514281/1‑139 (MQ=255)                                                                                                        aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                          <  1:1514281/140‑1 (MQ=255)                                                                                                        aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                          <  2:1134930/140‑1 (MQ=255)                                                                                                        aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt                          <  2:763549/140‑1 (MQ=255)                                                                                                           aTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTCTATTGATTAAATGAATGTCTATCTTCGTTTcc                       <  2:594389/140‑1 (MQ=255)                                                                                                                ccATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCaa                  <  2:1193292/140‑1 (MQ=255)                                                                                                                 cATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAac                 <  2:997884/140‑1 (MQ=255)                                                                                                                        aaCGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCAATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATg          <  2:1463927/140‑1 (MQ=255)                                                                                                                             cggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCa     <  1:1317922/140‑1 (MQ=255)                                                                                                                               gtgggTTTTTTTTGATCAATTTTAAAATAAAAACAAAGATCCGATAAAAATAAAAAAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCAtt   >  2:324050/1‑140 (MQ=255)                                                                                                                                tggtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTg  <  2:944434/140‑1 (MQ=255)                                                                                                                                tggtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTg  <  1:18604/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                  TGGAACGTGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTG  >  NC_000913/4101377‑4101643

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