Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,309,926	C→A	50.0%	noncoding (2634/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,309,926	0	C	A	50.0%	‑6.9 / 35.1	24	noncoding (2634/2902 nt)	PP_23SG	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (8/4);  new base A (6/6);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGC  >  NC_002947/5309819‑5310045                                                                                                            |                                                                                                                        gccttCAGCGGTTATCCCGTCCGTATATAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAAc                                                                               <  2:1860025/147‑1 (MQ=0) gccttCAGCGGTTATCCCGTCCGTATATAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTAC‑ACCACGGCAGATAgg                                                                                       <  2:1003934/138‑1 (MQ=0)                           atcgCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCt                                                     <  1:4315587/147‑1 (MQ=1)                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc                                             >  2:83155/7‑150 (MQ=0)                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc                                             >  1:2577022/7‑150 (MQ=0)                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc                                             >  2:1424112/7‑150 (MQ=0)                                   ccagcgaTGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGtct                                                                          <  1:5120757/115‑1 (MQ=0)                                   ccagcgaTGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGtct                                                                          >  2:5120757/7‑121 (MQ=0)                                   ccagcgaTGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTAcc                                             >  1:1215697/7‑150 (MQ=0)                                   ccagcgaTGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTAcc                                             >  1:1945981/7‑150 (MQ=0)                                   ccagcgaTGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTAcc                                             >  2:739549/7‑150 (MQ=0)                                     tgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAccgc                                           >  1:3078986/5‑148 (MQ=0)                                     tgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAccgc                                           <  2:2855878/146‑3 (MQ=0)                                     agcgaTGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTAccgc                                           >  1:4363884/5‑148 (MQ=0)                                     agcgaTGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTAccgc                                           <  2:4003112/146‑3 (MQ=0)                                                                     agcggTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGAcc           <  1:4415391/147‑1 (MQ=0)                                                                     agcggTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGAcc           >  2:2588718/4‑150 (MQ=0)                                                                     agcggTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGAcc           >  2:2588748/4‑150 (MQ=0)                                                                      gagatacGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCTACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGAcct          >  1:3431238/7‑149 (MQ=0)                                                                              tccatcccGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGc  >  1:4716168/7‑150 (MQ=0)                                                                              tccatcccGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGc  <  1:4767860/144‑1 (MQ=0)                                                                              tccatcccGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGc  <  1:4767989/144‑1 (MQ=0)                                                                              tccatcccGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGc  <  2:1832891/144‑1 (MQ=0)                                                                              tccatcccGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGc  >  1:1582440/7‑150 (MQ=0)                                                                                                            |                                                                                                                        GCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGC  >  NC_002947/5309819‑5310045

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