Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,390,910	T→C	9.1%	R576R (CGT→CGC)	rlmKL →	23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,390,910	0	T	C	9.1%	146.9 / 4.4	55	R576R (CGT→CGC)	rlmKL	23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (27/23);  new base C (3/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TCGCGAGCGGCAAAGCGGCACGCGCCAGTACGAACGTCAGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGACTGGGCACGGC  >  NC_002947/2390763‑2391055                                                                                                                                                    |                                                                                                                                                    tCGCGAGCGGCAAAGCGGCACGCGCCAGTACGAACGTCAGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTAt                                                                                                                                                   >  1:579084/1‑150 (MQ=255)      aGCGGCAAAGCGGCACGCGCCAGTACGAACGTCAGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑a                                                                                                                                            <  2:579084/150‑1 (MQ=255)       tcGGCAAAGCGGCACGCGCCAGTACGAACGTCAGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGTGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑At                                                                                                                                           <  1:532758/149‑1 (MQ=255)        cGGCAAAGCGGCACGCGCCAGTACGAACGTCAGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATc                                                                                                                                          <  2:413415/150‑1 (MQ=255)         ggCAAAGCGGCACGCGCCAGTACGAACGTCAGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTCTTCCTCGACCACCTCCCGATGCGTATCA‑GC‑ATcc                                                                                                                                         >  2:532758/1‑150 (MQ=255)                                gaccGTCAGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAgcgc                                                                                                                  >  1:704936/4‑150 (MQ=255)                                 aaCGTCAGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCt                                                                                                                 >  2:1177946/1‑150 (MQ=255)                                       aGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCa                                                                                                           >  2:569107/1‑150 (MQ=255)                                         gCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAAc                                                                                                         <  2:824049/150‑1 (MQ=255)                                         gCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAAc                                                                                                         >  2:451117/1‑150 (MQ=255)                                           tACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAgg                                                                                                                                >  1:872934/1‑125 (MQ=255)                                           tACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAgg                                                                                                                                <  2:872934/125‑1 (MQ=255)                                           tACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCt                                                                                                       >  1:502111/1‑150 (MQ=255)                                                aGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCt                                                                                                  <  2:818203/150‑1 (MQ=255)                                                aGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCt                                                                                                  <  2:502111/150‑1 (MQ=255)                                                   gCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCt                                                                                               >  2:525561/1‑150 (MQ=255)                                                      gTTTCCAGGAAGTGAACGAAGGTGGCGTCAGGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCt                                                                                                       <  2:1131664/139‑1 (MQ=255)                                                      gTTTCCAGGAAGTGAACGAAGGTGGCGTCAGGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCt                                                                                                       >  1:1131664/1‑139 (MQ=255)                                                                agggAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGc                                                                                  <  1:123475/147‑1 (MQ=255)                                                                    aaCGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCAcc                                                                              <  2:823289/150‑1 (MQ=255)                                                                       gAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTaca                                                                                            <  2:8260/133‑1 (MQ=255)                                                                       gAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTaca                                                                                            >  1:8260/1‑133 (MQ=255)                                                                               cGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGc                                                                   >  2:548497/1‑150 (MQ=255)                                                                                    aGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACC‑TCG‑CATATTCGTTTCCTGCTTTCCATCTCTAGGCTGCCTGCAAGCGCTTCCTCACCCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAgg                                                              <  1:1177946/150‑1 (MQ=255)                                                                                     gctgctGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAggg                                                             >  2:817476/1‑150 (MQ=255)                                                                                        gctgGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCgg                                                          <  2:704936/150‑1 (MQ=255)                                                                                            gTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTgcg                                                      >  2:884711/1‑150 (MQ=255)                                                                                                 ccTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAAc                                                                                                         <  1:136276/94‑1 (MQ=255)                                                                                                 ccTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAAc                                                                                                         >  2:136276/1‑94 (MQ=255)                                                                                                 ccTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCgcag                                                 >  2:624116/1‑150 (MQ=255)                                                                                                    cACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACg                                                                                    <  2:8320/112‑1 (MQ=255)                                                                                                    cACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACg                                                                                    >  1:8320/1‑112 (MQ=255)                                                                                                     aCCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGcacc                                             <  1:548497/150‑1 (MQ=255)                                                                                                      ccGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGcacca                                            >  2:334702/1‑150 (MQ=255)                                                                                                       cGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGcaccac                                           >  1:679425/1‑150 (MQ=255)                                                                                                         aCTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCAccacca                                         <  1:525561/150‑1 (MQ=255)                                                                                                                 gATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGAc                                 <  1:884711/150‑1 (MQ=255)                                                                                                                    aCCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTg                              >  1:191504/1‑150 (MQ=255)                                                                                                                      cGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTc                            >  1:171439/1‑150 (MQ=255)                                                                                                                           tGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGc                                                                          <  1:831422/99‑1 (MQ=255)                                                                                                                           tGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGc                                                                          >  2:831422/1‑99 (MQ=255)                                                                                                                                 tCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAacctacc                 >  2:1149335/1‑150 (MQ=255)                                                                                                                                      caccGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGAc            >  1:32509/1‑150 (MQ=255)                                                                                                                                      caccGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCACAATATACGCTGCACGCGGCCAAGGGCGGTGCGAGCAGCACCACCAGCGTCGACGTGACGAAAACCTACCTGGAc            >  1:812000/1‑150 (MQ=255)                                                                                                                                         cGCCCGATGCGTATGC‑AC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGt                                                                            >  2:568541/1‑83 (MQ=255)                                                                                                                                         cGCCCGATGCGTATGC‑AC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGt                                                                            <  1:568541/83‑1 (MQ=255)                                                                                                                                          gCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGACTggg        <  1:334702/150‑1 (MQ=255)                                                                                                                                            tggATGCGCATAC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCAcgc                                                                      <  1:180783/84‑1 (MQ=255)                                                                                                                                            tggATGCGCATAC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCAcgc                                                                      >  2:180783/3‑86 (MQ=255)                                                                                                                                            tggATGCGCATAC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGt                                                         >  1:476060/3‑99 (MQ=255)                                                                                                                                            tggATGCGCATAC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGt                                                         <  2:476060/97‑1 (MQ=255)                                                                                                                                            tggATGCGCATAC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGACTGGGCa      >  1:237244/3‑150 (MQ=255)                                                                                                                                                tGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGACTGGGCACGGc  >  1:369886/1‑150 (MQ=255)                                                                                                                                                tGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGACTGGGCACGGc  <  1:451117/150‑1 (MQ=255)                                                                                                                                                tGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGACTGGGCACGGc  <  1:1047965/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TCGCGAGCGGCAAAGCGGCACGCGCCAGTACGAACGTCAGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGACTGGGCACGGC  >  NC_002947/2390763‑2391055

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