Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,313,822	A→G	71.8%	noncoding (830/1522 nt)	PP_16SG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,313,822	0	A	G	71.8%	‑5.9 / 4.7	15	noncoding (830/1522 nt)	PP_16SG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (3/7);  total (5/10)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCATGTCAAGGCCTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGTCAACTTAATGCGTTAGCTGCGCCACTAAAATCTCAAGGATTCCAACGGCTAGTTGACATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAGT  >  NC_002947/5313673‑5313918                                                                                                                                                      |                                                                                                 tgaaTGTCAAGGCCTGGTAAGGTTATTCGCGTTGCTTCGAATTAAACCACATGCTACACCGCTTGTGCGGGCCCCAGTCAATTAATTTGAGTTATAACCGTGCGGCCGTACTCCCCAGGCGGTCAACTTAATGCGTTAGCTGAGACACTa                                                                                                  >  2:3592719/4‑150 (MQ=1)          agaccTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGTCAACTTAATGCGTTAGCTGCGCCACTAAAATCTCaa                                                                                         <  1:5106982/147‑1 (MQ=1)          agaccTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGTCAACTTAATGCGTTAGCTGCGCCACTAAAATCTCaa                                                                                         <  1:5107598/147‑1 (MQ=1)                                                                      gggccccGTCAATTCATTTTTGTTTTTACCCTGGGGGCGTTAGGCCCAGGGGGTCAAATTTATTGGTTTCGTTCGCCAACTAAATCTCAAGGATTCCAACGGCTAGTTGACATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGt                             <  2:1245774/147‑1 (MQ=0)                                                                                  tttaTTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGTCAACTTAATGCGTTAGCTGCGCCACTAAAATCTCAAGGATTCCAACGGCTAGTTGACATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACg                 >  1:877304/4‑150 (MQ=1)                                                                                  tttaTTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGTCAACTTAATGCGTTAGCTGCGCCACTAAAATCTCAAGGATTCCAACGGCTAGTTGACATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACg                 >  1:877554/4‑150 (MQ=1)                                                                                   tcctttGAGTTTTAATCTTGCGACCGTACTCCCCAGGCGGAATGTTTAATGCGTTAGCTGCGCCACCGAAGAGT‑AAACTCCCCGACGGCTAACATTCATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCt               >  1:997566/4‑150 (MQ=0)                                                                                   tcctttGAGTTTTAATCTTGCGACCGTACTCCCCAGGCGGAATGTTTAATGCGTTAGCTGCGCCACCGAAGAGT‑AAACTCCCCGACGGCTAACATTCATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCt               >  1:997567/4‑150 (MQ=0)                                                                                                aatcTTGCGACCGTACTCCCCAGGCGGAATGTTTAATGCGTTAGCTGCGCCACCGAAGAGT‑AAACTCCCCGACGGCTAACATTCATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAgc  <  1:4753406/147‑2 (MQ=0)                                                                                                aatcTTGCGACCGTACTCCCCAGGCGGAATGTTTAATGCGTTAGCTGCGCCACCGAAGAGT‑AAACTCCCCGACGGCTAACATTCATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAgc  <  1:5454707/147‑2 (MQ=0)                                                                                                aatcTTGCGACCGTACTCCCCAGGCGGAATGTTTAATGCGTTAGCTGCGCCACCGAAGAGT‑AAACTCCCCGACGGCTAACATTCATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAgc  <  1:564988/147‑2 (MQ=0)                                                                                                aatcTTGCGACCGTACTCCCCAGGCGGAATGTTTAATGCGTTAGCTGCGCCACCGAAGAGT‑AAACTCCCCGACGGCTAACATTCATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAgc  <  2:1583613/147‑2 (MQ=0)                                                                                                aatcTTGCGACCGTACTCCCCAGGCGGAATGTTTAATGCGTTAGCTGCGCCACCGAAGAGT‑AAACTCCCCGACGGCTAACATTCATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAgc  <  2:1583634/147‑2 (MQ=0)                                                                                                aatcTTGCGACCGTACTCCCCAGGCGGAATGTTTAATGCGTTAGCTGCGCCACCGAAGAGT‑AAACTCCCCGACGGCTAACATTCATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAgc  <  2:2317786/147‑2 (MQ=0)                                                                                                aatcTTGCGACCGTACTCCCCAGGCGGAATGTTTAATGCGTTAGCTGCGCCACCGAAGAGT‑AAACTCCCCGACGGCTAACATTCATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAgc  <  2:6935629/147‑2 (MQ=0)                                                                                                aatcTTGCGACCGTACTCCCCAGGCGGAATGTTTAATGCGTTAGCTGCGCCACCGAAGAGT‑AAACTCCCCGACGGCTAACATTCATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAgc  >  2:8251984/4‑149 (MQ=0)                                                                                                                                                      |                                                                                                 TGCATGTCAAGGCCTGGTAAGGTTCTTCGCGTTGCTTCGAATTAAACCACATGCTCCACCGCTTGTGCGGGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGTCAACTTAATGCGTTAGCTGCGCCACTAAAATCTCAAGGATTCCAACGGCTAGTTGACATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAGT  >  NC_002947/5313673‑5313918

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