Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,313,849	T→A	71.4%	noncoding (803/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,849	0	T	A	71.4%	‑6.2 / 6.9	14	noncoding (803/1522 nt)	PP_16SG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base A (3/7);  total (5/9)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

GGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGTCAACTTAATGCGTTAGCTGCGCCACTAAAATCTCAAGGATTCCAACGGCTAGTTGACATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAGTGTCAGTATCAGTCCAGGTGGTCGCCTTCGCCACTGGTGTTCCTTCCTATATCTACGCATTTCACCGC  >  NC_002947/5313742‑5313985                                                                                                            |                                                                                                                                         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)                                                                                               tctaaCGGCTAGTTGACATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAGTGTCAGTATCAGTCCAGGTGGTCGCCTTCGCCACTGGTGTTCCTTCCTATATCTACGCATTTCACCGc  <  2:4667527/147‑1 (MQ=1)                                                                                                            |                                                                                                                                         GGCCCCCGTCAATTCATTTGAGTTTTAACCTTGCGGCCGTACTCCCCAGGCGGTCAACTTAATGCGTTAGCTGCGCCACTAAAATCTCAAGGATTCCAACGGCTAGTTGACATCGTTTACGGCGTGGACTACCAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAGTGTCAGTATCAGTCCAGGTGGTCGCCTTCGCCACTGGTGTTCCTTCCTATATCTACGCATTTCACCGC  >  NC_002947/5313742‑5313985

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