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

CGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATG  >  NC_000913/4042844‑4043143                                                                                                                                                        |                                                                                                                                                         cGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCa                                                                                                                                                         <  1:552610/151‑1 (MQ=255)  gCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTTAc                                                                                                                                                        >  1:111751/1‑151 (MQ=255)     cGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTTAcc                                                                                                                                                       >  1:180085/1‑149 (MQ=255)        ccATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGTTCt                                                                                                                                                  <  1:476588/151‑1 (MQ=255)             ggTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GGAGCGTTTGCCGTTCttggg                                                                                                                                              >  1:175020/1‑146 (MQ=255)              gTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GGTGAAGTTTCCCCCGttttacg                                                                                                                                               >  2:584833/1‑144 (MQ=255)               tttGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GGAGATTTTACCGGTCTG‑Cgg                                                                                                                                             >  1:370957/1‑148 (MQ=255)                    aCGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTTACCGTTCTG‑CGGCTAc                                                                                                                                         >  2:140778/1‑147 (MQ=255)                                   aGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GGAGATTTCACCGGTCTG‑CGGGTACATTGGGAC‑TGCCACGCCgg                                                                                                                     >  2:267739/1‑151 (MQ=255)                                               tGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGTTTTG‑CGGCTACATGGGGAC‑AGCCACGCCGGGGATATtct                                                                                                            >  2:304262/1‑148 (MQ=255)                                                gTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTTACCGTTCTT‑CGGGTTCATTGGGAC‑CGCCCCCCCCGGGATAtt                                                                                                              >  1:392409/1‑145 (MQ=255)                                                gTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTg                                                                                                         <  1:139870/150‑1 (MQ=255)                                                  ttGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GGAGAGTTTTCCGTTCTG‑CGGGTACATGGAGAC‑TGCCACGCCGGGGAAATTctct                                                                                                          >  1:160957/1‑147 (MQ=255)                                                  ttGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTTACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTTTCTGgcg                                                                                                      >  2:158531/1‑151 (MQ=255)                                                  ttGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTGgc                                                                                                       <  1:137787/150‑1 (MQ=255)                                                   tGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCCGGGATATTCTCTgg                                                                                                        >  2:290568/1‑148 (MQ=255)                                                          gCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGTTCTC‑CGGCTACATGGAAAC‑AGCCCCGCCGGGAATAATCTCTGGCGCGGTgg                                                                                                >  2:738725/1‑149 (MQ=255)                                                              cacTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGAATTTACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGGATATTCTCTGGCGCGATGGTcc                                                                                             >  2:216481/1‑148 (MQ=255)                                                              cacTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAACGGTTCACCGTTCCC‑CGGCCTCATCGAGAC‑CGCCCCCCCGGGGATATTCTCTTGCGCGTTgggtcca                                                                                           >  1:622374/1‑145 (MQ=255)                                                                 tGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTTACCGGTCTG‑CGGCTACATGGGGAC‑TGCCCCGCCGGGGATATTTTCTGGCGCGATGGTCCCAtgtt                                                                                       >  2:594982/1‑151 (MQ=255)                                                                     aCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTTACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGGATATTTTCTGGCGCGGTGGTCCAAtgtttttt                                                                                   >  1:119293/1‑151 (MQ=255)                                                                       cTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTTAACGGTCTG‑CGGCTACATGGAGAC‑TGCCCCGCCGGGGAAATTTTCTGGCGCGGTGGTCCAATGTttgttg                                                                                  >  1:533139/1‑150 (MQ=255)                                                                               ttGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GGAGAATTTACCGCTCTG‑TTGCCTCATGGAGAC‑TGCCACGCCGGGGATAATCTCTGGCGCGGTGGTCCAATGTTTTTTGAACTgtgt                                                                          >  1:417808/1‑147 (MQ=255)                                                                                tGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGGTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGGATATTCTCTGGCGCGGAGGGCCAATGGTTGTTGGTCTgg                                                                            >  2:646441/1‑147 (MQ=255)                                                                                 gAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTTTCCGATTTGCCGGCT‑CACGGCGGT‑CGTCCCGCCCGGAATATTCTCCAGTGCGACGGTCCCATGTTTTTTGATCTgg                                                                            >  2:7198/1‑146 (MQ=255)                                                                                 gAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGgatgat                                                                       <  1:569695/151‑1 (MQ=255)                                                                                   aggCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGc                                                                     <  1:605237/149‑1 (MQ=255)                                                                                     gCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCAc                                                                   <  2:571164/151‑1 (MQ=255)                                                                                               tGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGG‑GAGGATTTTCCCGTT‑TG‑CCGGTTCATGGGGACAGGCCCCCCCGGGAA‑ATTCTCTGGGGCCCTGGGTCCCTATTTTTTTATTTTGAGCATGCACGTAATTGTcc                                                         >  1:653854/1‑151 (MQ=255)                                                                                                   aGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTTACCGTTCTG‑CGGCTACATGGGGAC‑TGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAATGTTTTTTGATCTGGATGGTGCACGTAATGGGCCa                                                        >  2:71868/1‑148 (MQ=255)                                                                                                   aGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAAAATTTACCCTTCTG‑CGGCTACATGGAGAC‑TGCCACACCGGGGAAATTCTTTGGCGCGATGGTCCAATGTTTGTTGATCTGGGTGAAGCACGTAATGGGCCAGcc                                                     >  2:450120/1‑151 (MQ=255)                                                                                                    gCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GGAGAATTCACCGTTCTT‑CGGCTACAAGGAGAC‑TGCCACGCCGGGAAAATTTTCTGGGGCGATGGTCCAATGTTTTTTGGTCTGGGTGATGCACGGAATGGGCCAGcc                                                     >  1:453535/1‑150 (MQ=255)                                                                                                                                        tGGC‑GG‑GGAGATTTCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATaa                  >  1:422179/1‑149 (MQ=255)                                                                                                                                             gg‑gAAGATTTCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCg             >  2:312304/1‑150 (MQ=255)                                                                                                                                                    aTTTCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCg      <  2:207657/151‑1 (MQ=255)                                                                                                                                                       tCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGa    <  2:111783/150‑1 (MQ=255)                                                                                                                                                        cACCGTTCTG‑CGGCTACTTGGAGAC‑TGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATg  >  2:60030/1‑151 (MQ=255)                                                                                                                                                        cACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATg  >  2:642259/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                         CGCCCGACCATCGGTTTGAACGAATATCTCATTGAGCCACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGC‑GG‑GAAGATTTCACCGTTCTG‑CGGCTACATGGAGAC‑TGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATG  >  NC_000913/4042844‑4043143

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