Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,644,370	(G)5→6	intergenic (‑87/+63)	rlmN ← / ← ndk	23S rRNA m(2)A2503 methyltransferase/tRNA m(2)A37 methyltransferase/nucleoside diphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,644,365	1	.	G	100.0%	94.7 / NA	27	intergenic (‑82/+68)	rlmN/ndk	23S rRNA m(2)A2503 methyltransferase/tRNA m(2)A37 methyltransferase/nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (17/10);  total (17/10)

TGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCC‑GGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGAGCGGCAGATTCG  >  NC_000913/2644229‑2644499                                                                                                                                          |                                                                                                                                       tGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCC‑ggg                                                                                                                                     <  1:331852/140‑1 (MQ=255)             cTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCAtt                                                                                                                          >  2:878168/1‑140 (MQ=255)                        aCGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCt                                                                                                               >  1:879031/1‑140 (MQ=255)                        aCGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAAGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCt                                                                                                               >  1:499429/1‑140 (MQ=255)                                          aaTTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGt                                                                                             <  1:878168/140‑1 (MQ=255)                                            ttGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGttt                                                                                           >  1:335691/1‑140 (MQ=255)                                                 cAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCAcg                                                                                      <  1:903449/140‑1 (MQ=255)                                                         tGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCAttt                                                                              >  2:333588/1‑140 (MQ=255)                                                             aTTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGa                                                                          <  2:383354/140‑1 (MQ=255)                                                                   gCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACg                                                                           <  2:934390/133‑1 (MQ=255)                                                                   gCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACg                                                                           >  1:934390/1‑133 (MQ=255)                                                                       cGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAAc                                                                <  1:284656/140‑1 (MQ=255)                                                                                  gTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCggg                                                     <  2:879031/140‑1 (MQ=255)                                                                                  gTTATGGCCCCTGGAAGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCggg                                                     <  2:499429/140‑1 (MQ=255)                                                                                    tATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGca                                                   >  2:542187/1‑140 (MQ=255)                                                                                     aTGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGcac                                                  >  2:62002/1‑140 (MQ=255)                                                                                      tGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGcaca                                                 >  1:28564/1‑140 (MQ=255)                                                                                        gcaccTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACAct                                               <  2:303563/137‑1 (MQ=255)                                                                                             tGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACActtcgcc                                          <  1:62002/140‑1 (MQ=255)                                                                                                    ttGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGTGGGCACACTTCGCCTTCGCCa                                   >  1:38357/1‑140 (MQ=255)                                                                                                    ttGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCa                                   >  2:149022/1‑140 (MQ=255)                                                                                                    ttGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCa                                   >  1:50444/1‑140 (MQ=255)                                                                                                    ttGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCa                                   >  1:200236/1‑140 (MQ=255)                                                                                                      gAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCaaa                                 >  1:504904/1‑140 (MQ=255)                                                                                                           aaGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAAt                            <  2:200236/140‑1 (MQ=255)                                                                                                              aaaCGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAg                         >  1:760884/1‑140 (MQ=255)                                                                                                                 cGCCCCGGTGAGCGGCTGCTCGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGa                      >  2:974735/1‑140 (MQ=255)                                                                                                                                     cGTCCGGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGAGCGGCAGATTCg  >  2:543655/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       TGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCC‑GGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCACGTTTGCACGCGGCATTTACGAAATTATTAACGGGTGCGCGGGCACACTTCGCCTTCGCCAAAGAAATAAGCGATTTCGCGAGCGGCAGATTCG  >  NC_000913/2644229‑2644499

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