Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,412,898	0	C	T	42.3%	19.2 / 20.1	26	A190T (GCC→ACC)	acrS	acrAB operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (11/4);  new base T (0/11);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 1.98e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGG  >  NC_000913/3412750‑3413046                                                                                                                                                      |                                                                                                                                                     ccACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCtt                                                                                                                                                    <  1:11743/151‑1 (MQ=255)     tAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGttt                                                                                                                                                >  2:537209/1‑151 (MQ=255)        aaCTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTtat                                                                                                                                              >  2:335587/1‑150 (MQ=255)                      gACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAAccc                                                                                                                                <  2:109938/150‑1 (MQ=255)                       aCGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGc                                                                                                                              <  2:455011/151‑1 (MQ=255)                         gAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTATTTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCa                                                                                                                            <  2:72969/151‑1 (MQ=255)                                       ggAAAGAATGTAAATTTAATTCACACTTAATTCTTTTGTTTGATTAATTAATTTTT‑TTATTTTTTCATCTAGCATAAACATTCTTAATTCTTAATCACCCAGGGGGGGGGTTTTTTTAAAAAAATCATAACCCGCCATATTCATTAACCAg                                                                                                              <  1:502327/151‑1 (MQ=255)                                                         cATGACACTTAATTCATTCGTTTGATGAATTAATTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTg                                                                                            >  2:452447/1‑151 (MQ=255)                                                          atcacatTTATTTCATTCTTTTAAAGAATTATTTTCT‑TTTTTTTTTCATCGGAAATAAACATTTTTAATCCGTTATCGACCAGAGGGGGGTTTTTTTTAAAAAAATCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCt                                                                                             <  1:694137/146‑1 (MQ=255)                                                                   aaTTCTTTCTTTTGATAAATTAATTTCT‑TTTTTTTTTCATCTGGCATAAACATTCTTAATACTTTATCGACCAGGGCGGGGGTTTTTTTAAAAAAATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCa                                                                                   <  2:659382/150‑1 (MQ=255)                                                                              ttaattaaataattTCCTTTATGTTTTAATCTGGCATAAACTTTTTTAATCCTTTTTCGCCCGGGGGGGGGGTTTTTTTAAAAATACCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCaataat                                                                           <  2:199783/139‑1 (MQ=255)                                                                              ttGATGAATTAATTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAatc                                                                       >  2:553460/1‑151 (MQ=255)                                                                                 atgaATTAATTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAatcatc                                                                    >  1:75481/1‑151 (MQ=255)                                                                                 atgaATTAATTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAatcatc                                                                    >  2:395848/1‑151 (MQ=255)                                                                                    atttatttTCT‑TTTTTTTTTCTTTTGGAAAAAACATTTTTAATTCGTTAACGCCCAGGGGGGGGTTTTGTTTAAAAAGATCAAAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAAtcatca                                                                   <  1:425104/147‑1 (MQ=255)                                                                                       taattTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCacaac                                                               >  1:409560/1‑150 (MQ=255)                                                                                                 tttGTTTTCATTTGCCATAAACATTATAAATTCTCTAGCACCCAGAGCGGGGTTTTTTTAAAAAAGTTAACAGCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATc                                                      <  1:265783/148‑1 (MQ=255)                                                                                                   tGTTTTCATCTGTAAAAAACATTCTAAATCCTTTATCGCCCAGGGCGGGGTTTTGTTAAAAAAGATAAAAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCg                                                     <  1:544519/149‑1 (MQ=255)                                                                                                   tGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCg                                                     >  2:689805/1‑149 (MQ=255)                                                                                                     ttttCATCGGGAAAAAACATTTTTAATACTTTATCGACCAGAGGGGGGTTTTTTTAAAAAAGATAAAACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAg                                                   <  1:439367/149‑1 (MQ=255)                                                                                                               acaaaaaaaTTTTTAATACTTTATCGACAAGGGGGGGGTTTTTTTTATAAAAATAATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTac                                       <  2:277553/143‑1 (MQ=255)                                                                                                                     aCATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAAcc                                 >  1:552098/1‑151 (MQ=255)                                                                                                                           ttAATTCTTTATCAACCAGGGCGGGTTTTTTTTAAAAAAAATAAAAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCt                                <  1:276572/146‑1 (MQ=255)                                                                                                                                           cAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAAt           >  2:405738/1‑151 (MQ=255)                                                                                                                                               gggggggTTTTTTAAAAAAAATAATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACtt       <  1:79678/149‑1 (MQ=255)                                                                                                                                                    ggCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCgg  >  2:458298/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCG‑TTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGG  >  NC_000913/3412750‑3413046

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