Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	578,811	0	T	G	47.4%	1.2 / 14.5	19	Q28P (CAG→CCG)	borD	prophage lipoprotein BorD
Reads supporting (aligned to +/- strand):  ref base T (5/5);  new base G (4/5);  total (9/10)
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.

AGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCT‑GGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCG  >  NC_000913/578569‑579034                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                agttgcccATAGCGATGGGCAATTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATGTACGAATGTTTGCT‑GTGTTTCTGTTTTAACAACATTTTCTGCTCCGCCACAAGTCTTGGCTGCATCGGCAGGTTTCTTCTTCCCAATGCCAGAAGCTAAGAAATGATGTGGGATGGTTTCCTTTGGTGCTACTGCTGCCGgtttgttt                                                                                                                                                                                                                         <  1:94048/245‑1 (MQ=255)     gACCATTTCTATGGGCAACTCATGCAATTATGTTGAGCAATACGCGCGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCT‑GTGTTTTTGGTTTAACACCATGTTCTGCGCCGCCACAAGTTGTAGCTGCAGCGACAGTTTTCTTCTACCCAATTCCAGAAACGAAGAAATGATGTGTTTTGGGTTCCTTTGGTTCTACTGCTGTCTGTTTGTTTTg                                                                                                                                                                                                                       <  1:150054/249‑1 (MQ=255)                  ggCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCT‑GGGTTTCTGTTTTCACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTg                                                                                                                                                                                                        >  1:86629/1‑251 (MQ=255)                                     gtGAGCAATACACCCGCTCTTCCATCGGAGTAAAAAAGCCTAAAGTAATAAAACCGAGTCATCCATTTACGAATGTTTGCTTGGGTTT‑TTTTTTAACAACATTTTCTGCGCCGCCACAAATTGTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGTTGGGTGATGTTTTCCTTTGGTTCTACTTCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGtaat                                                                                                                                                                                      <  1:149513/249‑1 (MQ=255)                                        aGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCT‑GGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGc                                                                                                                                                                                  <  1:50681/251‑1 (MQ=255)                                               cacacGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCT‑GGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGcc                                                                                                                                                                            >  1:76203/1‑250 (MQ=255)                                                   ngcgcTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCT‑GGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTATGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAgcgc                                                                                                                                                                       >  1:135311/2‑251 (MQ=255)                                                     cgcTTCCCGCGGCGTATAAATGCATAACGTAATACAACCGAGCAATCCATTTACGAATGCTTGCA‑GGGTTTCTGTTTTAACAACATTTTCTGCGCAGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGCAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAgcgc                                                                                                                                                                       >  1:141553/1‑249 (MQ=255)                                                       nttCCAGCGGAGTATAAATGACTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCT‑GGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGGAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGt                                                                                                                                                                    >  1:47771/2‑250 (MQ=255)                                                                gAGTATAAATGCCTAAATGAATAAAACCGAGCACTCCATTTAAGAATGTTTGCT‑GGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCGAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAACTGTCTGTTTAGCACATCCTGTAATAAGCAGGGCCAGAGCAGTAGCGAGTAGc                                                                                                                                                          >  1:134099/1‑251 (MQ=255)                                                                                      aaaaCCTAGCAATCCATTTACGAATGTTTGCT‑GGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATGCCAGAAACGAAGAAATGATTGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTGGAACAGTAAACGTCTGTTGATCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTTTTATTc                                                                                                                                     <  1:28061/250‑1 (MQ=255)                                                                                           cGAGCAATCCATTTACGAATGTTTGCG‑GGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATg                                                                                                                               <  1:117119/251‑1 (MQ=255)                                                                                                         aCGAATGTTTGCT‑GGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATGTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTc                                                                                                                  <  1:64462/250‑1 (MQ=255)                                                                                                                    cT‑GGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTGTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTg                                                                                                                                                                                                        <  1:20966/153‑1 (MQ=255)                                                                                                                                                                   ncTGCATCGCCAGTTTTCTTATGACCAATTCCAGCACCGAAGAAATGATGGGTGATGGTTTCCTTTGGTGATACTGCTGTATGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGACGTTCGCAGAATCGTATGTGAAGAAAATTAAACAAAACCTCAACCATGAGTTGAAATTTCatt                                                        >  1:48870/2‑250 (MQ=255)                                                                                                                                                                                       ctGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGCTTTGAACAGTAAACGTCTGTTGCGCACATCCTGTAATAAGAAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTAGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGt                                    >  1:117714/1‑251 (MQ=255)                                                                                                                                                                                          nccAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATg                                 >  1:17769/2‑251 (MQ=255)                                                                                                                                                                                                        caaataaaTGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCATCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCTTATTTGTAGAAAATTAAACAAACCCTAAAc                                                                           <  1:141044/190‑1 (MQ=255)                                                                                                                                                                                                                         aTGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGTATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCg  <  1:115715/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                AGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCT‑GGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCG  >  NC_000913/578569‑579034

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