Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,314,314	C→T	30.4%	noncoding (338/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,314,314	0	C	T	30.4%	‑1.9 / 14.1	23	noncoding (338/1522 nt)	PP_16SG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (5/11);  new base T (4/3);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

GTCGGTAACGTCAAAACAGCAAGGTATTAACTTACTGCCCTTCCTCCCAACTTAAAGTGCTTTACAATCCGAAGACCTTCTTCACACACGCGGCATGGCTGGATCAGGCTTTCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGCCATTAC  >  NC_002947/5314181‑5314413                                                                                                                                      |                                                                                                    ctctgtaaagtaaaAACAGCAAGGTATTAACTTACTGCCCTTCCTCCCAACATAAAGTGCTTTACAATCCGAAGACCTTCTTCACACACGCGGCATGGCTGGATCAGGCTTTCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTa                                                                                     <  2:4314339/146‑1 (MQ=0)                                                        agagCTTTACAACCCTAGGGCCTTCATCACTCACGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAgct                              >  1:1607834/4‑148 (MQ=1)                                                                        aaaaCCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTg              <  2:3517362/147‑1 (MQ=0)                                                                        aaaaCCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTg              >  2:2503381/4‑150 (MQ=0)                                                                        aaaaCCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTg              >  2:2503354/4‑150 (MQ=0)                                                                            ccgtcttcACACCCGCGGCATGGCTGGATCAGGCTTTCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGa          <  1:221046/147‑1 (MQ=1)                                                                               tcatcactcacGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGcc       >  1:353759/4‑150 (MQ=0)                                                                               tcatcactcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGcc       <  1:577680/147‑1 (MQ=1)                                                                               tcatcactcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGcc       >  1:976611/4‑150 (MQ=1)                                                                               tcatcactcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGcc       <  1:577660/147‑1 (MQ=1)                                                                               tcatcactcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGcc       <  1:1992990/147‑1 (MQ=1)                                                                               tcatcactcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCAGCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGcc       >  2:4382511/4‑150 (MQ=1)                                                                                    actcacGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCGTTAc  <  1:57466/147‑1 (MQ=0)                                                                                    actcacGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCGTTAc  >  1:4615606/4‑150 (MQ=0)                                                                                    actcacGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTATCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCGTTAc  <  2:1540004/147‑1 (MQ=0)                                                                                    actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc  <  1:1055266/147‑1 (MQ=1)                                                                                    actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc  >  1:5067203/4‑150 (MQ=1)                                                                                    actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc  >  2:2412892/4‑150 (MQ=1)                                                                                    actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc  <  1:4296356/147‑1 (MQ=1)                                                                                    actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc  <  2:398785/147‑1 (MQ=1)                                                                                    actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc  <  2:4304367/147‑1 (MQ=1)                                                                                    actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc  <  1:1786854/147‑1 (MQ=1)                                                                                    actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc  <  2:4578191/147‑1 (MQ=1)                                                                                                                                      |                                                                                                    GTCGGTAACGTCAAAACAGCAAGGTATTAACTTACTGCCCTTCCTCCCAACTTAAAGTGCTTTACAATCCGAAGACCTTCTTCACACACGCGGCATGGCTGGATCAGGCTTTCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGCCATTAC  >  NC_002947/5314181‑5314413

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