Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,620,859	0	C	T	55.6%	‑2.5 / 20.2	18	R556R (CGC→CGT)	rhsB	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base T (7/3);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GAGCCGTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAAC  >  NC_000913/3620790‑3621141                                                                      |                                                                                                                                                                                                                                                                                            gAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAATTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCAGGTGACCTGAACGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcgc                                                       >  1:408388/1‑299 (MQ=11)      gTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGATTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAg                                                   <  2:4940/298‑1 (MQ=14)      gTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTa                                                 >  1:452808/1‑300 (MQ=31)              cAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGAAACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGa                                             >  2:281133/1‑296 (MQ=14)               aattGCTTAGCTTCACCGACTGCTCCGTTTAGGTAACCCGGTGTGACCATTACCGCTTCGGGCAAACGACGGTGGTGCACCGCGAGGAAGGGTTGAGTCAGTGCCGCGCATGCGATAGCCGTGGACAGGTAATTGCC‑TGTAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCTTCATTGCCCCGGACTGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATAc                                         <  2:55885/297‑1 (MQ=11)               aGTTGCTGAGCTTCACCGACTGTTCCTGTTATGTAACCCGTTATGACCATGACCGTTTCGGGCAGATGACGGCGGTGCATCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAACCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAATACGGTATGAATACAACGTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCATAAACTGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATAc                                         <  2:55897/299‑1 (MQ=11)               aGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATAc                                         <  2:707099/299‑1 (MQ=31)                gTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAACACGa                                       >  2:442793/1‑300 (MQ=255)                   gctgaactccgcggACTGTTCTTGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACTGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCAGGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGCGGAAAGGCCGTCCGTACCACGCAGGGCTGGCTGACGCGCAGTATGGAATACAAtgc                                    <  2:356476/287‑1 (MQ=11)                             actgACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACgg                          <  2:635090/297‑1 (MQ=9)                                        cGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTAGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCCGTATGGAATACGATGCTTCCGGACGGGTCATCCGcat               >  2:616475/1‑298 (MQ=11)                                         ggTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGcc                <  1:281133/298‑1 (MQ=31)                                          gTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTCATTGCCGTG‑AAAGACACGCAGGGCCATGACACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGNGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGCTGCTGCCGGACGGGTCATCCGCCTGa             >  1:518236/1‑300 (MQ=11)                                          gTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGNGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGa             >  1:154168/1‑300 (MQ=9)                                          gTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGAAGGCAGCAGAAAAGGGACACAGTACGATGCGTGGNGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGACTGa             >  1:517080/1‑300 (MQ=11)                                             aTGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCa          >  2:414802/1‑300 (MQ=31)                                             aTGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCa          >  1:59882/1‑300 (MQ=9)                                                     cGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCAGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCAGTCCGTCCCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAAc  >  2:343302/1‑300 (MQ=14)                                                                      |                                                                                                                                                                                                                                                                                            GAGCCGTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAAC  >  NC_000913/3620790‑3621141

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