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

GTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTT  >  NC_000913/579271‑579645                                                                                                                                                                                                                                        |                                                                                                                                                gTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTAATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAg                                                                                                                              <  1:45787/251‑1 (MQ=255)                    acaggcaCTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTGGATAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGtt                                                                                                           <  1:102553/245‑1 (MQ=255)                     ggAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGtct                                                                                                          >  1:75445/1‑249 (MQ=255)                      gAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTTACGTTTTAGTTCAGAATAAAACAATTCACAGTTTa                                                                                                         <  1:5222/251‑1 (MQ=255)                                                cccTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTCTCATAGGACGTCTCTATAAGATGCGTATTTCTTGAAACTTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGCGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATACAACAATTCACAGTTTAAATc                                                                                                     >  1:101939/1‑229 (MQ=255)                                                                          ctTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAAGTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGc                                                        <  1:113068/248‑1 (MQ=255)                                                                                   aaCCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGATGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAAcc                                                             >  1:223355/1‑234 (MQ=255)                                                                                                           aGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTGGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGt                    <  1:214816/251‑1 (MQ=255)                                                                                                                cgtgtttCTTGAGAATTTAACATTTACAACCTTGTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATTTGAATATAATCGGTGGCTAGATAGTAAAGATAATGTGAGACGTTGGGACGTTTTAGTTCAGAATAAAACAATTcac                                                                                                               <  1:108618/151‑1 (MQ=255)                                                                                                                     ttCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAGTATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGGTGAAAGCTGTTCGCACGGGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGGGATACGAGGGCGCGTAGTTTGCAt           <  1:217388/250‑1 (MQ=39)                                                                                                                             aaTTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACTTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGtttt  <  1:151806/251‑1 (MQ=255)                                                                                                                                        tACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAAAGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATGTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTc                                        <  1:213351/202‑1 (MQ=255)                                                                                                                                             ncTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATCTGTATATTATATGTGACTAGATCGTAAATCTAATGTGAGACGTTGTGACGTTTTAATTCCGAATACAACAATTCACAGACTACATCTATTCGCACTTGCTCGAATATTTCTTTAAAACTAGCAACCTGAGCCATTGGTAAAACCTTCAATGTGATACGAGGGCGCGTAGTTTGCa            >  1:20014/2‑225 (MQ=255)                                                                                                                                                                     cGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGtttn              <  1:52393/199‑2 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                GTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTT  >  NC_000913/579271‑579645

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: ‑