Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,726,555	0	C	A	52.6%	‑4.5 / 27.7	19	noncoding (2630/2904 nt)	rrlG	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (1/8);  new base A (0/10);  total (1/18)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGG  >  NC_000913/2726531‑2726687                         |                                                                                                                                     tatgAGCAGCCCACCTCAGTTCTCCAGCGCACACGGCAGATAGGGACCGAACAGTCTCACGACGATCTAAAACCAGCTCGCGTACCACTTAAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGtga                  >  2:684054/4‑141 (MQ=1)       cacccccccTCAGATCTCCAGCGCCCACCGCAGACAGCGACCGAACTGTCACCCCACGTACTAAACCCAGCTCGCGTACCACTTCACATCGCGAACCGCCATACCCTTGGGACCTACTTCACCCCCAGGATGTGATGAGcc            <  2:1169812/138‑1 (MQ=1)               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  2:964881/135‑1 (MQ=1)               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  2:1311449/135‑1 (MQ=1)               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:1987103/135‑1 (MQ=1)               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:1987098/135‑1 (MQ=1)               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:1986890/135‑1 (MQ=1)               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:1815812/135‑1 (MQ=1)               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:1517586/135‑1 (MQ=1)               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:1392356/135‑1 (MQ=1)               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:1793006/135‑1 (MQ=1)               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:1792980/135‑1 (MQ=1)                 caattCTCCAACGCCCACGGCTGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTCGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:962264/138‑1 (MQ=1)                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  2:1115238/138‑1 (MQ=1)                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:1935958/138‑1 (MQ=0)                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:1003451/138‑1 (MQ=0)                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:2042961/138‑1 (MQ=0)                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:858792/138‑1 (MQ=0)                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:931481/138‑1 (MQ=0)                         |                                                                                                                                     TAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGG  >  NC_000913/2726531‑2726687

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