Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,429,160	0	C	G	53.3%	1.2 / 17.0	17	L38V (CTG→GTG)	stfR	putative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base G (2/6);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTC  >  NC_000913/1428916‑1429401                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  tATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACCATTCAGCTGAAAGCCAGACTTAACATCACCACGGTGGTGGTGAACACGGTGGnc                                                                                                                                                                                                                                              <  1:181883/250‑1 (MQ=255)                 gTGGCGACTGGCGTACTGACGGAGTCTTCGTTGGGGTCGGTGATAAATTCTTATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTTGGCTTTGGTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCGCAGGAAAACCGGTAC‑AGACCTGCACCATTCAGCGGAAAGCCAGACGTAACAGCGCCACGGTGGTGGTGAACACGGTGGGCTCAg                                                                                                                                                                                                                                          <  1:213237/238‑1 (MQ=255)                          ggCGCACTGACGGTTTCATCGGGGGTGTCGGTTATAAGTTCTGATGAGCCAGGTA‑CACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACAATCCAGCTGAAAGCAAAACTTAACAGCACCACGGTGGGGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGc                                                                                                                                                                                                                     <  1:182338/249‑1 (MQ=255)                                                   ggtCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTTCAGCATg                                                                                                                                                                                                          <  1:185238/235‑1 (MQ=255)                                                               gTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAgn                                                                                                                                                                                                                          <  1:59139/208‑2 (MQ=255)                                                                                      gtgtTATGACAGCCCGCCGGTTCAGGCGGGATTTTTTGGGGGGGGGATATGGCAGTACAGATTTCAGGGGCACTGAAAGAAGGCACAGGAAAACCGGGACCAAAAC‑GCACAATCCAGATGACAGCAAAACCGACCAGCACAACGGTGGTGGTGAACACGATGCCCTCAGAAAATCCGGATGACGCCGGGCGTTACAGCCTGGACGTTGAGCGCGGTAAGTACAGCGCTATTCTGTTGAGGGACGGATTccc                                                                                                                                                        >  1:100541/1‑251 (MQ=255)                                                                                                                 gggCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACCATTCAGCTGAAAGCGAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGTAGGACGGTCAGTACAgtgtcat                                                                                                                                                                             <  1:194167/203‑7 (MQ=255)                                                                                                                    tttttttG‑GGGGTGAATATGGCATGAAAGATTTCAGGTGTAATAAAAGACGGCACAGGAAA‑CCCGTACCAAAACTGCAAAAGCCAGTTGAAAGAAAAAAGCAACAGCACCAGGGTGGTGGTGAACACGCTGGCCGCAGAAAATCCGGATGAAGCCGGGCGTTCCAGCTTGGACGTTGAGTACTGTCAGTACATGTTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTa                                                                                                                          <  1:7013/249‑1 (MQ=255)                                                                                                                                                                    gACGGCACAGGAAAACCGGTAC‑AGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTTAACACGGTGGGTTCAGACAATCCGGATGAAGCCTGGCGTTACAGCATGGa                                                                                                                                                                                                        <  1:84132/125‑1 (MQ=255)                                                                                                                                                                                                  caccattcAGCTGAAAGCCAGACGTCNCAGCACNNCGNTGGTGGGGAACACGGTGGGNTCNNNGNNNCNGGATGNAGNNGGGCGTTACNGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTNCACCATCGCNNGCCGGGACCATCACCGTGTATNAAGNNTNNCNACCNggn                                                                                                       >  2:54309/8‑192 (MQ=255)                                                                                                                                                                                                  caccattcAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGc                                                              <  1:54309/227‑1 (MQ=255)                                                                                                                                                                                                           gCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCGGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTCCGGTCAGTACAGCGTTATTCTGTTGGGGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGGGCCATGACGGAGGGTGATGCCCGTCCGGAGGCACTGCGCAGTTTTg                                    >  1:162786/1‑251 (MQ=255)                                                                                                                                                                                                                aaGCAAAACGTCACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGa                                                                                                                                                                                                                                         >  1:82048/1‑49 (MQ=255)                                                                                                                                                                                                                         gTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGa                       >  1:73356/1‑250 (MQ=255)                                                                                                                                                                                                                                             gtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGAAGTACGGCCACTCCAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCCCGCCGGGACCATCACCGTGACTGAAGATTCACAACTGGGGACGCTGCATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGGGGAAGAGGTGGCGAGTAACGCGTc  >  1:114713/1‑251 (MQ=255)                                                                                                                                                                                                                                               gAACANNATNNNCTNNGAAAATCCGGATGAAGCCGGGCGNTACAGCATGGACGTTGAGTACGGTCAGTANANCNTTATNNTGTNGGTGGAAGGANTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACNCNGTACGCTGAATGATTTTCTCGGTGCCATGACNGAGGATGATGCCCGTCCGGAGGCANTGCGNNggntnnnnnnn                               >  2:166891/1‑210 (MQ=255)                                                                                                                                                                                                                                               gAACACGCTGACCTCAGAAAAGCCGGATGAAGCCGGGCGTTACAGCAGGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTGGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATGACCGTGTATGAAGATTCTCAACCCGGTACGCGGAATGATTTTCTCGTTGCCATGACGGAGGATGATGCCCGGCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAAcgc     <  1:166891/246‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  TATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTC  >  NC_000913/1428916‑1429401

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