Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,166,856	0	G	A	47.1%	1.8 / 15.9	17	noncoding (198/1542 nt)	rrsB	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base A (4/4);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACACGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCGGAT  >  NC_000913/4166613‑4166887                                                                                                                                                                                                                                                    |                                acacGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAAAGGGGGa                          <  1:427305/251‑1 (MQ=11)  cacGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAAAGGGGGAc                         <  1:472892/251‑1 (MQ=14)      tAATTCATTACGAAGTTTAATTCTTGGAGCGACAAACTTTTAAATTGCAGAGTTTGATCATGGCTCAGATTGAACGCTGACGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGAGGCGGACGGGTGAGTAATGTATGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGATAATACCGCATAACGTCGCAAGACCAAAGGGGGGGAccgtc                     >  1:354758/1‑248 (MQ=11)        aTTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGAACAAAAAGGGGGACCTTCgg                   >  2:120792/1‑251 (MQ=14)        aTTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTATAAAGTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAAAGGGGGACCGTCgg                   <  2:13855/251‑1 (MQ=9)        aTTCATTACGAAGATAAATTCTTTGAGCGTCAACCTTTTCAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGAAATGTCTGGGAAACTGCCTGATGGAGGGGGGTAACTACTGGAAACGGTAGCTAATACCGGATAACGTCGCAAGACCAAAGAGGGGGACCTTCgg                   >  2:29642/1‑251 (MQ=11)           cATTACGAAGTTTCATTCTTTGAGCGTCAAACTTCTAAATTGAAGAGTTTGATCATGGTTCAGATTGCACGCTGGCGGCAGGCCTAAAACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGCCGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCGGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAAAGGGGGACCTGCGGGcc                >  2:88452/1‑251 (MQ=11)               aCGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCtt            <  2:457764/251‑1 (MQ=31)                cGAAGTTTAATTCTTTGCGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAAAGGGGGACCTTCGGGCCTCTTg           >  2:569586/1‑251 (MQ=14)                cGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGGGGGGGGTAACTACTGGAAACGGTAGCTAATACCGCATAAAGTCGCCAGACCAAAAAGGGGGACCTTCGGGCCTCTTg           >  2:188368/1‑251 (MQ=9)                cGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTg           >  1:254579/1‑251 (MQ=31)                 gAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGc          >  1:36349/1‑251 (MQ=31)                 gAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAAAGGGGGACCTTCGGGCCTCTTGc          <  1:318487/251‑1 (MQ=14)                   aGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCa        <  1:111683/251‑1 (MQ=31)                      ttAATTCTTTGAGCGTCAAACTGTTAAATTGAAGAGTTTGATTATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACAGGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATATCGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATc      <  1:256117/250‑1 (MQ=9)                      ttAATCCTTTGAGCGTCAAACTTTTAAATAGAAGAGTTTGATCGTGGCTCAGATGGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAATCTTGCTTCTTTTGTGACGAGTGGCGTACGGGTGAGTAATGTCTTGGAAACGGCCTGATGGAGGGCGATAACTACTGGAAACGGTAGCTAATACCGTATAACGTCGCAAGTCCAAAGAGGGGGACCTTCGGGCCTCTTGCCATc      <  1:374604/250‑1 (MQ=11)                         aTTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCGGAt  <  1:247134/251‑1 (MQ=31)                                                                                                                                                                                                                                                    |                                ACACGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCGGAT  >  NC_000913/4166613‑4166887

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