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

CCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGG  >  NC_000913/2726540‑2726687                     |                                                                                                                                cccctttCAATTCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGAc         <  1:1962240/135‑1 (MQ=1)     tctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCg     >  1:637278/8‑141 (MQ=1)      ctcaactctcTAACGTCCACGGAATATAGGGACCTAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGTCCTTCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  2:209286/135‑1 (MQ=1)      ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  2:890118/135‑1 (MQ=1)      ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  2:53590/135‑1 (MQ=1)      ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:547010/135‑1 (MQ=1)      ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  2:480272/135‑1 (MQ=1)      ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  2:1877066/135‑1 (MQ=1)      ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  2:1726903/135‑1 (MQ=1)      ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:548015/135‑1 (MQ=1)      ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:2341400/135‑1 (MQ=1)        caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:1895883/138‑1 (MQ=0)        caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:517459/138‑1 (MQ=0)        caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:258604/138‑1 (MQ=0)        caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:2523114/138‑1 (MQ=0)        caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:1977741/138‑1 (MQ=0)        caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:1977641/138‑1 (MQ=0)        caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:1977583/138‑1 (MQ=0)        caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  2:917543/138‑1 (MQ=0)                     |                                                                                                                                CCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGG  >  NC_000913/2726540‑2726687

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