Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,043,030	0	A	G	52.6%	11.5 / 22.6	19	N206D (AAT→GAT)	srkA	Cpx stress response Thr/Ser protein kinase; MazF antagonist protein
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base G (10/0);  total (15/4)
Fisher's exact test for biased strand distribution p-value = 3.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-03

GCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGA  >  NC_000913/4042884‑4043175                                                                                                                                                   |                                                                                                                                                  gCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACAAGGAGACTGCCACGCCGGGGatat                                                                                                                                               >  1:174157/1‑151 (MQ=255)              ggACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTTTTCGGCCACCTGGAGACTGCCACGCCGGGGATATTTTCTTGcgcggt                                                                                                                                  >  2:518603/1‑149 (MQ=255)               gACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTCACCGTTCTGCGGCTACATGGGGACTGCCCCGCCGGGGATATTTTCTTGCGCGATg                                                                                                                                 >  1:133578/1‑151 (MQ=255)                    tACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTTTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTcc                                                                                                                             >  1:86255/1‑150 (MQ=255)                     acacTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCa                                                                                                                            <  1:510054/150‑1 (MQ=255)                                                     ccTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTTTGCGGGTACATGGAGACAGCCCCGCCGGGGATATTCTCTGGCGCGATGGTCCCAAGTTTGTTGAGCTGGGTGATGCACGGaaaag                                                                                            >  1:183967/1‑147 (MQ=255)                                                           aGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATg                                                                                             <  1:19615/143‑1 (MQ=255)                                                           aGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCGGGTACATGGGGACTGCCCCGCCGGGGATATTCTCTGGCGCGGTGGTCCAATGTTTGTTTATTTGGGTGATGCACGTAAttagatcgg                                                                                      >  2:19606/1‑142 (MQ=255)                                                           aGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTCACCCTTCTTCGGCTACATGGGGACTGCCACGCCCGGGATATTCTCTTGCGCGAAGGTCCCATGTTTTTTTGTCTGGATGATGCACCTAATGGTCCCGcc                                                                                     >  2:278026/1‑151 (MQ=255)                                                                cAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAAATTTCACCGTTCTGCGGCTAAATGGAAACTGCCACCCCCGGGATATTCTTTGGGGCGATGGTCCCATGTTTTTTTATCTGGGTGATGCACGTAATGGTCCAGCCGtt                                                                                  >  1:238615/1‑149 (MQ=255)                                                                           aTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTgg                                                                      <  1:106606/150‑1 (MQ=255)                                                                            ttGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATTGAGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGGCCAAAGGTTGTTGGTCTGGGTGAAGCACGTAATTGTCCAGCCGGTCAGGGTTTGTgg                                                                      >  1:634467/1‑149 (MQ=255)                                                                              gccgccGTTACAGCACACTGGCGGGGAGATTTTACCGGTCTGCGGCTACATGGAGACTGCCACGCCGGGGAAATTCTCTGGCGCGGTGGGCCAATGGTTTTTGATCTGGATGATGCACGTAATGGTCCAGCCGGTCAGGGTTtgtg                                                                       >  2:135083/1‑146 (MQ=255)                                                                               ccgccgTTACAGCACACTGGCGGGGAGATTTCACCGGTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCACAGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGGTGtt                                                                 >  1:524714/1‑151 (MQ=255)                                                                                                     ggAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCg                                             >  1:457393/1‑149 (MQ=255)                                                                                                           tttCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCgg                                     >  1:335492/1‑151 (MQ=255)                                                                                                                     cTGCGGCTACATGCAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTCGTTGATCTGGATGATGAACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTgg                           >  2:239796/1‑151 (MQ=255)                                                                                                                               aTGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACtattat                 >  2:163457/1‑151 (MQ=255)                                                                                                                                              ccGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATgaaga  <  2:321787/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGA  >  NC_000913/4042884‑4043175

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: ‑