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

TTCGACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCAC  >  NC_002947/5309750‑5310004                                                                                                                                            |                                                                                                                    ttagaaggCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCTTCTTTAGGCAAGTTTCCCTATTAGTTGCTTTCAGCGTTTATGTCTTACTTACATAGCTACCAGGAAATGCCACTGGAGTGTCAACAGGAACAACAGAGGTTCTTCCa                                                                                                           >  2:2755621/7‑150 (MQ=1)                                                                      gccttCAGCGGTTATCCCGTCCGTATATAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAAc                                      <  2:1860025/147‑1 (MQ=0)                                                                      gccttCAGCGGTTATCCCGTCCGTATATAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTAC‑ACCACGGCAGATAgg                                              <  2:1003934/138‑1 (MQ=0)                                                                                                atcgCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCt            <  1:4315587/147‑1 (MQ=1)                                                                                                        cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    >  1:2577022/7‑150 (MQ=0)                                                                                                        cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    >  2:1424112/7‑150 (MQ=0)                                                                                                        cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    >  2:83155/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)                                                                                                                                            |                                                                                                                    TTCGACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCAC  >  NC_002947/5309750‑5310004

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑