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

TGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATT  >  NC_000913/3404080‑3404352                                                                                                                                           |                                                                                                                                        tgaGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGtt                                                                                                                                        >  1:754088/1‑140 (MQ=255)    gTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTa                                                                                                                                     <  1:432060/140‑1 (MQ=255)      ttGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAc‑t                                                                                                                                  <  2:1070702/140‑1 (MQ=255)         ccGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTg                                                                                                                               <  1:1081679/140‑1 (MQ=255)         ccGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTg                                                                                                                               <  2:29612/140‑1 (MQ=255)          cGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGc                                                                                                                              <  2:171716/140‑1 (MQ=255)          cGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGGTTTAC‑TCTGc                                                                                                                              >  1:55425/1‑140 (MQ=255)             aTCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTTCGTGGCGGCCTGTGGATTGGCGGGATGATTTACTTCTGCCa                                                                                                                            >  2:1198568/1‑140 (MQ=255)               cccGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAAcc                                                                                                                         <  1:244853/140‑1 (MQ=255)               cccGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTAC‑TCTGCCAAcc                                                                                                                         >  1:672219/1‑140 (MQ=255)                ccGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCa                                                                                                                        <  1:1066006/140‑1 (MQ=255)                ccGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCa                                                                                                                        <  1:155287/140‑1 (MQ=255)                   tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACgg                                                                                                                     >  1:1174661/1‑140 (MQ=255)                   tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTC‑TTTGCCGACCACgg                                                                                                                     >  1:118872/1‑140 (MQ=255)                    ccTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCCGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGt                                                                                                                    <  1:611619/140‑1 (MQ=255)                       ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGTATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCAt                                                                                                                 <  2:1146183/140‑1 (MQ=255)                               aCAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTC‑TCTGCCCAACACCGGCAAGCCCtttt                                                                                                         >  2:541629/1‑139 (MQ=255)                                cAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTAC‑TCTTCCAACCACCGGCATGCCCTTTAt                                                                                                        >  2:228205/1‑140 (MQ=255)                                    tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCAt                                                                                                                 <  2:187054/127‑1 (MQ=255)                                    tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCAt                                                                                                                 >  1:187054/1‑127 (MQ=255)                                      tgggctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCg                                                                                                  <  2:55425/140‑1 (MQ=255)                                        ggctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCCGTGGTCTGGCGGGGGGGTTTTC‑TCTTCCAACCACCGGCCTGCCCTTTATTCTGCGTa                                                                                                >  2:687734/1‑140 (MQ=255)                                         gctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTAC‑TCTGCCAACCACGGTCATGCCCTTTATTTTGCGTaa                                                                                               >  2:963362/1‑140 (MQ=255)                                            ggggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGt                                                                                            <  2:1174661/140‑1 (MQ=255)                                             gggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTC‑TCTGCCAACCCCGGGCATGCCCTTTTTTTTTCGTAAtttt                                                                                           >  2:1087332/1‑137 (MQ=255)                                               gCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCg                                                                                         <  2:1156062/140‑1 (MQ=255)                                                    tGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGc                                                                                    <  2:1041047/140‑1 (MQ=255)                                                    tGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGc                                                                                    <  1:666644/140‑1 (MQ=255)                                                    tGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGc                                                                                    <  1:353594/140‑1 (MQ=255)                                                     gACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCa                                                                                   <  2:230178/140‑1 (MQ=255)                                                     gACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCa                                                                                   <  2:479461/140‑1 (MQ=255)                                                           gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGt                                                                             <  2:672219/140‑1 (MQ=255)                                                           gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGt                                                                             <  1:541629/140‑1 (MQ=255)                                                            ttGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTC‑TCTGCCAACCACGGGCATGCCCTTTTTTTTTCGTAATGTCCGTTTTCAAaggggg                                                                            >  2:1169023/1‑138 (MQ=255)                                                            ttGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTC‑TCTGCCAACCACGGGCATGCCATTTTTTTTGCGTAATGTCCGTTTTCAAgggggg                                                                            >  2:1158205/1‑138 (MQ=255)                                                              ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGa                                                                          <  1:16309/140‑1 (MQ=255)                                                                cGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGAtt                                                                        <  2:754088/140‑1 (MQ=255)                                                                     aaGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCAATTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTa                                                                   <  1:767823/140‑1 (MQ=255)                                                                       gTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAAt                                                                 <  1:1158205/140‑1 (MQ=255)                                                                                   gTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGa                                                     <  1:228205/140‑1 (MQ=255)                                                                                              aaTGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGcac                                          <  1:687734/140‑1 (MQ=255)                                                                                                         gcggcTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAg                               <  2:433620/140‑1 (MQ=255)                                                                                                            gcTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTC‑TCTTCCCACCACGGTCATGCCCTTTTTTCTTCGTAAAGTCCGGTTGCAAGGGGGGGGTTCAGTAATGACGCCCCCAGAACGCCGCGCACAAACCCGGCcgcgg                            >  2:637747/1‑139 (MQ=255)                                                                                                              tGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTC‑TCTGGCAACCCCCGGCATGCCATTTTTTCTGCGTAATGTCCGTTTGCAAGGGGGGGGTTCAGTAAAGACGCCCCCAGAACGCCGCGCACAAACCCGGCAGCGACt                          >  2:503608/1‑140 (MQ=255)                                                                                                                            tgGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTAc            <  2:118872/140‑1 (MQ=255)                                                                                                                               tCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCgg         <  1:199019/140‑1 (MQ=255)                                                                                                                                 tGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa       <  1:1087332/140‑1 (MQ=255)                                                                                                                                 tGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa       >  2:914985/1‑140 (MQ=255)                                                                                                                                  ggCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt      <  2:1187891/140‑1 (MQ=255)                                                                                                                                  ggCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt      <  2:1156547/140‑1 (MQ=255)                                                                                                                                   gCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATc     <  1:1169023/140‑1 (MQ=255)                                                                                                                                   gCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATc     >  2:604788/1‑140 (MQ=255)                                                                                                                                    cGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCg          >  2:1016839/1‑134 (MQ=255)                                                                                                                                    cGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCg          <  1:1016839/134‑1 (MQ=255)                                                                                                                                      ggtggtTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCAtt  >  1:733170/1‑140 (MQ=255)                                                                                                                                      ggtggtTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCAtt  >  2:958745/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                        TGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAC‑TCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATT  >  NC_000913/3404080‑3404352

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