Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,620	0	T	G	61.1%	1.3 / 11.5	18	A3A (GCA→GCC)	iclR	transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (0/7);  new base G (11/0);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 3.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-04

TTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCAC  >  NC_000913/4223476‑4223758                                                                                                                                                 |                                                                                                                                            ttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGAc                                                                                                                                        <  2:151707/150‑1 (MQ=255)      ccAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAc                                                                                                                                  <  1:309922/151‑1 (MQ=255)         gTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTTCCGCGGTTTGCGGGGAAGGGGGGGGCCCTGACAgg                                                                                                                               >  1:148294/1‑150 (MQ=255)                       acacTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTCGCGGGGATGGGGGGGGCCATGACAGTCTC‑ttttttttttt                                                                                                                >  2:182314/1‑146 (MQ=255)                       acacTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTGTTTTGCCGCGTTTCGCGGGGGTGGGGGGGACCATGACAGGCTC‑TTTTTTTTGtt                                                                                                                >  1:79030/1‑150 (MQ=255)                                  ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGCTTTGCGGGGAAGGGGGGGGGTCTGGCAAGCTCTTTTTTTTTGTATTGTGGaaa                                                                                                       >  2:174303/1‑150 (MQ=255)                                       cGGCAATCCCCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTCGCGGGGATGGGGGGGGCCCTGCCAGGCTC‑TTTTTTCTGTTTCTTGGAAAtcatttt                                                                                                >  2:73246/1‑151 (MQ=255)                                                          tttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGGGGGGCCATGACAGGCTC‑TTTTTTTTGTTTCGTGGGAATCATTTTTATttttttttttttctta                                                                             >  2:276654/1‑145 (MQ=255)                                                                       gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGGGATGGGGGCGGCCATGACAGGCTC‑TTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTTTTTTTGCTAATGCAATAGttt                                                                  >  1:144042/1‑148 (MQ=255)                                                                       gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCCGGTTTGCGGGGATGGGGGCGGCCCTGACCGGCTC‑TTTTTTTTGTTTCGTGGGAATCCTTTTTCTTTTTTTTTTTTGCTTATTCAATAGTTa                                                                  >  1:112669/1‑149 (MQ=255)                                                                           aaGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaac                                                            <  2:25558/151‑1 (MQ=255)                                                                                     cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATg                                                  <  2:130241/151‑1 (MQ=255)                                                                                       gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGGGTTTTCCGGGGAGGGGGGGGGCCCTGACAGGCTC‑TTTTTTTTTTTTTTTGGGAATCATTTTTTTTTTTTTTTTTTTCTTATTCAATAGTTTCTTAACTGATCCggttgg                                                >  1:185978/1‑146 (MQ=255)                                                                                             tCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGGCGCGTTTCGCGGGAATGGGGGCGACCCTGGGCGGCTC‑TTTTTTTTGTATCGTGGGAAACATTTTTATTTTTTTTTTTTGCTAATGCAATTGGTATTGTACTGATTCGATggggttatt                                          >  2:314374/1‑143 (MQ=255)                                                                                                cTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGCGGCCCTGACAGGCTC‑TTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTTTTTTTGCTAATTCAATAGTTACTTAACTTATTCGATGAGTTTATGttt                                       >  1:10943/1‑150 (MQ=255)                                                                                                                       tttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggtg                <  2:400360/151‑1 (MQ=255)                                                                                                                                   ttCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCTTTTTTATTGTTAGCTAATGCAATAGTTTCTGAGCTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc    <  2:101259/151‑1 (MQ=255)                                                                                                                                     cgcgGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcac  <  2:144327/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                            TTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCAC  >  NC_000913/4223476‑4223758

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