Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,726,555	0	C	A	55.6%	‑4.2 / 23.5	18	noncoding (2630/2904 nt)	rrlG	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/8);  new base A (0/10);  total (0/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 = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGTGATGCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGG  >  NC_000913/2726499‑2726687                                                         |                                                                                                                                     cccgTGACGCGCCCCCCCCCGTCCTCTCGCACTCGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATa                                                  <  1:1082952/138‑1 (MQ=1)                            cggaCTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGa                       <  2:602118/138‑1 (MQ=1)                                               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:1123780/135‑1 (MQ=1)                                               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:180411/135‑1 (MQ=1)                                               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:180541/135‑1 (MQ=1)                                               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:1827587/135‑1 (MQ=1)                                               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:519132/135‑1 (MQ=1)                                               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:758845/135‑1 (MQ=1)                                               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  1:763920/135‑1 (MQ=1)                                               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  2:1690532/135‑1 (MQ=1)                                               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  2:1899806/135‑1 (MQ=1)                                               ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa    <  2:527921/135‑1 (MQ=1)                                                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCTCTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  2:718822/138‑1 (MQ=0)                                                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAACTGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:371910/138‑1 (MQ=0)                                                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  2:98012/138‑1 (MQ=1)                                                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:617295/138‑1 (MQ=0)                                                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:371939/138‑1 (MQ=0)                                                 caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg  <  1:1817540/138‑1 (MQ=0)                                                         |                                                                                                                                     CCAGTGATGCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGG  >  NC_000913/2726499‑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: ‑