Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,042,993	0	C	T	57.7%	2.4 / 22.9	26	F193F (TTC→TTT)	srkA	Cpx stress response Thr/Ser protein kinase; MazF antagonist protein
Reads supporting (aligned to +/- strand):  ref base C (10/1);  new base T (15/0);  total (25/1)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-03

CCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGC  >  NC_000913/4042846‑4043135                                                                                                                                                    |                                                                                                                                               cccGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGGTTTTAcc                                                                                                                                             >  2:525558/1‑151 (MQ=255)              ttGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCgg                                                                                                                                    >  1:480620/1‑147 (MQ=255)                       tatCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTAACGTTCTTCGGCTACATGGAGACt                                                                                                                       >  1:615769/1‑151 (MQ=255)                                 aGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTTCGGCTACATGGAGACTGCCCCCCCgg                                                                                                             >  1:673750/1‑151 (MQ=255)                                        cAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATTGAGACTGCCACGCCGGGGATAtt                                                                                                      >  2:792768/1‑151 (MQ=255)                                         aaGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACAGTTCTGCGGCTACATGGGGACTGCCACGCCGGGAGTAtt                                                                                                      >  1:798630/1‑150 (MQ=255)                                         aaGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCCCGCCGGGGATAtt                                                                                                      >  2:15941/1‑150 (MQ=255)                                                 tGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTGgcg                                                                                              >  2:353524/1‑150 (MQ=255)                                                           acacTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTTTGCGGCTACATGGAGACTGCCACCCCCGGGATATTCTCTTGCGCGGTGGTcc                                                                                     >  1:615644/1‑149 (MQ=255)                                                              cTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTCACCCTTCTTCGGCTACATGGGGACCGCCACGCCGGGGATATTTTCTGGCGCGATGGGccca                                                                                   >  1:49084/1‑146 (MQ=255)                                                                   aCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTTTGCGGCTAACTGGGGATTGCCCCACCGGGGATATTTTTTTGCGGGGAGGGCCAAtgtttttt                                                                           >  2:674285/1‑151 (MQ=255)                                                                        ccGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTTTGCGGCTACATGGAGACCGCCCCGCCGGGGATATTTTCTGGCGCGATGGTCCCAtgtttgtt                                                                           >  2:250287/1‑146 (MQ=255)                                                                             ttGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTAAATGGAGACTGCCCCGCCGGGAATATTTTCTGGCGCGGTGGTCCAATGTTTGTTGATCTGggtg                                                                 >  2:377262/1‑151 (MQ=255)                                                                              tGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTCTGCGGGGACATGGGGACTGGCACCCCGGGGATATTTTCTGGCGCGGTGGGCCAAtgtttttttgttt                                                                      >  1:275498/1‑140 (MQ=255)                                                                                aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGTAGAATTTAACGTTCTGCGGCTCCATGGAGAGTGCCCCGCCGGGAATATTTCTTGGGGCGATGGTCCAATATTTGTTGATTTGggtg                                                                 >  1:278069/1‑148 (MQ=255)                                                                                  aGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGGCCAATGTTTGTTTATCt                                                                      >  2:774795/1‑141 (MQ=255)                                                                                             tGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTCTGCGGCTACATGGGGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGTCCCATGTTTGTTGATCTGGATGATGCACGGAATGGTcc                                                 >  2:243069/1‑151 (MQ=255)                                                                                                aaGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCa                                                >  2:788359/1‑149 (MQ=255)                                                                                                                    gccgccGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATg                           <  2:270923/150‑1 (MQ=255)                                                                                                                    gccgccGTTACAGCACACTGGCGGGAAGATTTCACCCGTCTGCGGCTACATGGGGACAGCCACCCCGGGGATATTTTCTGGCGCGATGGTCCAATGTTTGTTGTTCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGGTTTTTgg                              >  2:596733/1‑147 (MQ=255)                                                                                                                       gccgTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCCACATTGAGACAGCCCCGCCGGGAATATTTTCTCGCGCGGTGGGCCCATGTTTGTTGGTCTGGGTGATGGACGGAATGGTCCAGCCGGTCAGGATTTGTgg                              >  1:125851/1‑144 (MQ=255)                                                                                                                                gCACACTGGCGGGAAGATTTTACCGTCCTGCGGCTACCTGGGGGATGCCACCCCCGGGATATTCTCTGGCGCGATGGGCCACAGTTTTTTGATCTGGGTGGTGGACGGGATGGGCCAGCCCTTTCGGGTTTGGGGGTGTTTCTCAATgg                >  1:180909/1‑149 (MQ=255)                                                                                                                                     cTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATaa          >  2:271023/1‑150 (MQ=255)                                                                                                                                         cGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCg     >  1:803108/1‑151 (MQ=255)                                                                                                                                           ggAAGATCTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCg     >  1:532289/1‑149 (MQ=255)                                                                                                                                             aaGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagc  >  2:393486/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                               CCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGC  >  NC_000913/4042846‑4043135

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