Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,390,913	G→A	9.8%	M577I (ATG→ATA)	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,913	0	G	A	9.8%	140.7 / 5.2	54	M577I (ATG→ATA)	rlmKL	23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (25/22);  new base A (3/2);  total (29/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.84e-01

AGCGGCAAAGCGGCACGCGCCAGTACGAACGTCAGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGACTGGGCACGGCGTAACCT  >  NC_002947/2390768‑2391062                                                                                                                                                  |                                                                                                                                                        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                                                                                                              >  1:1131664/1‑139 (MQ=255)                                                 gTTTCCAGGAAGTGAACGAAGGTGGCGTCAGGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCt                                                                                                              <  2:1131664/139‑1 (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                                                                                                                >  2:136276/1‑94 (MQ=255)                                                                                            ccTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAAc                                                                                                                <  1:136276/94‑1 (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                                                                                 >  2:831422/1‑99 (MQ=255)                                                                                                                      tGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGc                                                                                 <  1:831422/99‑1 (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                                                                             >  2:180783/3‑86 (MQ=255)                                                                                                                                       tggATGCGCATAC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCAcgc                                                                             <  1:180783/84‑1 (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:1047965/150‑1 (MQ=255)                                                                                                                                           tGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGACTGGGCACGGc         >  1:369886/1‑150 (MQ=255)                                                                                                                                           tGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGACTGGGCACGGc         <  1:451117/150‑1 (MQ=255)                                                                                                                                                  gc‑gc‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCGACGGCCACCGTAAAAAAACCGAGAGGCGGTGCGCGCAGCACCACCAGCGTCGCCCTGTCGAAAACCTACCTGGACTGGGCACGGCGTAACCt  >  1:449178/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                        AGCGGCAAAGCGGCACGCGCCAGTACGAACGTCAGGCTACCGAGGGCCGTTTCCAGGAAGTGAACGAAGGTGGCGTCAAGCTGCTGGTCAACCTCACCGACTACCTGGATACCGGGCTGTTCCTCGACCACCGCCCGATGCGTATGC‑GC‑ATCCAGCGCGAGGCTGCCGGCAAGCGCTTCCTCAACCTGTTCTGCTACACCGCCACGGCCACCGTGCACGCGGCCAAGGGCGGTGCGCGCAGCACCACCAGCGTCGACCTGTCGAAAACCTACCTGGACTGGGCACGGCGTAACCT  >  NC_002947/2390768‑2391062

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