Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	581,443	0	C	T	61.1%	0.8 / 13.9	18	intergenic (+64/‑91)	nohD/tfaD	putative DNA‑packaging protein NohD/putative tail fiber assembly protein TfaD
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base T (6/5);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTT  >  NC_000913/581226‑581685                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   gATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTTAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGTCGGGCTTCGCGCACTGTTCAGGCCGGAGCCACAGACCGCCGTCGAATGGGCGGATGCTAATTACTATCTCCCg                                                                                                                                                                                                                   <  1:34612/251‑1 (MQ=255)                                                     gATTTCCTGAAACGGGATATCATCAAAGCCAGGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCggg                                                                                                                                                               <  1:90311/251‑1 (MQ=255)                                                      aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATTAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCGCTCTGTTCAGGCCGGAGCCACAGGCCGCTTTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACTAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGc                                                                                                                                                              <  1:158160/251‑1 (MQ=38)                                                                                                            aaCTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCGCACTGTTCAGGCCGGAGCCACAGACCGCTGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGTGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTc                                                                                                        <  1:214442/251‑1 (MQ=255)                                                                                                            aaCTGATACCGGGGTTACTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATCTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCTGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATg                                                                                                                                                   >  1:119821/1‑208 (MQ=255)                                                                                                                       gggTTGCTGAGTGCATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGAGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCg                                                                                             >  1:59516/1‑251 (MQ=255)                                                                                                                               gAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAAGGCGAGGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTc                                                                                      <  1:107374/250‑1 (MQ=255)                                                                                                                                                     ggTTAACAGGCTGCGGCATTTTGTCGGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTGTGAGCGGGCCATCATGAATGCGATGAGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCGGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATg                                                               <  1:105683/251‑1 (MQ=255)                                                                                                                                                         anNAGGCTGNNGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCAGGAGCCACAGACCGCTGTTGAATGGGCGGATGCTAATTACTATCTCCCg                                                                                                                                                                                                                   <  2:16255/99‑1 (MQ=255)                                                                                                                                                             ggCTGAGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCTGTTGAATGGGAGGATGCTAATTACTATCTCCCGCAAGAATCCGCATACCAGGAAGGGCGGTGGGAAACACTGACCTTTCAGCGGGCCCTCATGAATGCGATGGGCAGCGACTACATCAGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACttt                                                       >  1:115021/1‑251 (MQ=255)                                                                                                                                                              gCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTc                                                                                      >  1:2391/1‑219 (MQ=255)                                                                                                                                                                                      gCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCTGTTGAATGGGCGGATGCTAATTTCTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACCCTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaac                                                                                 >  1:80420/1‑199 (MQ=255)                                                                                                                                                                                      gCTTCGCTCACTGTTCAGGACGGAGCCACAGACCGCTGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACAccc                              >  1:168375/1‑251 (MQ=255)                                                                                                                                                                                                     cAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGcc                                                             <  1:74466/205‑1 (MQ=255)                                                                                                                                                                                                       ggCCGGAGCCACAGACCGCTGTTGACTGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGAGCTGGGAAACACTGCACTTTCAGCGGGCCATCATGAATGCGCTgg                                                                                                                                             >  1:157965/1‑123 (MQ=255)                                                                                                                                                                                                       ggCCGGAGCCACAGACCGCTGTTGAATGNGCGGATGCTANTTANNATCTNCNGNAAGAATCCGCATACCAGGAAGGGCGCTGGNAAACACTGCCCTTTCAGCGGGCCNNCANNNATNNGATgg                                                                                                                                             <  2:157965/123‑1 (MQ=255)                                                                                                                                                                                                              gCCACAGACCGCTGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTAtctggttgccgacggatggtg      >  1:111989/1‑231 (MQ=255)                                                                                                                                                                                                                  caGACCGCTGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTAtctggttgccgacggatggtgatgc  <  1:103180/251‑25 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   GATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTT  >  NC_000913/581226‑581685

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