Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	578,427	0	G	C	50.0%	‑6.3 / 20.9	18	T107T (ACG→ACC)  R34P (CGG→CCG)	rzpD rzoD	putative prophage endopeptidase RzpD putative prophage lysis lipoprotein RzoD
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base C (5/4);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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.

TAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGC  >  NC_000913/578243‑578663                                                                                                                                                                                         |                                                                                                                                                                                                                                             tatcGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTCCGCTCGATGCCCACAACACGACGGAGTTAGCTGCTGCGAAAACTGAAACTGATGCTCTGCGTGATGCTGTTGCCGCTGGTCGTCGTCGGCTGCACATCCAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATa                                                                                                                                                                                                                              >  1:84188/4‑201 (MQ=11)            ncAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAATTACACGCAGGAGTTAGCTGATGCGAAAGCTGAAACTGCTGCTCTGCGTGCTGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTTGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTg                                                                                                                                                                                             >  1:20963/2‑223 (MQ=11)                   cTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATc                                                                                                                                                          >  1:36587/1‑251 (MQ=40)                                 cttCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTCAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCCCCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGGAGACCCCGCTGAACGGGCTTATTTCACCCTCAGAGAGAGGCTGATCagatcgggaaagcc                                                                                                                                            >  1:157952/3‑238 (MQ=11)                                                gctgcGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAAccn                                                                                                                              <  1:103074/250‑2 (MQ=255)                                                          aTGCAAAATACACGAAGGAGTTAGCGGATGCTAAAGCTAAAAAGGATGCTCTGCGTGATGATTTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaac                                                                                                                                          <  1:9582/228‑1 (MQ=14)                                                                                    naTGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATagagttgccc                                                                                         >  1:44391/2‑245 (MQ=255)                                                                                                gAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCAACGTCGGTTGCACATCAAAGCAGCCCGTCAGCCAGTGCGTGCAGCAAACACCGCCTCCGGCGTAGATACTGCAGCCTACCCCCCACTGGCAAACACCGCTGAACGCGCTCATTTCACCCTCCGAGAGAGGCTGATCACTATGCAAAAACAACTGGAACGAACCCAGAAGTATATTAATGAGCAATACAGATAGAGTTGCCCATATCGATGGGc                                                                             >  1:129923/1‑251 (MQ=18)                                                                                                                 tgatgaTGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACGGCTGAACGGGATTATTTCACCCGCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGTAACCCAGAAGTATATTAATGAGCAg                                                                                                          <  1:53881/205‑1 (MQ=255)                                                                                                                  gatgatGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAattattg                                                            <  1:146204/250‑2 (MQ=255)                                                                                                                  gatgatGTTGCCGCTGGTCGTAGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACACTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCCGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCCTATCGATGGGCAACTCATGCAattattgt                                                           >  1:173380/1‑249 (MQ=21)                                                                                                                                            ctnnnnntcanaaGCGGTCGGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaaca                                                                                                                                            <  2:104102/133‑1 (MQ=17)                                                                                                                                                      naaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaaca                                                                                                                                            >  1:104102/2‑134 (MQ=255)                                                                                                                                                        aGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGATCAATACACACGCGCTTCCAGCGGAGTATAAATGcc                     <  1:46503/251‑1 (MQ=255)                                                                                                                                                                  cagtcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGCTTATTTCACCCTCAGAGAGAgg                                                                                                                                                                >  1:27384/1‑102 (MQ=255)                                                                                                                                                                    gtcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATAGTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATg                       <  1:51073/237‑1 (MQ=255)                                                                                                                                                                              gAAGCCACCACGGCCGCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCaac                                                                                                                                                                            <  1:163523/78‑3 (MQ=255)                                                                                                                                                                                   ccccacGGCCTCCGGCGTGGATAATTCAGCCTCCCCCCGACTGGCAGACACCGCGGAACGGTATTATTTCACCCTCAGAGAGAGGCTGATCATTATGCAAAAACAACTGGAATGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGATCTGACCATATCGATGGGCATCGCATGCAAGTATTTTGATCAATACACACGCTCTTCCAGCGGATTATAAATGCCTAAAGTAATAAAACCGAGc  <  1:68748/241‑1 (MQ=38)                                                                                                                                                                                         |                                                                                                                                                                                                                                             TAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGC  >  NC_000913/578243‑578663

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