Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,310,082	A→C	43.1%	noncoding (2478/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,082	0	A	C	43.1%	‑5.4 / 43.1	35	noncoding (2478/2902 nt)	PP_23SG	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (7/13);  new base C (2/13);  total (9/26)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

CAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGAC  >  NC_002947/5309933‑5310209                                                                                                                                                      |                                                                                                                                cacatctcaaaCGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCACGCGTACCCCTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAaca                                                                                                                                 <  2:1039809/147‑1 (MQ=0)         aatcGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTc                                                                                                                         >  1:3309900/4‑150 (MQ=1)            cgcccACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGat                                                                                                                      <  1:2153106/147‑1 (MQ=1)            cgcccACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGat                                                                                                                      <  1:2846902/147‑1 (MQ=1)            cgcccACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGat                                                                                                                      <  1:3486470/147‑1 (MQ=1)                                                          caacTCACGTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCtttt                                                                        <  2:3673770/147‑1 (MQ=1)                                                          caacTCACGTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCtttt                                                                        <  2:3674680/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                                                             >  1:1746696/4‑150 (MQ=1)                                                                     ccgcTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGc                                                             <  2:5864205/147‑1 (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:3759565/144‑1 (MQ=1)                                                                          ttaatgggCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATgg                                                        <  2:7640347/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                          <  2:2081096/144‑1 (MQ=1)                                                                                                        cgactacAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATGCGGAACCACCGGATCACTaa                          <  1:748232/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               >  2:2519133/4‑150 (MQ=1)                                                                                                                   cctggATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTc               <  1:2536936/147‑1 (MQ=1)                                                                                                                   cctggATGTGATGAGCCGACAGCGAGGTGCCACACACCGCCGTCGATATGAACTCTTGGGCGGTTTCTGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTc               <  2:5649341/147‑1 (MQ=1)                                                                                                                       gttgcgatgaGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGAGTGTTGGGGGTGATCAGCCTGTTATCCCCGGCGTAGCTTTTATCCGTTGAGCGATGGGCCTTGCACGCGGGACCACCGGATGACTATGTCCGTCTTTCGtct           >  2:6580598/6‑148 (MQ=0)                                                                                                                         tgcgatgaGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGtctct         <  1:5657304/147‑5 (MQ=1)                                                                                                                                agtcGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGAc  <  1:1060911/147‑1 (MQ=1)                                                                                                                                agtcGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGAc  <  2:4584520/147‑1 (MQ=1)                                                                                                                                agtcGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGAc  <  2:5724984/147‑1 (MQ=1)                                                                                                                                                      |                                                                                                                                CAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGAC  >  NC_002947/5309933‑5310209

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