Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,634,676	0	A	G	60.0%	‑6.1 / 32.5	30	D70D (GAT→GAC)	tfaQ	putative tail fiber assembly protein TfaQ
Reads supporting (aligned to +/- strand):  ref base A (8/4);  new base G (7/11);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATTATG  >  NC_000913/1634432‑1634889                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                      cTGCAGCATCCTGAAGCTGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCATGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTc                                                                                                                                                                                                                  >  1:102672/1‑250 (MQ=14) cTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACAGAGTTCAGAATTAAATAACGCGTCGACGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTc                                                                                                                                                                                                                  >  2:197705/1‑250 (MQ=14)          ccTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATg                                                                                                                                                                                                         <  1:224576/250‑1 (MQ=11)                     gcAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtc                                                                                                                                                                                             <  2:3987/251‑1 (MQ=31)                                aCTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGCATC‑AGAACAAGTTTTTCTGATTCCGTCTCATTCACCCAGGATGTGACGTTACACTTCAGAAACTCCACTTACGGCGACAACCAGGCAACATTTTACGGTAACGAACCGAGTTCAGAACTAAACAACGCGTCAAAGGAAGCAACGGCCTAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTg                                                                                                                                                                                   >  1:17951/1‑250 (MQ=2)                                   ggCTAGCTGCATCCGGCTTTTTTATGTTTATTACGCCTCCCGGATCCGGCACA‑GTTTTTCTGCTATCGTATCCTTAACCCAGGCAGTGCCGTTCCCATTCTGATATTCCCCTCCCGGCGATAACCCGGTAAAATTTTCCGGTAACGGACAGGGATCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTATTCAACGAGATGCCACGATGCCTCATAACTGTTGaaa                                                                                                                                                                                >  2:150297/1‑250 (MQ=255)                                             aTCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCAcaa                                                                                                                                                                      >  1:216693/1‑248 (MQ=11)                                                 ggattttttttGTTTCTTCCTCCTACCGGATCCGGTACA‑GTTTTGCTGCTCCCGTATCCTTCCCCCAGGCTGGGCCGTTCCACTTCTGAATCTCCCATTCCGACGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTTCCCCGACGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagcca                                                                                                                                                                 <  2:714/246‑1 (MQ=11)                                                     tttttttGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAagccagc                                                                                                                                                               <  2:192312/248‑1 (MQ=11)                                                          ttGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAat                                                                                                                                                         <  1:70168/250‑1 (MQ=11)                                                                  tcCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCg                                                                                                                                                 <  1:217080/250‑1 (MQ=11)                                                                                      cA‑GTTTATCTGCTTCCGCATCCTTCACGCAGGCTGAGCCGTTCCACTTCTGTTATTCCCCTCCCGGCGATACCCCGGTAACATTTTCCGGTAACGGACCGACTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGCTGCCACGATGCCTCATCACTGTTGAAAACAACCCCAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTg                                                                                                                            >  2:104920/1‑251 (MQ=11)                                                                                                 gCTTCCGTATCCTTCACCCAGGCTGGTCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAAAAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGAGGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATGTCGGTACTGTTTgcaggcag                                                                                                                     <  2:255492/248‑1 (MQ=11)                                                                                                      cGTATCCTTCACCCAGGCTGTGCCGTTCCCCTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCCCGTCATAGACGGTTTTACCCCGATGGTCTTCACCGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACCAAGCCCGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTAt                                                                                                             >  1:74626/1‑251 (MQ=11)                                                                                                           ccTTCCCCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCGGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAgg                                                                                                         <  2:20559/250‑1 (MQ=11)                                                                                                                      gCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGc                                                                                             >  1:10989/1‑251 (MQ=31)                                                                                                                                                    cccGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATAGCGGTACTGTTTGCAGGCAGACCGGTATGATGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTcc                                                               <  1:58149/250‑1 (MQ=11)                                                                                                                                                     ccGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTcc                                                               <  2:28465/250‑1 (MQ=31)                                                                                                                                                       ggCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCCACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGc                                                            >  1:135197/1‑251 (MQ=14)                                                                                                                                                            tAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCcagc                                                        >  2:189541/1‑250 (MQ=31)                                                                                                                                                             aaCCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCcagcag                                                      <  1:138615/251‑1 (MQ=11)                                                                                                                                                                cAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGCAATAAATAACGCGACGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGAtt                                                   >  2:45272/1‑251 (MQ=11)                                                                                                                                                                               ggTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGt                                    >  2:74371/1‑251 (MQ=31)                                                                                                                                                                                         cGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGgttgt                          <  1:157286/251‑1 (MQ=31)                                                                                                                                                                                                   aataaataACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTAATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTc                >  1:262102/1‑251 (MQ=11)                                                                                                                                                                                                          aaCGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGAGGGTCGTCAACGAGATGCCACTATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAAGAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGc         <  1:230922/251‑1 (MQ=11)                                                                                                                                                                                                               gTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGCGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCAttt    >  1:218132/1‑250 (MQ=11)                                                                                                                                                                                                                 cGCCGGAAGCCACGTCATAGACGGTTTTTCCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATTATg  <  1:184081/251‑1 (MQ=11)                                                                                                                                                                                                                 cGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATTATg  <  2:45491/251‑1 (MQ=11)                                                                                                                                                                                                                 cGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATTATg  >  1:8486/1‑251 (MQ=11)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                      CTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATTATG  >  NC_000913/1634432‑1634889

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