Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,309,858	A→T	47.0%	noncoding (2702/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,858	0	A	T	47.0%	‑5.5 / 14.2	17	noncoding (2702/2902 nt)	PP_23SG	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base T (4/3);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACACCCAACCTATCAACGTCGTAGTCTTCGACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAAC  >  NC_002947/5309722‑5309988                                                                                                                                         |                                                                                                                                   acccacCCAACCTATCAACGTCGTAGTCTTCGACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCGTAGATGCTTTCAGCGGTTATCTCTTCCGACCATAGCTACCCGGCAATGCCACTGGCGt                                                                                                                       <  2:244265/147‑1 (MQ=0)                             tttgACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCa                                                                                           >  2:3157443/4‑150 (MQ=1)                             ttagACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCa                                                                                           >  2:4577448/4‑150 (MQ=1)                                     cccataaGGGGATTCAAGATCCCAGAGAGATCTCATCTTGAGGCAAGTTTCCCGCATAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTACGTCCACTCCGGTc                                                                                   >  2:720086/7‑150 (MQ=1)                                                                 gaactcATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCtctcaa                                                       <  2:6834304/147‑1 (MQ=1)                                                                         ctagAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCtctcaaatctcaaa                                               <  2:6687002/147‑1 (MQ=1)                                                                                                  gccttCAGCGGTTATCCCGTCCGTATATAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAAc                      >  2:7953280/4‑150 (MQ=0)                                                                                                  gccttCAGCGGTTATCCCGTCCGTATATAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAAc                      <  1:627228/147‑1 (MQ=0)                                                                                                  gccttCAGCGGTTATCCCGTCCGTATATAGCTAACCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAAc                      >  2:3337634/4‑150 (MQ=0)                                                                                                       cagcacttATCCCTTCCGCACATAGCTACCCAGCGATGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGtct                 <  2:6230250/144‑1 (MQ=0)                                                                                                       cagcacttATCCCTTCCGCACATAGCTACCCAGCGATGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGtct                 <  2:6230280/144‑1 (MQ=0)                                                                                                       cagcacttATCCCTTCCGCACATAGCTACCCAGCGATGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGtct                 <  1:7072608/144‑1 (MQ=0)                                                                                                               atcccgtCCGTATATAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGt         >  2:7348559/7‑150 (MQ=0)                                                                                                                tccctTCCGCACATAGCTACCCAGCGATGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGtt        <  2:4545261/147‑1 (MQ=0)                                                                                                                      ccgtataTAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACt                     >  2:6464269/7‑131 (MQ=0)                                                                                                                      ccgtataTAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACt                     <  1:6464269/125‑1 (MQ=0)                                                                                                                      ccgtataTAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAAc  <  1:6886888/144‑1 (MQ=0)                                                                                                                                         |                                                                                                                                   ACACACCCAACCTATCAACGTCGTAGTCTTCGACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAAC  >  NC_002947/5309722‑5309988

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