Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,997,114	0	T	G	51.6%	31.4 / 40.7	31	D234A (GAC→GCC)	yecC	putative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (11/4);  new base G (16/0);  total (27/4)
Fisher's exact test for biased strand distribution p-value = 4.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-03

TATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAGATCGCCCGGTCCGCAACATCCCGGGCAAAGCTCATTTCGTGCGTCACAATCACCATCGTGCGCTTTTCCTGCGCCAGCTGACGGATGGTGTT  >  NC_000913/1996967‑1997262                                                                                                                                                    |                                                                                                                                                     tattaATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGGGCGGCTTCTCGGGGTCgg                                                                                                                                                   >  1:629783/1‑151 (MQ=255)     aaTGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGGCGGCaa                                                                                                                                                >  1:764707/1‑150 (MQ=255)           ttATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGGGGGGGTGGTCGGGGGCGGCaaaaaaag                                                                                                                                          >  2:541420/1‑148 (MQ=255)              tGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGCGAGGGTGCTCGGGGGCGGCCAATAACGCttt                                                                                                                                      >  1:249061/1‑151 (MQ=255)               gCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGGCGGCAAATAACGCcttt                                                                                                                                     >  2:396733/1‑151 (MQ=255)                     ttGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGGCGGCAAATAACGCTTTTTCgggc                                                                                                                               >  1:494812/1‑149 (MQ=255)                         tAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGCGAGGCTGCTCGGGGGGGGCAAATAAAGCTTTTGCGGccccc                                                                                                                            >  2:443675/1‑150 (MQ=255)                           aGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGGCGGCAAAAAAAGCTTTTTCGGccccc                                                                                                                            >  1:573146/1‑148 (MQ=255)                            gCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCt                                                                                                                        <  2:592498/151‑1 (MQ=255)                                       aGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGGGAGGGTGCTCGGGGTGGGCAAAAAAAGCTTTTGCGGccccc                                                                                                                            >  2:735053/1‑136 (MQ=255)                                           ttGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCcccccc                                                                                                          >  1:532716/1‑150 (MQ=255)                                            tGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGGCGGCAAATAACGCTTTTTCGGCCCCCTGCTCCACTATCCGcccc                                                                                                          >  1:355071/1‑149 (MQ=255)                                              aaaTGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGTATGCGAGGCAACTCGGGGGAGCCGCATAACGCTTTTGCGGGGCCCCTCTCGACTATTcccccc                                                                                                          >  1:70211/1‑147 (MQ=255)                                              aaaTGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGGGAGGGTGCTCGGGGGGGGCAAATAACGCTTTTTCGGCCCCCTGCTCGGCTATCCGCCCCTgg                                                                                                       >  2:53256/1‑150 (MQ=255)                                                  ggCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCa                                                                                                   <  1:475824/150‑1 (MQ=255)                                                          gCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAGATCg                                                                                          <  2:629746/151‑1 (MQ=255)                                                            cccGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGGTGCTCGGGGGCGGCCAATAACGCTTTTTCGGGCCCCCTCTCGAATATTCCCCCCCGGTCCCTAAAAATCGcc                                                                                        >  2:526175/1‑151 (MQ=255)                                                                  aGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGGGAGGCTGCTCGGGGGCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCCCCCCCGGTCCATAAAAATCGCCCGGGcc                                                                                  >  1:566448/1‑151 (MQ=255)                                                                                      ttttttCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTTCCTGGCTGTTCGGGGGCGGGAAAATAAACTTTTGCGGGCCCCCGCTCGGCTATCCGCCCCTGGTCCATAAAAATCGCCCGGTCCGCCACATCCCGGGCaaaac                                                               >  2:162132/1‑148 (MQ=255)                                                                                          ttCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGGCGGCAAATAAAGCTTTTTCGGGCCCCCGCTCGACTATCCCCCCCTCGGCCATAAAGATCGCCCGGGCCGCAACATCCCGGGCACAGCTCAtt                                                          >  2:68560/1‑151 (MQ=255)                                                                                               attattGCAGCAGAAACTTCTCGAGGAACTGGCGGCTGCGAGGCTGCTCGGGGGGGGCACATCACCCTTTTTCTGCCCCCTCCTCGACTATTCCCCCCTGGGCCCTAAAAATCACCCCCGCCACAACAACCCCGGCCAACCTccattct                                                       >  2:425906/1‑143 (MQ=255)                                                                                                 tattGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGCGGCTTCTCGGGGGCGGCAAATAACGCTTTTTCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAAATCGCCCCGTCCGCAACATCCCGGGCCAAACTCCttt                                                         >  2:415958/1‑145 (MQ=255)                                                                                                 tattGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCCCCCCTGGTCCATAAAAATCCCCCGGTCCCCAACATCCCCGGCAAAGCTCATTTCGTg                                                     >  1:350974/1‑149 (MQ=255)                                                                                                 tattGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGCGAGGCTGATCGGGGTCGGCAAAAAACGCTATTGGGGCCCCCTGCCCGACTAACCGCCCCCGGTCCATAAAAATCGCCCGGTCCGCAACATCCCGGGGAAAGCCCATTTTGTGCg                                                   >  1:706515/1‑151 (MQ=255)                                                                                                                           cTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAGATCGCCCGGTCCGCAACATCCCGGGCAAAGCTCATTTCGTGCGTCACAATCACCATCGTGCGCTTTTcc                         >  1:10260/1‑151 (MQ=255)                                                                                                                                 ggTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGG‑CCCCTGCTCGACTATCCGCCCCTGGTCCATAAAGATCGCCCGGTCCGCAACATCCCGGGCAAAGCTCATTTCGTGCGTCACAATCACCATCGTGCGCTTTTCCTGCGcc                   >  1:648966/1‑150 (MQ=255)                                                                                                                                  gTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAGATCGCCCGGTCCGCAACATCCCGGGCAAAGCTCATTTCGTGCGTCACAATCACCATCGTGCGCTTTTCCTGCGcc                   >  2:67674/1‑150 (MQ=255)                                                                                                                                   tGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAGATCGCCCGGTCCGCAACATCCCGGGCAAAGCTCATTTCGTGCGTCACAATCACCATCGTGCGCTTTTCCTGCGcc                   >  2:718174/1‑149 (MQ=255)                                                                                                                                      gAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAGATCGCCCGGTCCGCAACATCCCGGGCAAAGCTCATTTCGTGCGTCACAATCACCATCGTGCGCTTTTCCTGCGCCAGCTg              >  2:651757/1‑151 (MQ=255)                                                                                                                                          ctgctCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAGATCGCCCGGTCCGCAACATCCCGGGCAAAGCTCATTTCGTGCGTCACAATCACCATCGTGCGCTTTTCCTGCGCCAGCTGACgg          >  1:561310/1‑151 (MQ=255)                                                                                                                                                  ggTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAGATCGCCCGGTCCGCAACATCCCGGGCAAAGCTCATTTCGTGCGTCACAATCACCATCGTGCGCTTTTCCTGCGCCAGCTGACGGATGGTGtt  <  1:709723/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAGATCGCCCGGTCCGCAACATCCCGGGCAAAGCTCATTTCGTGCGTCACAATCACCATCGTGCGCTTTTCCTGCGCCAGCTGACGGATGGTGTT  >  NC_000913/1996967‑1997262

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