Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,526	0	T	G	45.5%	36.8 / 14.1	33	T329P (ACC→CCC)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (3/15);  new base G (15/0);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.05e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  NC_000913/366377‑366670                                                                                                                                                        |                                                                                                                                                 aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGggtg                                                                                                                                                 <  1:407161/151‑1 (MQ=255)        ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGttttt                                                                                                                                           <  2:92599/150‑1 (MQ=255)         gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCTGTT‑TTGGGGGCACGGGGGGgttt                                                                                                                                            >  2:586540/1‑149 (MQ=255)                 aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCAc                                                                                                                                   <  1:724165/149‑1 (MQ=255)                  aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                                  >  1:444289/1‑149 (MQ=255)                     agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                                                                                             <  1:173111/151‑1 (MQ=255)                           cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGGGGCCGGGGGGGTTTTT‑TTTTTACCAGTGAGACggg                                                                                                                       >  1:393357/1‑151 (MQ=255)                                   caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTA‑TTGGGGGCCAGGGTGGGTTTTCTTTTCACCAGTGAGAAGGGCAACAGc                                                                                                                >  1:697717/1‑150 (MQ=255)                                      tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGAt                                                                                                            >  1:244142/1‑151 (MQ=255)                                            cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCt                                                                                                      <  2:244133/151‑1 (MQ=255)                                            cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGGGGGTTTTTTTTTTACCCGTGAGAAGGGCAAAAGCTGGTTGccc                                                                                                       >  1:55347/1‑150 (MQ=255)                                                     cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGcc                                                                                              <  2:44832/150‑1 (MQ=255)                                                                 tCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTC‑TTTTGGGGCAGGGGGGGTTTTCTTTTTGCCAGTGGTAGGGGCAACAGGTAATTGCCCTTTCCCGCCCGCCCCtcgg                                                                                    >  2:641472/1‑145 (MQ=255)                                                                  ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTA‑TTGGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCAACAGCCGATTGCCCCTCACCCCCTCGCCCTcagag                                                                                  >  2:228753/1‑149 (MQ=255)                                                                     gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGGGCCCGGGGGGGTTTTCTTTTTACCAATGAGACGGGCAACAGCTGGTTGGCCTTCACCGCCCGGccccgggggaggt                                                                              >  2:522660/1‑140 (MQ=255)                                                                        gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTC‑GTTGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAGAGGGGCGACAGCTGATTGCCCTTTCCCGCCTGGCCCCGCGAGGGtt                                                                              >  2:357335/1‑147 (MQ=255)                                                                        gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTC‑GTA‑TTTGGCGCCCGGGGGGGTTTTCTTTTTACCAGAGAGACGGGCAACAGCTGGTTGGCCCTCACCCCCCGGCCCTGagagag                                                                                >  2:267470/1‑145 (MQ=255)                                                                              ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                    <  2:296875/151‑1 (MQ=255)                                                                                      gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGA‑ATTGGGGGCAGGGGGGGTTTTTTTTTCACCAGGGAGACGGGCAACAACTGATTGCCCTTCACCCCCTGGCCCTGGGAGAGGTGCACCAAGCGGGCCACg                                                             >  1:527089/1‑150 (MQ=255)                                                                                            tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCCTCACCCCCCGGCCCTGAGAGAGGTGCAGCAAGCGGGCCAccccgggtt                                                      >  2:239534/1‑143 (MQ=255)                                                                                             gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GGC‑TTTGGTGCCCGGGGGGTTTTTTTTTTTACCAGTGAGACGGGGCACAGCCGATTTCCCCTCCCCCCCTGGCCCTGCGAGAGTTTCCGCAAGCGGGCCCCGCTGGttt                                                      >  2:293744/1‑150 (MQ=255)                                                                                                ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GGA‑TTTGGCGCCAGGGGGGTTTTTTTTTTTTCCAGGGAGACGGGGCAAAACTGATTGGCCTTCCCCCCCTGGCCCTGAGAGAGTTGCAGCAAACGGGCCACGCCgggt                                                       >  2:480787/1‑144 (MQ=255)                                                                                                     gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                             <  1:370658/151‑1 (MQ=255)                                                                                                          ggCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                                        <  1:418512/151‑1 (MQ=255)                                                                                                               aCGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATc                                   <  1:702289/151‑1 (MQ=255)                                                                                                                     ggggAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGttt                             <  2:144994/151‑1 (MQ=255)                                                                                                                         agagGCGGTTTTC‑GTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCCCGGCGACGGGCGACAGCCGGTTGCCCTTCACCGCCTGGCCCCGAGGGAGTTGGAGCAAGCGGGCCACCCTGGTTTGTCCCCCCCGGCGGAAATCCTTTTTgtt                          >  2:74245/1‑148 (MQ=255)                                                                                                                            ggCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                      <  1:63335/151‑1 (MQ=255)                                                                                                                                   tGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg               >  1:122943/1‑151 (MQ=255)                                                                                                                                     c‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc             >  2:720756/1‑151 (MQ=255)                                                                                                                                             gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat        <  1:401101/150‑1 (MQ=255)                                                                                                                                              ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat      <  1:357348/151‑1 (MQ=255)                                                                                                                                                 gcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa    <  1:35724/150‑1 (MQ=255)                                                                                                                                                  ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa  <  1:33657/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                 AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  NC_000913/366377‑366670

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