Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,628,614	A→G	7.7%	intergenic (‑84/+23)	rlmJ ← / ← putP	23S rRNA (adenine(2030)‑N(6))‑methyltransferase/sodium/L‑proline transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,628,614	0	A	G	7.7%	141.8 / 2.6	53	intergenic (‑84/+23)	rlmJ/putP	23S rRNA (adenine(2030)‑N(6))‑methyltransferase/sodium/L‑proline transporter
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base G (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

AGCACGATGTGTTTGAGGACGTCGGCGTGGTTGCCGGCGTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCAGAACAATGGCGATGCTGGCGAACAGGAAGCCCGGGATGA  >  NC_002947/5628469‑5628758                                                                                                                                                   |                                                                                                                                                 aGCACGATGTGTTTGAGGACGTCGGCGTGGTTGCCGGCGTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGgc                                                                                                                                               <  3:184681/149‑1 (MQ=255)     cGATGTGTTTGAGGACGTCGGCGTGGTTGCCGGCGTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGcgccg                                                                                                                                           <  5:110775/149‑1 (MQ=255)           gtTTGAGGACGTCGGCGTGGTTGCCGGCGTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAAt                                                                                                                                      >  8:209726/1‑148 (MQ=255)           gtTTGAGGACGTCGGCGTGGTTGCCGGCGTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAAt                                                                                                                                      >  2:320888/1‑148 (MQ=255)               gAGGACGTCGGCGTGGTTGCCGGCGTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCaa                                                                                                                                  <  6:63811/148‑1 (MQ=255)                                 gCCGGCGTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCg                                                                                                               >  4:29011/1‑149 (MQ=255)                                  ccGGCGTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCg                                                                                                               <  5:41227/148‑1 (MQ=255)                                  ccGGCGTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGt                                                                                                              <  3:29011/149‑1 (MQ=255)                                    ggCGTGGAAGGCGTGACGATAGTTCATGGCAACTTCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGt                                                                                                              >  1:43960/1‑148 (MQ=255)                                       gTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCGCCTGCTTTGCGGGTGCCGCCAAtc                                                                                                                                     <  4:167924/121‑2 (MQ=38)                                       gTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCGCCTGCTTTGCGGGTGCCGCCAAtc                                                                                                                                     <  4:173294/121‑2 (MQ=38)                                       gTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCGCCTGCTTTGCGGGTGCCGCCAAtc                                                                                                                                     >  3:173294/1‑120 (MQ=38)                                       gTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCGCCTGCTTTGCGGGTGCCGCCAAtc                                                                                                                                     >  3:167924/1‑120 (MQ=38)                                        tGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATAcg                                                                                                        >  1:6195/1‑149 (MQ=255)                                          gAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATAcgc                                                                                                       <  3:111854/148‑1 (MQ=255)                                               cGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATc                                                                                                 <  7:209726/149‑1 (MQ=255)                                                gTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCAt                                                                                                  >  4:285320/1‑147 (MQ=255)                                                   aCGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGc                                                                                              <  2:43960/148‑1 (MQ=255)                                                   aCGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGc                                                                                              <  1:249552/148‑1 (MQ=255)                                                     gATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGc                                                                                           >  5:280772/1‑149 (MQ=255)                                                         gTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTc                                                                                        >  1:77463/1‑148 (MQ=255)                                                          ttCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCg                                                                                       >  8:194472/1‑148 (MQ=255)                                                             aTGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAgc                                                                                   >  8:39340/1‑149 (MQ=255)                                                                 cAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCtt                                                                               <  1:320888/149‑1 (MQ=255)                                                                  aaC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCtt                                                                               >  8:241783/1‑148 (MQ=255)                                                                  aaC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCtt                                                                               >  5:314871/1‑148 (MQ=255)                                                                   aC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGa                                                                             <  3:285320/149‑1 (MQ=255)                                                                      tCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGAc                                                                            >  7:170898/1‑148 (MQ=255)                                                                                 cggcAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGc                                                                  >  8:217591/1‑147 (MQ=255)                                                                                 cggcAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCg                                                                 >  1:208543/1‑148 (MQ=255)                                                                                  ggcAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGaa                                                               <  3:114806/149‑1 (MQ=255)                                                                                   gcAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAg                                                              >  1:43455/1‑149 (MQ=255)                                                                                    cAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAgg                                                             <  7:39340/149‑1 (MQ=255)                                                                                       gTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTTATTGGCGGCACCCGCAAAGCTGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATAcgc                                                                                                       <  5:12770/104‑1 (MQ=255)                                                                                       gTTTACCTTGCATTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATAcgc                                                                                                       >  6:12770/1‑104 (MQ=255)                                                                                        tttACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAgg                                                                                                                                                >  1:100269/1‑62 (MQ=255)                                                                                        tttACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAgg                                                                                                                                                <  2:100269/62‑1 (MQ=255)                                                                                                gccttgcAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCt                                                  <  7:241783/148‑1 (MQ=255)                                                                                                   ttgcAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATAcgc                                                                                                       >  3:137528/1‑92 (MQ=255)                                                                                                   ttgcAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATAcgc                                                                                                       <  4:137528/92‑1 (MQ=255)                                                                                                       aGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGa                                                                             >  6:324874/1‑114 (MQ=255)                                                                                                       aGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGa                                                                             <  5:324874/114‑1 (MQ=255)                                                                                                        ggCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCa                                         <  6:314871/149‑1 (MQ=255)                                                                                                        ggCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCa                                         <  2:43455/149‑1 (MQ=255)                                                                                                           ggAGTCAGGGCTGGCTGCTGATTGGCGGCAACCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCAGaa                                      <  8:284649/149‑1 (MQ=255)                                                                                                                            ctgATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCAGAACAATGGCGATGCTGGc                      <  8:170898/148‑1 (MQ=255)                                                                                                                              gATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCAGaa                                      <  5:80106/130‑1 (MQ=255)                                                                                                                              gATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCAGaa                                      >  6:80106/1‑130 (MQ=255)                                                                                                                                ttGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCAGAACAATGGCGATGCTGGCGAACa                 <  2:6195/149‑1 (MQ=255)                                                                                                                                     ggcACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCAGAACAATGGCGATGCTGGCGAACAGGAAg            >  7:236991/1‑149 (MQ=255)                                                                                                                                           cGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCAGAACAATGGCGATGCTGGCGAACAGGAAGCCCgg       >  5:109996/1‑148 (MQ=255)                                                                                                                                            gCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCAGAACAATGGCGATGCTGGCAAACAGGAAGCCCGGga     >  5:176533/1‑149 (MQ=255)                                                                                                                                               aaGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCAGAACAATGGCGATGCTGGCGAACAGGAAGCCCGGgatga  <  2:208543/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AGCACGATGTGTTTGAGGACGTCGGCGTGGTTGCCGGCGTGGAAGGCGTGACGATAGTTCATGGCAAC‑TCCTGCGAAGGCGGCAAGTTTACCTTGCCTTGCAGGCGGAGTCAGGCCTGGCTGCTGATTGGCGGCACCCGCAAAGCAGGCGCCGCCAATTCAACGGGTGTTACTTGTCGGCGTGATACGCCGCATCAGCCGCTTCGAAGCGCTTGACCATCGTCTGCGAAGGGCTGCCCAGCTTGCTCACCAGAACAATGGCGATGCTGGCGAACAGGAAGCCCGGGATGA  >  NC_002947/5628469‑5628758

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