Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,404,218	0	T	G	21.3%	105.9 / 10.4	47	G245G (GGT→GGG)	yhdH	acrylyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base T (6/31);  new base G (10/0);  total (16/31)
Fisher's exact test for biased strand distribution p-value = 1.55e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCA  >  NC_000913/3404080‑3404350                                                                                                                                            |                                                                                                                                     tgaGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGtt                                                                                                                                     >  2:248168/1‑140 (MQ=255)    gTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTa                                                                                                                                  <  1:942090/140‑1 (MQ=255)     tttGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTAc                                                                                                                                 <  2:939171/140‑1 (MQ=255)     tttGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGgttttt                                                                                                                                 >  2:243790/1‑138 (MQ=255)      ttGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTAct                                                                                                                                <  2:1027758/140‑1 (MQ=255)       tGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTActc                                                                                                                               >  1:1081880/1‑140 (MQ=255)        gCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTActct                                                                                                                              <  2:367787/140‑1 (MQ=255)         ccGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTg                                                                                                                             <  1:248168/140‑1 (MQ=255)             aTCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTCTTccca                                                                                                                         >  2:751795/1‑138 (MQ=255)              tCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTCTGcccaa                                                                                                                        >  2:816074/1‑137 (MQ=255)               cccGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTACTCTGCCAAcc                                                                                                                       >  1:834432/1‑140 (MQ=255)                ccGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCa                                                                                                                      <  2:834432/140‑1 (MQ=255)                       ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGG‑GGGCGGGTTTTATCTGCCAACCACGGTCAt                                                                                                               >  2:900053/1‑140 (MQ=255)                       ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                               <  2:697265/140‑1 (MQ=255)                       ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                               <  1:751795/140‑1 (MQ=255)                        tGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTCTTCCAACCACGGGCATg                                                                                                              >  2:113535/1‑140 (MQ=255)                        tGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTCTGCCAACCACGGTCATg                                                                                                              >  2:733288/1‑140 (MQ=255)                          gAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGcc                                                                                                            >  1:219360/1‑140 (MQ=255)                          gAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGcc                                                                                                            >  1:184524/1‑140 (MQ=255)                                cAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTAt                                                                                                      <  2:219360/140‑1 (MQ=255)                                cAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTAt                                                                                                      <  2:184524/140‑1 (MQ=255)                                 aaGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTAtt                                                                                                     <  1:70828/140‑1 (MQ=255)                                 aaGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTAtt                                                                                                     <  1:113535/140‑1 (MQ=255)                                   gTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGTGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCt                                                                                                   <  2:205520/140‑1 (MQ=255)                                    tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTg                                                                                                  >  1:19891/1‑140 (MQ=255)                                        ggctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTCTGCCCACCACGGTCCTGCCATTTTTTTTGCGTa                                                                                              >  2:834099/1‑140 (MQ=255)                                             gggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTCTGCCCAACACGGGCATGCCCTTTTTTTTGCGTAAtttt                                                                                         >  2:532649/1‑137 (MQ=255)                                              ggCAATTGATACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTCTGCCCACCACCGGCATGCCCTTTTTTTTGCGTAAAGTcc                                                                                        >  2:277105/1‑140 (MQ=255)                                               gCAATTGACACCGTTGGCGACAAAGTGTTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCg                                                                                       <  1:174276/140‑1 (MQ=255)                                                    tGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGc                                                                                  <  1:786016/140‑1 (MQ=255)                                                      acacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCaa                                                                                <  1:796090/140‑1 (MQ=255)                                                      acacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCaa                                                                                <  1:710320/140‑1 (MQ=255)                                                           gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGt                                                                           <  2:19891/140‑1 (MQ=255)                                                           gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGt                                                                           <  1:277105/140‑1 (MQ=255)                                                           gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCAGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGt                                                                           <  1:759393/140‑1 (MQ=255)                                                            ttGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTCTGGCAACCACCGGCATGCCATTTTTTTTTCGGAAAGTTCGGTTGCAAaggggg                                                                          >  2:630995/1‑138 (MQ=255)                                                            ttGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGGGGGTTTTCTCTGCCAACCACGGGCATGCCATTTTTTTTGCGTAATGTTCGTTTGCAAgggggg                                                                          >  2:769550/1‑138 (MQ=255)                                                             tGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGGGGCCTGTGGTCTGGCGGGGGGGTTTTATTAACCCAACC‑CGGGCATGCCCTTTATTTTGCGTGATGGTCGTTTTGAAGgggggg                                                                         >  2:759393/1‑140 (MQ=255)                                                                 gACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTc                                                                     <  1:274922/140‑1 (MQ=255)                                                                    aaaGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGGGGATTCAGt                                                                  >  1:231299/1‑140 (MQ=255)                                                                         gCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGa                                                             <  1:900053/140‑1 (MQ=255)                                                                          cTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGAc                                                            <  1:327633/140‑1 (MQ=255)                                                                                         gcgcAAATGAATTACGGCGGCTGCGTGGCGGCCTGGGGTCTTGC‑GGGCGGGTTTTCTCTGCCCACCCCCGGCACGGCATTTATTCTTCGGAATGTCCCTTTGCAAGGGGTGGGTTTAGTAAATACCGCACCCAAAcgccc                                             >  2:1181762/1‑139 (MQ=255)                                                                                              aaTGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGcac                                        <  1:769550/140‑1 (MQ=255)                                                                                                            gcTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGa                          <  1:532649/140‑1 (MQ=255)                                                                                                                       gcCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGa               <  1:733288/140‑1 (MQ=255)                                                                                                                            tgGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTAc          <  2:1081880/140‑1 (MQ=255)                                                                                                                                 tGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa     <  1:936910/140‑1 (MQ=255)                                                                                                                                 tGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa     <  1:630995/140‑1 (MQ=255)                                                                                                                                   gC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATc   >  2:6117/1‑140 (MQ=255)                                                                                                                                   gC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATc   <  2:231299/140‑1 (MQ=255)                                                                                                                                   gC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATc   <  1:674015/140‑1 (MQ=255)                                                                                                                                    c‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCa  >  2:506459/1‑140 (MQ=255)                                                                                                                                       ggtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt    >  1:777698/1‑136 (MQ=255)                                                                                                                                       ggtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt    <  2:777698/136‑1 (MQ=255)                                                                                                                                         tggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa     >  1:480300/1‑133 (MQ=255)                                                                                                                                         tggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa     <  2:480300/133‑1 (MQ=255)                                                                                                                                            |                                                                                                                                     TGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGC‑GGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCA  >  NC_000913/3404080‑3404350

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