Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,310,185	A→T	42.9%	noncoding (2375/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,185	0	A	T	42.9%	‑5.0 / 22.5	21	noncoding (2375/2902 nt)	PP_23SG	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (4/8);  new base T (2/7);  total (6/15)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAGCGCACCTTCGTACTCCTCCGTTACTCTTTG  >  NC_002947/5310036‑5310306                                                                                                                                                      |                                                                                                                          cgactacAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTaa                                                                                                                           <  1:748232/144‑1 (MQ=1) cgactacAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTaa                                                                                                                           <  2:2081096/144‑1 (MQ=1)  gactacAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTAAg                                                                                                                          >  2:9040615/6‑150 (MQ=1)   gcgtCAGCCCCAGGATGGGCTGAGCCGCCAGCGAGGTGCCACCCCCCGCCGTCGAGATGAACTCTTGGGCGGCATCCGCCTGTTATCCCCGGAGTACCTTTTAGCCGATGAGCGCGGGCCCTTCCATACAGAACCACCGGATCACTAAGa                                                                                                                         <  1:7420054/147‑1 (MQ=1)   actacAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTAAgc                                                                                                                         <  1:2790585/146‑2 (MQ=1)   actacAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTAAgc                                                                                                                         <  2:8007705/146‑2 (MQ=1)    ctccAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGAc                                                                                                                        <  2:8438695/147‑1 (MQ=1)            cctggATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTc                                                                                                                <  1:2536936/147‑1 (MQ=1)            cctggATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTc                                                                                                                >  2:2519133/4‑150 (MQ=1)            cctggATGTGATGAGCCGACAGCGAGGTGCCACACACCGCCGTCGATATGAACTCTTGGGCGGTTTCTGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTc                                                                                                                <  2:5649341/147‑1 (MQ=1)                gttgcgatgaGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGAGTGTTGGGGGTGATCAGCCTGTTATCCCCGGCGTAGCTTTTATCCGTTGAGCGATGGGCCTTGCACGCGGGACCACCGGATGACTATGTCCGTCTTTCGtct                                                                                                            >  2:6580598/6‑148 (MQ=0)                  tgcgatgaGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGtctct                                                                                                          <  1:5657304/147‑5 (MQ=1)                         agtcGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGAc                                                                                                   <  2:4584520/147‑1 (MQ=1)                         agtcGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGAc                                                                                                   <  2:5724984/147‑1 (MQ=1)                         agtcGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGAc                                                                                                   <  1:1060911/147‑1 (MQ=1)                                                             gatctctTGGGCGGTATCAGCCTGTTATCCCCGGAGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTta                                                               >  2:1919709/4‑150 (MQ=1)                                                                               agacTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTc                                             <  1:7651898/147‑1 (MQ=1)                                                                                             ggcgTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCTGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAgc                               <  2:8243483/147‑1 (MQ=0)                                                                                             ggcgTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCt                                   <  2:4385741/143‑1 (MQ=0)                                                                                             ggcgTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCt                                   >  1:4385741/4‑146 (MQ=0)                                                                                               aggaCCTTTTAGCCGTTGCGCGACGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAgcgc                             <  2:5133511/147‑1 (MQ=0)                                                                                                                gaccGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATAATCTACGTCCGATTTCCGACCGGTCTGAGCGCTCCTTCGTACTCCTCCg            >  2:8216657/4‑150 (MQ=0)                                                                                                                          ccgtCCAGACAGCACCCCCGGATCACGAAGACCTACTTTCGTACCGGCTCGACGTGTATGTCTCGCAGTCAAGCGCGCGTTTGCCTTCATACTCTCCGACCGATTTCCGACCGGTCTGAGCGCACCTTCGTACTCCTCCGGTACTCttgg  <  2:3505395/147‑3 (MQ=1)                                                                                                                                                      |                                                                                                                          CGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAGCGCACCTTCGTACTCCTCCGTTACTCTTTG  >  NC_002947/5310036‑5310306

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