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

TGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCA  >  NC_000913/3404080‑3404350                                                                                                                                            |                                                                                                                                      tgaGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCAACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGtt                                                                                                                                      >  2:201472/1‑140 (MQ=255)     tttGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGgttttc                                                                                                                                  >  2:104164/1‑138 (MQ=255)       tGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTtctc                                                                                                                                >  2:72106/1‑140 (MQ=255)         ccGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTg                                                                                                                              <  1:448367/140‑1 (MQ=255)              tCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGTGGGGAGGGCTGTGGTTTGGGCGGGGGGTTTTTcttgtccac                                                                                                                          >  2:63966/1‑133 (MQ=255)              tCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTATTCTGCCAAc                                                                                                                         >  2:559719/1‑140 (MQ=255)                   tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTTCTCTGCCAACCACgg                                                                                                                    >  2:626784/1‑140 (MQ=255)                   tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTACTCTGCCAACCACgg                                                                                                                    >  1:815708/1‑140 (MQ=255)                       ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                <  1:797379/140‑1 (MQ=255)                        tGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                <  2:804304/139‑1 (MQ=255)                        tGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                >  1:804304/1‑139 (MQ=255)                          gAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGTTTTACTCTGCCAACCACGGTCATGcc                                                                                                             >  1:304808/1‑140 (MQ=255)                            aaaaCAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCATCCACGGTCATGCCAt                                                                                                           <  2:411548/140‑1 (MQ=255)                            aaaaCAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCAt                                                                                                           <  1:104164/140‑1 (MQ=255)                            aaaaCAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCAt                                                                                                           <  2:231081/140‑1 (MQ=255)                             aaaCAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTTCTCTGCCAACCACGGTCATGCCAtt                                                                                                          >  2:797379/1‑140 (MQ=255)                                    tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCTCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTTCTCTTCCCACCACCGGCATGCCCttttttttt                                                                                                   >  2:478534/1‑137 (MQ=255)                                    tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGt                                                                                                                   >  1:204284/1‑124 (MQ=255)                                    tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGt                                                                                                                   <  2:204284/124‑1 (MQ=255)                                            ggggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGt                                                                                           <  1:652439/140‑1 (MQ=255)                                             gggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTTCCCACCACGGGCATGGCACTTATTCTGCGTAAAGT‑c                                                                                         >  2:652439/1‑140 (MQ=255)                                             gggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTTCTTTGCCAACCACGGTCATTCCCTTTTTTTTTCGTAATTT‑c                                                                                         >  2:344875/1‑140 (MQ=255)                                              ggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑cc                                                                                        >  1:827216/1‑140 (MQ=255)                                                  aTTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGGGGCCTGTGGTCT‑GGCGGGGGGGTTTAATCTTCCCACCACCGGCCAGCCCTTTTTTTTTCGGAATGT‑CCGgtt                                                                                    >  2:966435/1‑140 (MQ=255)                                                  aTTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCATTTTTTCTTCGTAATGT‑CCGctt                                                                                    >  1:648719/1‑140 (MQ=255)                                                    tGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGc                                                                                  <  1:881258/140‑1 (MQ=255)                                                    tGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGc                                                                                  <  1:559719/140‑1 (MQ=255)                                                          gagacaggACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAgggg                                                                            <  1:103665/133‑1 (MQ=255)                                                            ttGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGt                                                                           <  2:11811/139‑1 (MQ=255)                                                            ttGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGGTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGt                                                                           >  1:11811/1‑139 (MQ=255)                                                             tGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATG‑CATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGa                                                                        <  1:28486/140‑1 (MQ=255)                                                             tGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGGTTTACTCTGCCAAACCCCGGCATGCCCTTTTTTTTGCGTAATGG‑CCGTTTTCAAGgggggg                                                                         >  1:837264/1‑140 (MQ=255)                                                              ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGt                                                                           <  2:633261/137‑1 (MQ=255)                                                              ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGt                                                                           >  1:633261/1‑137 (MQ=255)                                                              ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGa                                                                        <  2:815708/140‑1 (MQ=255)                                                              ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGa                                                                        <  1:344875/140‑1 (MQ=255)                                                                cGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTTCTCTTCCAAACCCCGTCCCTCCATTTATTTTGCGGAATGTACCTTTTGTAggggggggtt                                                                       >  2:645771/1‑138 (MQ=255)                                                                  aCAAAGTGCTGGCAAAAGTGCTGGCCCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCa                                                                    <  2:154485/140‑1 (MQ=255)                                                                        tGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCCTTTATTTTTCGGAATGT‑CCCGTTGCAAGGGGGGGGTTCAGTaaag                                                              >  2:122277/1‑138 (MQ=255)                                                                         gCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGa                                                             <  2:47692/140‑1 (MQ=255)                                                                         gCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTTCTCTGCCCACCACGGTCATGCCCTTTATTTTTCGTAATGT‑TCGGTTGGCAGGGGGGGGTTTAGTAAAGa                                                             >  2:1001592/1‑140 (MQ=255)                                                                         gCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCCGCGTGGCGGCCTGTGGTCT‑GGCCGGTGGTTTTACCCTGCCGAACACCGGCAAGCCATTTATTCTGGGTAAAGT‑TCGGTTGGCAGGGGGGGGTTTAGGAATGa                                                             >  2:545922/1‑140 (MQ=255)                                                                          cTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGAc                                                            <  1:626784/140‑1 (MQ=255)                                                                          cTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGAc                                                            <  1:645771/140‑1 (MQ=255)                                                                          cTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTTCTCTGCCCACCACGGGCATGCCCTTTATTCTTCGTAATGT‑CCGTTTGCAAGGGGGGGGTTCAGTAATGAc                                                            >  2:103732/1‑140 (MQ=255)                                                                            ggCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGTTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGc                                                          <  1:1001592/140‑1 (MQ=255)                                                                            ggCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGc                                                          <  1:545922/140‑1 (MQ=255)                                                                            ggCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGCGGGCTGTGGTCT‑GGCGGGGGGTTTTTCTCTTCCAACCCCGGTCATGCCATTTTTTCTGCGGAATGT‑CCCGTTTCAAGGGGTGGATTTAGTACTgtttc                                                          >  2:103665/1‑136 (MQ=255)                                                                             gCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGGGGGTTTACTCTGCCAACCACGGGCAAGCCAATTTTTCTGCGTAATGT‑CCGGTTTCAAGGGGGGGGTTTAGTAATGACGcc                                                         >  2:821047/1‑140 (MQ=255)                                                                                                  aaTTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGc                                    <  1:201472/140‑1 (MQ=255)                                                                                                                        ccTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGGGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGAt              >  1:132318/1‑140 (MQ=255)                                                                                                                            tgGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTAc          <  2:181306/140‑1 (MQ=255)                                                                                                                            tgGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTAc          <  1:966435/140‑1 (MQ=255)                                                                                                                                 t‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa     >  1:449084/1‑140 (MQ=255)                                                                                                                                   ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt    <  1:63966/140‑1 (MQ=255)                                                                                                                                   ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt    <  2:304808/140‑1 (MQ=255)                                                                                                                                    gCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATc   >  1:290188/1‑140 (MQ=255)                                                                                                                                     cGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCa  <  1:72106/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                      TGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCT‑GGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGT‑CCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCA  >  NC_000913/3404080‑3404350

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: ‑