Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,410,053	(T)9→8	intergenic (+47/‑80)	rlmB → / → yjfI	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,410,045	0	T	.	100.0%	115.9 / NA	25	intergenic (+39/‑88)	rlmB/yjfI	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (5/20);  total (5/20)

TGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGA  >  NC_000913/4409895‑4410160                                                                                                                                                       |                                                                                                                    tGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                                                     >  2:482898/1‑151 (MQ=255)           gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                                           >  1:486682/1‑148 (MQ=255)              gAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCta                                                                                                        <  1:290332/150‑1 (MQ=255)                        gCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAg                                                                                             <  2:648508/151‑1 (MQ=255)                         cATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACAGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCTGTTTTGGATGG‑TTTTTTTTGTCTATAGCTgg                                                                                                 >  1:665074/1‑146 (MQ=255)                             ccGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCTGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAAt                                                                                        <  1:109075/151‑1 (MQ=255)                              cGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                                                       <  1:282113/151‑1 (MQ=255)                                  ggCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                                    <  2:410856/150‑1 (MQ=255)                                  ggCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                                   <  1:672573/151‑1 (MQ=255)                                      ggAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                                                 <  2:475143/149‑1 (MQ=255)                                      ggAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTg                                                                               <  2:414187/151‑1 (MQ=255)                                        aaGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                                             <  1:603130/151‑1 (MQ=255)                                        aaGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                                             <  1:687844/151‑1 (MQ=255)                                        aaGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                                             <  1:239059/151‑1 (MQ=255)                                                  tcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                   >  1:390163/1‑151 (MQ=255)                                                   cccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                   <  1:318636/150‑1 (MQ=255)                                                        aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                                             <  1:68057/151‑1 (MQ=255)                                                          cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                                          <  1:461829/151‑1 (MQ=255)                                                           gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                                          <  2:206220/151‑1 (MQ=255)                                                             ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat                                                         <  1:674967/150‑1 (MQ=255)                                                                ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTa                                                      >  2:157668/1‑150 (MQ=255)                                                                    gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                                                 <  2:665321/151‑1 (MQ=255)                                                                                              ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt                       <  2:196095/151‑1 (MQ=255)                                                                                                  cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGaa                   <  2:542735/151‑1 (MQ=255)                                                                                                   gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAAt                  >  2:242890/1‑151 (MQ=255)                                                                                                              tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGacgc        >  1:678055/1‑150 (MQ=255)                                                                                                                      cAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGa  <  2:472010/148‑1 (MQ=255)                                                                                                                                                       |                                                                                                                    TGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGA  >  NC_000913/4409895‑4410160

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