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%	15.5 / 24.4	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 = 7.68e-05

AGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAA  >  NC_000913/4042887‑4043164                                                                                                                                                 |                                                                                                                                       aGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATtct                                                                                                                                 <  2:603474/151‑1 (MQ=255)       ttGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTTGGCGGCTACATGGAGGCTGCCCCGCCCGGGATAtt                                                                                                                                   >  1:525253/1‑143 (MQ=255)           ggACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGAATTTAACCCTT‑TGCGGCGTAATGGAGCCTGCCACGGCGGGGATAGTCTCTTGcgcggt                                                                                                                       >  2:69171/1‑149 (MQ=255)               gCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑AATTTACCGTTCTGGGGGTACATGGGGACCGCCACGCCGGGGATATTTTCTTGCGCGGTgg                                                                                                                     >  2:524558/1‑149 (MQ=255)                    acTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAG‑ATTTCACCCTTCTGCGGCTACATGGGGACCGCCCCGCCGGGGATATTCTCTGGCGCGGTGGTccca                                                                                                                >  1:610453/1‑147 (MQ=255)                        aTACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑AATTCACCGGTTTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTTGGGCGGTGGGCCAAtgttt                                                                                                           >  1:5286/1‑150 (MQ=255)                                    ttGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTTTCTTGCGCGATGGGCCCATGTTTGTTTATCTgg                                                                                                 >  1:210964/1‑148 (MQ=255)                                         aGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCa                                                                                         <  2:70138/151‑1 (MQ=255)                                                  ccTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTAACGGTCTTCGGCTACCTGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGTCCCATGTTTGTTTTTTTGGGTGATGCACGTa                                                                                     >  1:683854/1‑146 (MQ=255)                                                   cTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTTTCCGTTTTTGGGCTACATGGGGACTGCCCCGCCGGGGATATTTTCTGGCGCGGTGGGCCAATGTTTGTTGATTTGGATGATGCACGTaa                                                                                    >  1:407426/1‑146 (MQ=255)                                                     gAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGCAG‑ATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGAAATTTTCTGGCGCGGGGGGTCCATGTTTTTTGTTCTGGATGATGCACGTAATTGTcc                                                                              >  2:195916/1‑150 (MQ=255)                                                                                   tACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCCACATAGAGACTGCCACGCCGGGGATATTCTCTGGCGCGGTGGTCCAATGTTTGTTGATCTGGGGGATGCACGTAATGGTCCAGCCGTTCAGGGTTTTTGGGTGTTGCTCa                                                 >  2:210955/1‑149 (MQ=255)                                                                                                                gTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAAc                   <  2:399451/151‑1 (MQ=255)                                                                                                                   cTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGCTGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTgg                >  2:693197/1‑151 (MQ=255)                                                                                                                     gCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGa               >  1:213751/1‑150 (MQ=255)                                                                                                                          tACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACtat         >  1:322626/1‑151 (MQ=255)                                                                                                                             aTGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACtattat      >  2:514888/1‑151 (MQ=255)                                                                                                                             aTGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACtattat      >  1:142289/1‑151 (MQ=255)                                                                                                                                 agaCTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGaa  <  2:99046/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                       AGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAG‑ATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAA  >  NC_000913/4042887‑4043164

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