Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,042,993	0	C	T	42.3%	15.7 / 10.2	26	F193F (TTC→TTT)	srkA	Cpx stress response Thr/Ser protein kinase; MazF antagonist protein
Reads supporting (aligned to +/- strand):  ref base C (9/6);  new base T (11/0);  total (20/6)
Fisher's exact test for biased strand distribution p-value = 2.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

TCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGC‑GATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAT  >  NC_000913/4042854‑4043142                                                                                                                                            |                                                                                                                                                       tCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAGGAATTTACCGTTCTGc                                                                                                                                              >  2:500509/1‑150 (MQ=255)   ggTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGGTCt                                                                                                                                                >  1:584834/1‑146 (MQ=255)         aaCGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCGACATg                                                                                                                                     >  2:871822/1‑151 (MQ=255)         aaCGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGGTCTGCGGCCACAt                                                                                                                                      >  1:321730/1‑150 (MQ=255)              atatCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGAATTCACCGGTCTGCGGCTACAtaaaga                                                                                                                                 >  1:188443/1‑145 (MQ=255)                    cATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGc                                                                                                                             <  1:546151/148‑1 (MQ=255)                      ttGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAAAGTTCACCATTCTGCTGGTCCAAGGATACTGc                                                                                                                             >  2:878103/1‑146 (MQ=255)                           ccACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCggg                                                                                                                    <  2:996110/150‑1 (MQ=255)                              cGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCGGGCACATGGAGAATGCCCCGcccg                                                                                                                     >  2:909075/1‑144 (MQ=255)                               gCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCGGGTACAAGGAGACTGCCACCCCggg                                                                                                                    >  2:716308/1‑146 (MQ=255)                                       tttGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAGTTCACCCTTCTTCCGCTACATGGGGACTGCCCCGCCCGGGATATTTTCTg                                                                                                         >  2:287710/1‑149 (MQ=255)                                          gAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAGTTTAACGCTCTGCGGCTACATGGAGACCGCCCCGCCGGGAATAtttt                                                                                                            >  2:639076/1‑141 (MQ=255)                                            ggACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTGgcg                                                                                                      >  1:978402/1‑147 (MQ=255)                                                       tGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGGCACATGGAGACTGCCACCCCGGGAATATTCTCTGGGGC‑GATGGTCCAAcgtt                                                                                      >  1:436411/1‑151 (MQ=255)                                                                          aGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTGGCGC‑GGTGGTCCCATGTTTGTTGATCTGggtg                                                                        >  1:738386/1‑146 (MQ=255)                                                                          aGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCCCCCTTCTGCGGCTACATGGAGACTGCCCCGCCGGGGATATTCTCTGGCGC‑GGTGGTCCAATGTTTTTTGATCTTgatg                                                                        >  1:161656/1‑146 (MQ=255)                                                                          aGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGC‑GATGGTCCAATGTTTGTTGATCTGGATGATGCa                                                                   <  1:525256/151‑1 (MQ=255)                                                                              gCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTTCCGTTCTTCGGCCACATGGAGACTGCCCCGCCGGGGAAATTCTCTTGCGC‑GGTGGTCCCATGTTTGTTTTTCTGggtg                                                                        >  2:702326/1‑142 (MQ=255)                                                                                       aaaGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGGTTTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTTGCGC‑GGTGGTCCCATGTTTGTTGATCTGGGTGATGCACGTAATGGTcc                                                        >  2:394581/1‑149 (MQ=255)                                                                                        aaGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTTTGCGGCCTCATGGAGACTGCCCCGCCGGGGATATTCTCTTGCGCGGTTGGTCCAA‑GTTTTTTTATCTTGATGATTTACGTAATTGTCCAGc                                                     >  1:1007636/1‑151 (MQ=255)                                                                                        aaGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGC‑GGTGGGCCAAAGTTTGTTGATCTGGGTGATGCACGTAATGGtac                                                        >  1:647754/1‑146 (MQ=255)                                                                                                 aTGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCGGGTACATGGGGACCGCCACGCCCGGGATATTCTCTTGCGC‑GATGGTCCAATGTTTGTTTTTCTGGATGATGCAAGTCATGGGCCAACCGTTCAgg                                             >  2:880685/1‑150 (MQ=255)                                                                                                     gCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTGGCGC‑GATGGTCCCATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGGttt                                         >  1:110511/1‑150 (MQ=255)                                                                                                            gccgccGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGC‑GATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATg                                  <  2:161611/150‑1 (MQ=255)                                                                                                                                  gggAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGC‑GATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGa           <  1:744769/151‑1 (MQ=255)                                                                                                                                          ttCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGC‑GATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCgg    <  2:470335/150‑1 (MQ=255)                                                                                                                                            cACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGC‑GATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAt  <  2:738299/150‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                       TCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGC‑GATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGAT  >  NC_000913/4042854‑4043142

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