Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,727,050	T→G	100%	noncoding (2135/2904 nt)	rrlG ←	23S ribosomal RNA of rrnG operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,727,050	0	T	G	54.1%	‑6.2 / 2.5	7	noncoding (2135/2904 nt)	rrlG	23S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (0/4);  total (1/6)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-02

CACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCAAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGG  >  NC_000913/2726955‑2727141                                                                                                |                                                                                            cgctttCCGGAACGCGGATCACGGTACAATTTTAGAACAGCAAACATTATAGGGTGGTATTTCAAGTTCGGCTCCATGCAGACTGGCGTCCACACTTCaa                                                                                         <  1:3107566/95‑1 (MQ=1)                                 tAGAACATCAAACATTAAAGGGTTGTATTTCAAGGTCGGCTCCATGCAGACTGGCTTCTTAATGACACGGCCACCAACCAGATCTACACATCAAGGCTCa                                                         <  1:3000595/100‑1 (MQ=0)                                                      gtcgtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCTACACGTCAAAGCCTCCCTCCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAg                                    <  1:1305626/97‑1 (MQ=1)                                                       tgttATTTCAAGGTCGGCTCCATGCAGACTGGCGTTCTAATGACAAGGCGTCCCCCGAGACCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGt                                   <  1:3201839/97‑1 (MQ=0)                                                       tgctATTTCAAGGTCGGCTCCATGCAGACTGTCGTACTCACTTCATAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCtata                                     <  1:70836/95‑1 (MQ=1)                                                        ggaaTTTCAAGGTCGGATCAATGCAGACTGGCGTCCACACTTCAAAGACTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTa                                  >  1:790760/4‑100 (MQ=1)                                                                                         tctacacGTCAAAGCCTCCCACCTAACCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCgg  <  1:149843/96‑1 (MQ=0)                                                                                                |                                                                                            CACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCAAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGG  >  NC_000913/2726955‑2727141

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