Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,310,228	A→G	43.7%	noncoding (2332/2902 nt)	PP_23SG ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,310,228	0	A	G	43.7%	‑4.7 / 14.3	16	noncoding (2332/2902 nt)	PP_23SG	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (2/7);  new base G (3/4);  total (5/11)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAGCGCACCTTCGTACTCCTCCGTTACTCTTTGGGAGGAGACCGCCCCAGTCAAACTACCCACCATACACT  >  NC_002947/5310092‑5310344                                                                                                                                         |                                                                                                                     atgtggaCTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTAAGCCCGACTTTCGTCCCTGCTCGACTTGTAGGTCTCGCAGTCAAGCTCCCTTGTGcc                                                                                                         <  1:4416009/144‑1 (MQ=1)  tatgggctctTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGcca                                                                                                        <  2:4217627/144‑2 (MQ=0)   gtggaCTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTAAGCCCGACTTTCGTCCCTGCTCGACTTGTAGGTCTCGCAGTCAAGCTCCCTTGTGCCtt                                                                                                       <  1:3096634/146‑1 (MQ=1)    tggaCTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTAAGCCCGACTTTCGTCCCTGCTCGACTTGTAGGTCTCGCAGTCAAGCTCCCTTGTGCCttt                                                                                                      <  2:2159042/147‑1 (MQ=1)    tggaCTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTAAGCCCGACTTTCGTCCCTGCTCGACTTGTAGGTCTCGCAGTCAAGCTCCCTTGTGCCttt                                                                                                      <  2:2159138/147‑1 (MQ=1)                       aggcTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATAGAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTc                                                                                   <  2:953395/147‑1 (MQ=1)                                     ggggTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTAAGCCCGACTTTCGTCCCTGCTCGACTTGTAGGTCTCGCAGTCAAGCTCCCTTGTGCCTTTACACTCTACGAATGATTTCCAAccatt                                                                           >  1:3836799/4‑141 (MQ=1)                                     ggggTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTAAGCCCGACTTTCGTCCCTGCTCGACTTGTAGGTCTCGCAGTCAAGCTCCCTTGTGCCTTTACACTCTACGAATGATTTCCAAccatt                                                                           <  2:3836799/141‑4 (MQ=1)                                     ggggTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTAAGCCCGACTTTCGTCCCTGCTCGACTTGTAGGTCTCGCAGTCAAGCTCCCTTGTGCCTTTACACTCTACGAATGATTTCCAACCATTCTGAgg                                                                     <  1:1238194/147‑2 (MQ=0)                                     ggggTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTAAGCCCGACTTTCGTCCCTGCTCGACTTGTAGGTCTCGCAGTCAAGCTCCCTTGTGCCTTTACACTCTACGAATGATTTCCAACCATTCTGAgg                                                                     >  2:441582/4‑149 (MQ=0)                                     ggcgTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGCCCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAgc                                                                     <  1:2522104/147‑1 (MQ=0)                                     ggcgTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAgc                                                                     <  1:2808383/147‑1 (MQ=0)                                                                                                 ccgaCTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTACccacc         >  2:981829/4‑150 (MQ=0)                                                                                                 ccgaCTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTACccacc         >  2:981859/4‑150 (MQ=0)                                                                                                 ccgaCTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCGGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTACccacc         <  2:1453266/147‑1 (MQ=0)                                                                                                        tcgtctcTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTACCCACCATacacg  >  1:497697/7‑149 (MQ=0)                                                                                                                                         |                                                                                                                     ATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAGCGCACCTTCGTACTCCTCCGTTACTCTTTGGGAGGAGACCGCCCCAGTCAAACTACCCACCATACACT  >  NC_002947/5310092‑5310344

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