Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,607,596	C→T	A35A (GCG→GCA)	rsmC ←	16S rRNA m(2)G1207 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,607,596	0	C	T	100.0%	114.1 / NA	35	A35A (GCG→GCA)	rsmC	16S rRNA m(2)G1207 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/18);  total (17/18)

CAGCGTATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGCGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACGCAGCAAGACTTCACTTGCCGGGGTAAATGCAGACATGGAATGCTCCTCAATTGATACTGGCGGC  >  NC_000913/4607461‑4607728                                                                                                                                        |                                                                                                                                     cAGCGTATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCgg                                                                                                                                  <  2:656676/140‑1 (MQ=255)   gCGTATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCgg                                                                                                                                  >  2:578393/1‑138 (MQ=255)   gCGTATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCgg                                                                                                                                  <  1:578393/138‑1 (MQ=255)   gCGTATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGc                                                                                                                                >  1:401924/1‑140 (MQ=255)   gCGTATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGc                                                                                                                                >  2:511522/1‑140 (MQ=255)   gCGTATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGc                                                                                                                                >  2:731051/1‑140 (MQ=255)    cGTATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCa                                                                                                                               >  1:128245/1‑140 (MQ=255)    cGTATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCa                                                                                                                               >  1:536231/1‑140 (MQ=255)     gTATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGCCCACGGTATGTAAACGTGCGGGCAg                                                                                                                              >  1:300452/1‑140 (MQ=255)      tATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCgg                                                                                                                                  >  2:473744/1‑135 (MQ=255)      tATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCgg                                                                                                                                  <  1:473744/135‑1 (MQ=255)       aTCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGt                                                                                                                            <  1:400967/140‑1 (MQ=255)         cGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCa                                                                                                                          <  1:20763/140‑1 (MQ=255)                        cATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCg                                                                                                           <  2:100760/140‑1 (MQ=255)                          tCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCGGGCCGCGGTATCTAAATGTGCGGGCAGGTCCCCCTGTAAGTCTCCGGc                                                                                                         >  2:339936/1‑140 (MQ=255)                          tCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGAGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTTCGGGCAGGTCATCCTGTAAGTCTCCGGc                                                                                                         >  1:652992/1‑140 (MQ=255)                             gccgTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCaaa                                                                                                      <  2:401924/140‑1 (MQ=255)                                  ggCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGaa                                                                                                 >  2:648247/1‑140 (MQ=255)                                             tGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGt                                                                                      <  2:652992/140‑1 (MQ=255)                                             tGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGt                                                                                      <  1:432139/140‑1 (MQ=255)                                               aaaCGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTc                                                                                    <  2:520443/140‑1 (MQ=255)                                                 aCGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCtttcgaaca                                                                                  <  1:648247/140‑7 (MQ=255)                                                 aCGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGa                                                                                  <  2:541003/140‑1 (MQ=255)                                                  cGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGaa                                                                                 <  1:731051/140‑1 (MQ=255)                                                     gCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAatc                                                                              >  2:461865/1‑140 (MQ=255)                                                     gCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAAC‑TGCGGGCAGGTCATCCTGTAAGTCTCAGGCAAACAGAATACGGCTTTGTTCGAAatca                                                                             >  1:18045/1‑140 (MQ=255)                                                                                              gAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACGCAGCAAGACTTCACTTGCCGGGGTAAATg                                     <  2:536231/140‑1 (MQ=255)                                                                                                 ttGCTGGGTATGAGCTCGGCTGACCTCGGTATCTAAAGGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACGCAGCAAGACTTCACTTGCCGGGGTAAATGCAg                                  <  1:339936/140‑1 (MQ=255)                                                                                                  tGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACGCAGCAAGACTCCACTTGCCGGGGTAAATGCAGa                                 <  1:511522/140‑1 (MQ=255)                                                                                                   gCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAACTCATCACTGTGAAGCATCAAGACTTGACTTGCCGGGGTAAATGCAGAc                                >  2:302711/1‑140 (MQ=255)                                                                                                   gCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACGCAGCAAGACTTCACTTGCCGGGGTAAATGCAGAc                                >  2:711985/1‑140 (MQ=255)                                                                                                              gCACGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACGCAGCAAGACTTCACTTGCCGGGGTAAATGCAGACATGGAATGctc                     <  1:461865/140‑1 (MQ=255)                                                                                                                 cGGCTGGCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACACAGCAAGAATTCACCTTCCGGGGTAAATGCACACATGGAGTActcctc                  >  2:351377/1‑140 (MQ=255)                                                                                                                       gCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACGCAGCAAGACTTCACTTGCCGGGGTAAATGCAGACATGGAATGCTCCTCAATTGa            >  2:98624/1‑140 (MQ=255)                                                                                                                       gCCGCGGTATCTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACGCAGCAAGACTTCACTTGCCGGGGTAAATGCAGACATGGAATGCTCCTCAATTGa            >  2:624768/1‑140 (MQ=255)                                                                                                                                tCTAAACGTGCGGGCCGGTCATTCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACGCAGCAAGACTTCACTTGCCGGGGTAAATGCAGACATGGAATGCTCCTCAATTGATACTggcgg   <  1:351377/140‑1 (MQ=255)                                                                                                                                 cTAAACGTGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACGCAGCAAGACTTCACTTGCCGGGGTAAATGCAGACATGGAATGCTCCTCAATTGATACTggcggc  <  2:128245/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     CAGCGTATCGCAATCTGCGACGTCATCCGCCGTGGCGACCAGACTGAAACGGGCGTTATCCCCCATCTGGCGGCTTAATACCTGCCAGTGGTGGAATTGCTGGGTATGAGCACGGCTGGCCGCGGTATCTAAACGCGCGGGCAGGTCATCCTGTAAGTCTCCGGCAAACAGAATACGGCTTTGTTCGAAATCATCACTGTGACGCAGCAAGACTTCACTTGCCGGGGTAAATGCAGACATGGAATGCTCCTCAATTGATACTGGCGGC  >  NC_000913/4607461‑4607728

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