Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	579,099	0	A	.	55.6%	5.9 / 15.9	10	intergenic (‑206/+85)	borD/ybcV	prophage lipoprotein BorD/DUF1398 domain‑containing protein YbcV
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base . (1/4);  total (2/8)
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.37e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGAT  >  NC_000913/578894‑579314                                                                                                                                                                                                              |                                                                                                                                                                                                                        ggTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTAGGTCAGGTGCGATGAATCGTCATTGTATTCCCGTATGAACTAGGTCCACAGCCGTGACGGTGAACTTCGCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCaa                                                                                                                                                                            <  1:96034/250‑1 (MQ=255)                    tttGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACgg                                                                                                                                                         <  1:118554/250‑1 (MQ=255)                          gTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCGAAACAATGAGTTGAAATTTCATAGTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAgc                                                                                                                                                                                                     <  1:11142/201‑1 (MQ=255)                             cgccgaatcatatttatAGAATATTAAACAAACCCTAAACAATTATTTGATATTTCATATTGTTAATATTTATTAATGTTTGCCATGTTCTATGCATCTTCATTGTATTCCCGTATTAACTATGTCCACAGCCCTTACGGGGCACTTCTCTGCGGGAGTGTACGGGAATAATTAAAACATATTCACACATGTTTTAGCGCGTTCATGTATTGTATTATGCCAACCCCCCGGTGCTTCCACGgatgaaccc                                                                                                                                                 <  1:13247/234‑4 (MQ=255)                              gCAGAATCGTATGTGTAGAAAATTAAACAAACCCTACACATTGGGTTGAAATTTCCTATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCCTTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGCAATTCTTTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAAAGCCCCGGTGCTGACACGGAAGAAACCgg                                                                                                                                               >  1:168228/1‑250 (MQ=255)                                            gtAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATGAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGnc                                                                                                                                                                              <  1:136029/206‑1 (MQ=255)                                                                              aaaTTTCATATTGTTAATATTTATTAATGTATGCCAGGTACCATGAATCGTCCTTGTATTGCCGGATTAACTATGTACACAGCCTTGACGGGCAACTTCTCTGAGGGCGTGTCCGGGAATACTTGAAAACGATGCACACCGGGGTTAGTGCGTACATGTATTGTATTATGca                                                                                                                                                                              >  1:176046/1‑171 (MQ=255)                                                                                                                            aTCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGAGTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATg                                                 <  1:119496/251‑1 (MQ=255)                                                                                                                                                                             cTTCTCTGCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGAt  <  1:76256/248‑1 (MQ=37)                                                                                                                                                                                            ggccGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTAGTGCATTATGc                                                                                                                                                                               <  1:117177/58‑1 (MQ=255)                                                                                                                                                                                                              |                                                                                                                                                                                                                        GGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGAT  >  NC_000913/578894‑579314

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