Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,597	0	C	T	57.9%	‑3.5 / 10.7	19	R11K (AGA→AAA)	iclR	transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (3/5);  new base T (11/0);  total (14/5)
Fisher's exact test for biased strand distribution p-value = 4.82e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

CGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTG  >  NC_000913/4223452‑4223744                                                                                                                                                  |                                                                                                                                                     cGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTGTTTTGcc                                                                                                                                                  >  1:78667/1‑150 (MQ=255)  gTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGGcg                                                                                                                                                 >  1:146945/1‑150 (MQ=255)     gAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGtt                                                                                                                                             >  1:130212/1‑151 (MQ=255)                         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)                                                                                                                                                  |                                                                                                                                                     CGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTG  >  NC_000913/4223452‑4223744

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