Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	579,385	0	A	G	56.2%	3.2 / 13.3	17	N70N (AAT→AAC)	ybcV	DUF1398 domain‑containing protein YbcV
Reads supporting (aligned to +/- strand):  ref base A (1/6);  new base G (5/4);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCG  >  NC_000913/579141‑579624                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                ccAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCtt                                                                                                                                                                                                                                           >  1:77958/1‑251 (MQ=255)                     cGGAAGAAACCGGACGTTATGATTTAGCGTGGGAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAAc                                                                                                                                                                                                                       >  1:66785/1‑251 (MQ=255)                      ggAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTTGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACCTTTACAAcc                                                                                                                                                                                                                      <  1:111454/251‑1 (MQ=255)                      ggAAGAAACCGAACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAAcc                                                                                                                                                                                                                      >  1:130883/1‑251 (MQ=255)                                     ttATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATCTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTAc                                                                                                                                                                                                                          <  1:126927/232‑1 (MQ=255)                                                                         gAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATGTGGAAGCCATCGGAAAACTCCTGCTTTAGCAAGATGTTCCCTGTATTGCTGAAATGTGATTTCTCATGATTTCAACCGATCATAGGACTTTTCTATAAGATGCTTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAAtgtg                                                                                                                                                                    <  1:65140/250‑1 (MQ=38)                                                                                                                                                aTCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTc                                                                                            <  1:2911/251‑1 (MQ=255)                                                                                                                                                                                   ccTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATAttt                                                         >  1:155594/1‑251 (MQ=255)                                                                                                                                                                                           gCTGAAATNTNNNTNNTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATNTCTTGAGAATTTAACNTTNACAANCTTNTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGc                                                                                                                                                     >  2:117032/1‑151 (MQ=255)                                                                                                                                                                                           gCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGCATTTAACATTTACAACCTTTTTACGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATATGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTaaa                                                  >  1:28245/1‑250 (MQ=255)                                                                                                                                                                                           gCTGAAATGTGATTTCTCTTGATTTCAACCTAGCATAGGACGTTGCTATAAGATGCGTATTTCTGGAGAAGTTAACATTTACAACCTTTTGAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGc                                                                                                                                                     <  1:117032/151‑1 (MQ=255)                                                                                                                                                                                              gAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATTCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTGATCGTTTTCTATCACGATGGGAAGATTATCTGTTGTTAGATAGTAAATATAATGTGAGTCGTTGTGACGTTTTAGTTCATAATAAACCAATTCACAGTCTAAATCttt                                                                             <  1:159827/220‑1 (MQ=255)                                                                                                                                                                                                            ctTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGGTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAAGGGCAACCTGAgn                                   <  1:40947/248‑2 (MQ=255)                                                                                                                                                                                                                               aGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTGACAACCTTTTTAAGTCCTTTTATTAGCACGGCGTTATCGTTTTCTAACACAATGTGAATATTATCGGTGGCTAGATAGTAAATATAATGTGAGACAGTGTGACGTGTGAGTTCAGAATAAAGCAAGTCGCAGTTTAAGTCTTTTCGCACTTGATCGAATATTTCTTTAAAAAGGGCAACCTGAGCCATTGGTAAAACCTTCCAtgtg             <  1:149884/251‑1 (MQ=255)                                                                                                                                                                                                                                   cGTTTCTATAAGCTGTGTCTTTCTTGAGAATTTAGCATTTACAACCTTTTTAAGTCATTTTATTAACACGGAGTTATCGTTTTCTAACACAATGTGCATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATAc         >  1:142199/1‑251 (MQ=255)                                                                                                                                                                                                                                   cGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGGGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTGAGTTCAGAAGAAAACAAGTCACAGTTTAAATCTTTTCGCTCGTTATCGATTATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATAc         <  1:123605/251‑1 (MQ=255)                                                                                                                                                                                                                                          ataAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGACTAGATAGTAAATATAATGTGAGACGGTGTGACGTTTTAGTTCAGAATAAAACAAGTCACAGGCTAAATCTTTTCGCACTTGATCGAATAGTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGgcg  <  1:167549/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                CCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCG  >  NC_000913/579141‑579624

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