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

ACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCG  >  NC_000913/4042863‑4043132                                                                                                                                   |                                                                                                                                            aCGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTCCCCCTTTTTCGGccaacagg                                                                                                                         >  1:682492/1‑144 (MQ=255)       atCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTCACCGTTCTGCGGCTACATGGAGACt                                                                                                                    >  1:612677/1‑150 (MQ=255)         ctcATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTCTGGGGGCACATGGAGACCGcc                                                                                                                 >  2:113501/1‑151 (MQ=255)                   cACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGCAGATTTTACCGGTCTTTGGCTACATAGTAATTGCCCCTCCgggg                                                                                                        >  1:631594/1‑149 (MQ=255)                   cACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGaa                                                                                                       <  2:870836/151‑1 (MQ=255)                     cGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCCCCATGGAGACTGCCACCCCggg                                                                                                         >  2:263328/1‑147 (MQ=255)                                            acTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGGCTGCCACCCCGGGGATATTTTCTGGCCCGGTGGGCCAAt                                                                               >  2:58874/1‑150 (MQ=255)                                                aTACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGGGGCGACATTGAGACTGCCCCGCCGGGGATATTCTCTGgcgggcggggcctcagt                                                                             >  2:328375/1‑134 (MQ=255)                                                  aCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGGGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAAtgtttgtt                                                                        >  1:297669/1‑151 (MQ=255)                                                      tCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCGGCTACATGGAGGCTGCCACGCCCGGGAGATTCTCTGGGGCGCTGGGCCAAtggttgtt                                                                        >  2:302486/1‑147 (MQ=255)                                                       ccGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTTTCTGGCGCGATGGTCCAATGTTTGTTGATCt                                                                   >  1:52964/1‑151 (MQ=255)                                                               aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCCTTCTGCGGCTACATGGGGACTGCCACGCCGGGGATATTTTTTGGCGCGATGGGCCAATGTTTGTTGATCTGgatg                                                              >  1:751787/1‑148 (MQ=255)                                                               aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTTTTCGGGCACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGGTGGTCCCATGTTTTTTGATCTTggtgat                                                            >  2:15227/1‑150 (MQ=255)                                                               aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCCGTCTGCGGCTACATGGAGACTTCCCCGCCGGGGATATTTTCTTGCGCGATGGTCCCATTTTTTTTGAt                                                                     >  2:477630/1‑141 (MQ=255)                                                                          ccTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCCACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGt                                                >  1:814462/1‑151 (MQ=255)                                                                           cTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTAACGTTCTGCGGCTACCTGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCCCGTAAAgg                                                 >  2:661369/1‑149 (MQ=255)                                                                                aGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTCACCGTTCTTCGGCTACATGGAGACTGCCCCGCCGGGGATATTTTCTGGCGCGATGGTCCAATGTTTGTTTGTCTGGATGATGCACGGAATGGTCCAGcc                                          >  2:210599/1‑151 (MQ=255)                                                                                 gCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATTGTCCAGcc                                          >  1:333466/1‑150 (MQ=255)                                                                                     cAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCCCGCCGGGGATATTTTCTGGCGCGGTGGGCCCATGTTTTTTGTTCTGGATGATGCACGTAATGGTCCAGCCGTTc                                      >  2:473678/1‑150 (MQ=255)                                                                                        aTGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTTTGCGGCTACATGGAGACTGCCACGCCGGGAATATTTTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGGTGATGGACGTAATGGTCCAGCCGTTCAgg                                   >  1:855189/1‑150 (MQ=255)                                                                                        aTGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGGGACTGCCACGCCGGGGATATTTTCTGGCGCGGTGGTCCCATGTTTGTTGATCTGGGTGATGCCCGTAATGGTCCAGCCGTTCAgg                                   >  1:25148/1‑150 (MQ=255)                                                                                              tGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCGTTCTGCGGCTACATGGAGAATGCCCCGCCGGGGATATTCTCTTGCGCGGGGGGCCCATGTTTGTTGTTCTGGGTGATGCACGTAATTGTCCAGCCGGTCAGGGTTtgt                             >  2:772573/1‑150 (MQ=255)                                                                                                aTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTTGCGCGATGGGCCAAAGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGGTTTTTgg                           >  1:180183/1‑150 (MQ=255)                                                                                                 ttGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCCTTCTGCGGCCACATGGAGACTGCCCCGCCGGGGATATTTTCTGGCGCGATGGGCCAATGTTTGTTGGTCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTgg                           >  2:115857/1‑149 (MQ=255)                                                                                                        cgTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCCGCTACATGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGTCCCATGTTTGTTGATCTGGGTGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTc                  >  1:590342/1‑151 (MQ=255)                                                                                                                         gggAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCg  >  2:812278/1‑150 (MQ=255)                                                                                                                                   |                                                                                                                                            ACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCG  >  NC_000913/4042863‑4043132

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