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%	79.2 / NA	19	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 . (7/12);  total (7/12)

CGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACT  >  NC_000913/4409924‑4410112                                                                                                                          |                                                                    cGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCt                                                   >  1:803170/1‑139 (MQ=255)  gATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttt                                                                 >  2:248697/1‑125 (MQ=255)  gATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttt                                                                 <  1:248697/125‑1 (MQ=255)   aTGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTgg                                                 >  1:137091/1‑139 (MQ=255)        tGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGa                                            <  1:448613/139‑1 (MQ=255)                       ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                             <  1:789600/139‑1 (MQ=255)                           aaCGTTTCGGCTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGa                                            <  1:201189/120‑1 (MQ=255)                           aaCGTTTCGGCTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGa                                            >  2:201189/1‑120 (MQ=255)                               tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                     <  1:133001/139‑1 (MQ=255)                                ttCGGTTGCGACCGGAATTTGGTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                    <  1:494148/139‑1 (MQ=255)                                ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttt                                                                >  1:596019/1‑96 (MQ=255)                                ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttt                                                                <  2:596019/96‑1 (MQ=255)                                       gCGACCGGATTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta             <  2:562472/139‑1 (MQ=255)                                            cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt        <  2:137091/139‑1 (MQ=255)                                                  ttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt  <  2:803170/139‑1 (MQ=255)                                                                   tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGa                                            >  1:646911/1‑80 (MQ=255)                                                                   tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGa                                            <  2:646911/80‑1 (MQ=255)                                                                   tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                             <  1:395617/95‑1 (MQ=255)                                                                   tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                             >  2:395617/1‑95 (MQ=255)                                                                               gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttt                                                                >  1:202799/1‑49 (MQ=255)                                                                               gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttt                                                                <  2:202799/49‑1 (MQ=255)                                                                                                 tGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat         <  1:137815/85‑1 (MQ=255)                                                                                                 tGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat         >  2:137815/1‑85 (MQ=255)                                                                                                          aTCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                              <  2:497482/55‑1 (MQ=255)                                                                                                          aTCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                              >  1:497482/1‑55 (MQ=255)                                                                                                                          |                                                                    CGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACT  >  NC_000913/4409924‑4410112

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