Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,310,055	T→C	79.5%	noncoding (2505/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,055	0	T	C	79.5%	0.8 / 11.9	35	noncoding (2505/2902 nt)	PP_23SG	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base C (5/22);  total (8/27)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTC  >  NC_002947/5309911‑5310196                                                                                                                                                 |                                                                                                                                              tcttACTAGGAGCAGCCCATCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGACTAACGAAGTTCAAAACACAGATCTCGTACCACTATAAATGGAGAACAGCCATACCCTTGGGACCGGATTCAGCCCCAGGATGtgatga                                                                                                                                          >  2:6164102/4‑150 (MQ=1)               gctccTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGt                                                                                                                            <  2:8387243/147‑1 (MQ=1)               gctccTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGt                                                                                                                            <  2:4924099/147‑1 (MQ=1)               gctccTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGt                                                                                                                            <  1:5103618/147‑1 (MQ=1)               gctccTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGt                                                                                                                            <  2:4924346/147‑1 (MQ=1)                       cacatctcaaaCGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCACGCGTACCCCTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAaca                                                                                                                    <  2:1039809/147‑1 (MQ=0)                               aatcGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTc                                                                                                            >  1:3309900/4‑150 (MQ=1)                                  cgcccACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGat                                                                                                         <  1:2846902/147‑1 (MQ=1)                                  cgcccACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGat                                                                                                         <  1:3486470/147‑1 (MQ=1)                                  cgcccACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTACCGCTTTAATGGGCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGat                                                                                                         <  1:2153106/147‑1 (MQ=1)                                                                                caacTCACGTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCtttt                                                           <  2:3674680/147‑1 (MQ=1)                                                                                caacTCACGTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCtttt                                                           <  2:3673770/147‑1 (MQ=1)                                                                                    tcacgTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGTGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATcc                                                       >  2:4705111/4‑150 (MQ=1)                                                                                    tcacgTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATcc                                                       >  1:6163179/4‑150 (MQ=1)                                                                                    tcacgTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATcc                                                       <  2:235258/147‑1 (MQ=1)                                                                                    tcacgTACCGCTTTAAATGGCGAACAGCCATACCCTTGGGACCTGCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATcc                                                       <  2:235276/147‑1 (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:1238828/144‑1 (MQ=1)                                                                                                ttaatgggCGAACAGCCCAACCCTTGGGACCGACTACAGCCCCAGGATGCGATGAGCCGACATCGAGGTGCCAAACCTCCCCGTCGATGTGGACTCTTGGGGGAGATAAGCCTGTTATCCCCGGGGTAGCTTTTATCCGTTGAGTGATgg                                           <  1:1238696/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              >  1:5077992/4‑150 (MQ=1)                                                                                                                             cctgCTCCAGCCCCAGGATGCGATGAGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTa              >  2:5153013/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)                                                                                                                                                 |                                                                                                                                              TCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTC  >  NC_002947/5309911‑5310196

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