Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,310,134	C→G	34.9%	noncoding (2426/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,134	0	C	G	34.9%	‑3.9 / 32.8	37	noncoding (2426/2902 nt)	PP_23SG	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (6/18);  new base G (4/9);  total (10/27)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAGCGC  >  NC_002947/5309990‑5310279                                                                                                                                                 |                                                                                                                                                  caacTCACGTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCtttt                                                                                                                                              <  2:3674680/147‑1 (MQ=1) caacTCACGTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCtttt                                                                                                                                              <  2:3673770/147‑1 (MQ=1)     tcacgTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGTGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATcc                                                                                                                                          >  2:4705111/4‑150 (MQ=1)     tcacgTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATcc                                                                                                                                          <  2:235276/147‑1 (MQ=1)     tcacgTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATcc                                                                                                                                          <  2:235258/147‑1 (MQ=1)     tcacgTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATcc                                                                                                                                          >  1:6163179/4‑150 (MQ=1)            ccgcTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGc                                                                                                                                   <  2:5864205/147‑1 (MQ=1)            ccgcTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGc                                                                                                                                   >  1:1746696/4‑150 (MQ=1)            ccgcTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGc                                                                                                                                   <  1:988398/147‑1 (MQ=1)                 ttaatgggCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGTGATgg                                                                                                                              <  1:1238696/144‑1 (MQ=1)                 ttaatgggCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGTGATgg                                                                                                                              <  1:1238828/144‑1 (MQ=1)                 ttaatgggCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATgg                                                                                                                              <  2:7640347/144‑1 (MQ=1)                 ttaatgggCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATgg                                                                                                                              <  2:3759565/144‑1 (MQ=1)                                    acgcTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAaccac                                                                                                           <  2:1955606/147‑1 (MQ=1)                                              cctgCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTa                                                                                                 >  2:5153013/4‑150 (MQ=1)                                              cctgCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTa                                                                                                 >  1:5077992/4‑150 (MQ=1)                                               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        >  1:4385741/4‑146 (MQ=0)                                                                                                                                           ggcgTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCt        <  2:4385741/143‑1 (MQ=0)                                                                                                                                             aggaCCTTTTAGCCGTTGCGCGACGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAgcgc  <  2:5133511/147‑1 (MQ=0)                                                                                                                                                 |                                                                                                                                                  CAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAGCGC  >  NC_002947/5309990‑5310279

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