Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	177,119	0	G	A	45.2%	19.8 / 19.0	31	S169F (TCC→TTC)	yadS	UPF0126 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (11/6);  new base A (0/14);  total (11/20)
Fisher's exact test for biased strand distribution p-value = 1.50e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATGATCCCGCCGCCAACGCCAGTAATGACGCCCA  >  NC_000913/176977‑177260                                                                                                                                               |                                                                                                                                               aaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑tgtg                                                                                                                                       >  1:118443/1‑151 (MQ=255) aaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑tgtg                                                                                                                                       <  2:55192/151‑1 (MQ=255)    tATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑tgtg                                                                                                                                       <  1:534767/148‑1 (MQ=255)                      tctTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTAGCGTGGAc                                                                                                                   >  1:351238/1‑150 (MQ=255)                                      tcATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGataa                                                                                                     >  1:356245/1‑148 (MQ=255)                                                  aaCGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCATCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGGAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGt                                                                                      >  2:115549/1‑151 (MQ=255)                                                     gTCGTTAGCTAAAGATGCCAACGAAAAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTTTTTCCAGTGGTACGAAAAA‑TGTGTAATAAGCCGTGGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTa                                                                                     <  1:544706/149‑1 (MQ=255)                                                        ggTAGTTAAAGATGCCAACGAATAGCCGCCAGCCAAATCAATAGCGTCACGACCATGCCCATCATACTGGTTTTTTCCAGTGGTACGAAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGAttt                                                                                <  2:333090/151‑1 (MQ=255)                                                        ggTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGAttt                                                                                >  2:88008/1‑151 (MQ=255)                                                                aaGATCCCAACAAATAGCCGCCAGCCGAATAAATAGCGTCACGCCCATGCCCATCATACTGGCTGTTTCCAGTGGTACGAAAAA‑TGTGTAATAAGCGGTAGCGGGAAAAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACg                                                                         <  1:567324/150‑1 (MQ=255)                                                                        aaCGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGAAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTaaa                                                                     <  2:67161/146‑1 (MQ=255)                                                                        aaCGAATAGCCGCCAGCCGAATAAATAGCGTCACGCCCATGCCCATCATACTGGTTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTGGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTaaaa                                                                    <  2:498658/147‑1 (MQ=255)                                                                           gAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATggg                                                             >  2:473312/1‑151 (MQ=255)                                                                                 ccgccAGCCGAATCAATAGCGTCACGCCCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑GGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGAttt                                                        <  2:45971/150‑1 (MQ=255)                                                                                 cccccAGCCGAATAAATAGCGTCACGCCCATCCCCATCATACGGGTTTTTTCCAGTGGTACGAAAAA‑TGTGTAATAAGCCGTAGCGGGGAAAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGAt                                                          <  1:555993/145‑1 (MQ=255)                                                                                  cgccAGCCGAATCAATAGCGTCCCGACCATGCCCATCATCCTGGCTTTTTCCAGTGGTACGAAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTcg                                                      <  2:576742/151‑1 (MQ=255)                                                                                  cgccAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTcg                                                      >  2:37711/1‑151 (MQ=255)                                                                                   gccAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGCTTTCTGTACGTAAAATCATGGGGATTTc                                                       >  1:481906/1‑149 (MQ=255)                                                                                        ccGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGcc                                                >  1:380654/1‑151 (MQ=255)                                                                                         cGAATAAATAGTGTCCCGACCATGCCCATCATCTTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAACCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCggg                                                  <  2:46244/148‑1 (MQ=255)                                                                                                       cACGACCATGCCCATCATACGGGTTGTTCCCAGTGGTACGAAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATg                                 <  2:387875/151‑1 (MQ=255)                                                                                                         cgcccATCCCCATCATACGGGCTGTTTCAAGTGGTACGAAAAA‑TGTGAAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATg                                 <  2:176660/146‑1 (MQ=255)                                                                                                          gACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATGATc                              >  1:525410/1‑151 (MQ=255)                                                                                                                      ctgttgggCTGTTCCCATGGGTACGAAAAT‑TGTGTAATAGGCGGTAGCGGGGAAAATCCCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATGATCccg                           <  2:78403/136‑1 (MQ=255)                                                                                                                        tACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTATGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATGATCCCGCCGCCAACGcc                >  1:289432/1‑151 (MQ=255)                                                                                                                             gtctgttCCCGGGGGACGGAAAATTGTTAAATAGCCCGTAGGGGGGAAAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATGATCCCGCCGCCAAc                   <  1:410386/139‑1 (MQ=255)                                                                                                                               tttttcgtttGGTACGAAAAT‑TGTATAATACGCGTTGGCGTGGAAAAGACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCCCGAATGATCCCGCCGCCAACGCCAGTAAt          <  1:148253/141‑1 (MQ=255)                                                                                                                               tGTTTCAAGTGGTACGAAAAA‑TGTGAAATAAGCCGTAGGGGGAACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATGATCCCGCCGCCAACGCCAGTaa           <  2:514539/149‑1 (MQ=255)                                                                                                                                 tttCCAGTGGTACGAAAAA‑TGTGTAATACCCCGTGGCGTGAACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATGATCCCGCCGCCAACGCCAGTAAt          <  1:33718/148‑1 (MQ=255)                                                                                                                                  ttCCAGGGTTACGAAAAA‑TGTGTAAAAGGCCGTAGCGTGAACAATCCCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATGATCCCGCCGCCAACGCCAGTAATGAc       <  2:351303/150‑1 (MQ=255)                                                                                                                                         ggTAGGAAAAA‑TGTGTAATAAGCCGTGGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATGATCCCGCCGCCAACGCCAGTAATGACGCCCa  <  2:243895/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                               AACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAA‑TGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGATTTCTGTACGTAAAATCATGGGGATTTCGCGGGCCAGAACATCACGAATGATCCCGCCGCCAACGCCAGTAATGACGCCCA  >  NC_000913/176977‑177260

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