Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,636,503	0	T	C	55.9%	1.5 / 42.1	34	intergenic (‑136/‑253)	nohA/ynfO	putative prophage DNA‑packaging protein NohA/DUF3950 domain‑containing protein YnfO
Reads supporting (aligned to +/- strand):  ref base T (7/8);  new base C (9/10);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAAC  >  NC_000913/1636266‑1636752                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                          gCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCgg                                                                                                                                                                                                                                               >  2:164833/1‑250 (MQ=255)       gCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAAttt                                                                                                                                                                                                                                        <  1:164833/251‑1 (MQ=255)                                    gCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTc                                                                                                                                                                                                             <  2:232482/249‑1 (MQ=255)                                             gAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTa                                                                                                                                                                                                  <  2:130572/251‑1 (MQ=255)                                                             cGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCGGCCCCGTAACCTGTCGGATCAcc                                                                                                                                                                                   <  1:90768/250‑1 (MQ=11)                                                                            agccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTaaa                                                                                                                                                                   <  2:64122/251‑1 (MQ=255)                                                                                   tgctTTTTGTTGCCTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATa                                                                                                                                                             >  2:18611/1‑250 (MQ=11)                                                                                    gctTTTTGTTGACTTACATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGaca                                                                                                                                                             <  2:75405/249‑3 (MQ=11)                                                                                          tgttgACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTat                                                                                                                                                     >  1:171706/1‑251 (MQ=255)                                                                                             tgACTTCCATTGTTCATTCCACGGCCAAAACCAGAGAAAGGAAACGACAGAGGCACAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAAAATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATACTGATTatcat                                                                                                                                                  >  1:116898/1‑251 (MQ=255)                                                                                                                           aCAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCCCCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGGAACCTGTCGGAGCACCGGAAAGGACCCGTAAAGTGATAATGATTATCAGCTACAGATCTCAACGTGCTTGGAGGCCATc                                                                                                                    <  2:23547/251‑1 (MQ=14)                                                                                                                              gagaAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAACTAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCaa                                                                                                                  >  2:121183/1‑250 (MQ=255)                                                                                                                                gaAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCTCAACGTGCGTGGAGGCCATCAAAcc                                                                                                               <  1:205743/251‑1 (MQ=255)                                                                                                                                 aaaGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGACAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAAcc                                                                                                               >  2:199016/1‑250 (MQ=255)                                                                                                                                       aaCGACAGAGGCCAAAAATCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCa                                                                                                         >  2:218106/1‑250 (MQ=255)                                                                                                                                           aCAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAaata                                                                                                     >  1:204108/1‑250 (MQ=255)                                                                                                                                                  cccAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTa                                                                                             <  2:3695/251‑1 (MQ=255)                                                                                                                                                      aaaGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGAc                                                                                         <  2:74065/251‑1 (MQ=255)                                                                                                                                                          cTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGc                                                                                       >  2:131152/1‑249 (MQ=255)                                                                                                                                                                ttCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATc                                                                                <  1:157105/250‑1 (MQ=255)                                                                                                                                                                     cACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTTATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTa                                                                          <  2:234647/251‑1 (MQ=255)                                                                                                                                                                      aCCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTaa                                                                         <  1:225233/251‑1 (MQ=255)                                                                                                                                                                          gTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAAtt                                                                       <  2:4116/249‑1 (MQ=255)                                                                                                                                                                                tCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCTCCGGAAAgg                                                                                                                                                                            <  2:159687/142‑1 (MQ=11)                                                                                                                                                                                tCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAgg                                                                                                                                                                            >  1:159687/1‑142 (MQ=255)                                                                                                                                                                                        tttCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACtt                                                        <  1:110066/250‑1 (MQ=255)                                                                                                                                                                                                ggTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTaaaa                                               <  1:24741/251‑1 (MQ=255)                                                                                                                                                                                                      ttAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCa                                                                                                         <  2:38360/187‑1 (MQ=255)                                                                                                                                                                                                      ttAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCa                                                                                                         >  1:38360/1‑187 (MQ=255)                                                                                                                                                                                                              aaaCATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACa                                  >  1:78816/1‑250 (MQ=255)                                                                                                                                                                                                               aaCATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCAATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTACAAACAACTTCAGACAAt                                >  1:151379/1‑251 (MQ=255)                                                                                                                                                                                                                                aagaagaaCGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTg               >  1:188511/1‑251 (MQ=255)                                                                                                                                                                                                                                   aagaaCGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTACTTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAAt            >  1:31993/1‑251 (MQ=255)                                                                                                                                                                                                                                             aTGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAAc  >  2:149802/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                          GCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAAC  >  NC_000913/1636266‑1636752

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