Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,726,652	0	A	G	50.0%	‑4.5 / 27.0	20	noncoding (2533/2904 nt)	rrlG	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (3/7);  new base G (0/10);  total (3/17)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAAC  >  NC_000913/2726526‑2726695                                                                                                                               |                                            cggaCTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGa                               <  2:602118/138‑1 (MQ=1)                    ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa            <  1:180411/135‑1 (MQ=1)                    ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa            <  1:180541/135‑1 (MQ=1)                    ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa            <  1:1123780/135‑1 (MQ=1)                    ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa            <  1:1827587/135‑1 (MQ=1)                    ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa            <  1:519132/135‑1 (MQ=1)                    ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa            <  2:527921/135‑1 (MQ=1)                    ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa            <  2:1899806/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)                      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)                           ctcaaacgcCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCa     >  1:831561/7‑141 (MQ=0)                              caacgcCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAac  >  1:742338/4‑141 (MQ=1)                              caacgcCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAac  >  1:742301/4‑141 (MQ=1)                                                                                                                               |                                            CGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAAC  >  NC_000913/2726526‑2726695

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