Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	807,343	0	C	G	66.7%	26.9 / 21.0	24	intergenic (‑62/+90)	ybhC/ybhB	outer membrane lipoprotein YbhC/putative kinase inhibitor
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base G (9/7);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAA  >  NC_000913/807047‑807636                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                      gcgcAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGTGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTg                                                                                                                                                                                                                                                                                                    >  2:68758/1‑300 (MQ=255) gcgcAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTg                                                                                                                                                                                                                                                                                                    >  1:750863/1‑300 (MQ=255)    caggcaggGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggca                                                                                                                                                                                                                                                                                                 <  1:203728/300‑4 (MQ=255)          gggTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataa                                                                                                                                                                                                                                                                                           >  2:235321/1‑291 (MQ=255)          gggTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGCTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGANATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataa                                                                                                                                                                                                                                                                                           >  2:631735/1‑291 (MQ=255)             tAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGNGTTGAGCTACAGGCGGTCAGCGTCACNCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaa                                                                                                                                                                                                                                                                                        <  2:577470/300‑13 (MQ=255)                   gggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcatt                                                                                                                                                                                                                                                                                  >  2:722206/1‑282 (MQ=255)                   gggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcatt                                                                                                                                                                                                                                                                                  >  1:309916/1‑282 (MQ=255)                     gAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcattgc                                                                                                                                                                                                                                                                                <  1:805328/300‑21 (MQ=255)                       aGGATTCCACGCTGCAGCACCTGGTGTCATGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCNATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTANACTAAGTTggcattataaaaaagcattgctt                                                                                                                                                                                                                                                                              <  2:625160/300‑23 (MQ=255)                          aTTCCACGCTGCAGCACCTGGTGTCAGGGATNCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTNGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattatcaaaaagcattgcttatc                                                                                                                                                                                                                                                                           >  2:615741/1‑275 (MQ=255)                          aTTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAACATTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGAGCGGGCGGGGGGGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcattgcttatc                                                                                                                                                                                                                                                                           <  2:229035/300‑26 (MQ=255)                             ccACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcattgcttatcaat                                                                                                                                                                                                                                                                        <  2:227825/300‑29 (MQ=255)                             ccACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCCAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGNAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGCTCAGCGTCACGCCAAAAGCCAATGCCNGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcattgcttatcaat                                                                                                                                                                                                                                                                        <  2:693548/300‑29 (MQ=255)                               aCGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcattgcttatcaattt                                                                                                                                                                                                                                                                      >  2:497342/1‑270 (MQ=255)                               aCGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcattgcttatcaattt                                                                                                                                                                                                                                                                      >  1:778431/1‑270 (MQ=255)                                                                                                                                                                                                                                                                     tcattttacgtttctcgttcagCTTTTTTATACTAACTTGAGCGAAACGTGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTNAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCTGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTAccc                                <  2:207819/279‑1 (MQ=255)                                                                                                                                                                                                                                                                     tcattttacgtttctcgttcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTAccc                                <  2:504496/279‑1 (MQ=255)                                                                                                                                                                                                                                                                       attttacgtttctcgttcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCgg                              >  1:699385/20‑300 (MQ=255)                                                                                                                                                                                                                                                                                   tcgttcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGt                  <  1:237094/293‑1 (MQ=255)                                                                                                                                                                                                                                                                                            ttttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCTAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTAcgcg         <  1:838441/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                             tttttATACTAACTTGAGCGAAACGGGAAGGNAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCNAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGGACGCGt        >  2:619942/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                   taCTAACTTGAGCGAAACGTGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCTAGGCGCTTGCCAGNGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCaa  <  2:200430/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                   taCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCaa  <  2:371609/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                      GCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAA  >  NC_000913/807047‑807636

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