Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,920,373	G→T	A51S (GCT→TCT)	rsmF →	16S rRNA m(5)C1407 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,920,373	0	G	T	100.0%	54.3 / NA	16	A51S (GCT→TCT)	rsmF	16S rRNA m(5)C1407 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (9/7);  total (9/7)

CAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTGCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATATTCAGGAAGCCAGTTCAATGTTGCCTGTTGCCGCCTTGTTTGCTGACGGTAATGCACCACAGCGGGTGATGGATGTCGCTGCCGCGCCTG  >  NC_000913/1920137‑1920607                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           cAAATAATGCCCGCTGCTCCTGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCATGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGGCGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCAGTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCaa                                                                                                                                                                                                                              <  2:56943/251‑1 (MQ=255)   aataatGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAAtt                                                                                                                                                                                                                            <  1:219652/251‑1 (MQ=255)                   ccGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTGTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATg                                                                                                                                                                                                              <  1:100265/249‑1 (MQ=255)                        gTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGAc                                                                                                                                                                                                       >  1:95588/1‑251 (MQ=255)                                   ggAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCa                                                                                                                                                                                             >  1:130906/1‑250 (MQ=255)                                                                                  tGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGGTTGAACGCGACAATGa                                                                                                                                                 >  1:122315/1‑247 (MQ=255)                                                                                              aaCACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTg                                                                                                                                 >  1:124755/1‑251 (MQ=255)                                                                                                     cGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGTCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTATGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGgtcgt                                                                                                                           >  1:34429/1‑247 (MQ=255)                                                                                                       tttatttaCCGGACGCCTTTACGACACAACTGCGCGAAGCGATGCCTTCCACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTAcc                                                                                                                        >  2:13917/1‑251 (MQ=255)                                                                                                          aTTTCCCGGACGCCTTTCTGACACAACTGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTAccgcc                                                                                                                     <  2:49394/251‑1 (MQ=255)                                                                                                           tttCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTAccgccg                                                                                                                    >  2:236010/1‑251 (MQ=255)                                                                                                                          tCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCTCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATGCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATTAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGc                                                                                                      <  1:209774/250‑1 (MQ=255)                                                                                                                                  aaaTGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTtat                                                                                             >  2:173213/1‑251 (MQ=255)                                                                                                                                                                        gatTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATATTCAGGAAGCCAGTTCAATGTTGCCTGTTGCCGCCtt                                                       <  2:183924/251‑1 (MQ=255)                                                                                                                                                                                                                  cgcgTTAATACGCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATATTCAGGAAGCCAGTTCAATGTTGCCTGTTGCCGCCTTGTTTGCTGACGGTAATGCACCACAGCGGGTGATGGATGTCGc             >  1:80668/1‑251 (MQ=255)                                                                                                                                                                                                                             gCTGAAAATCTCCGTTTCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATATTCAGGAAGCCAGTTCAATGTTGCCTGTTGCCGCCTTGTTTGCTGACGGTAATGCACCACAGCGGGTGATGGATGTCGCTGCCGCGCCTg  <  2:124755/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           CAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCGTTGCGCCGCAGCATTCGCGTTAATACGCTGAAAATCTCCGTTGCTGATTTCCTGCAATTAACCGCTCCTTATGGCTGGACGCTTACGCCAATTCCGTGGTGTGAAGAAGGTTTCTGGATTGAACGCGACAATGAAGATGCATTGCCATTGGGTAGTACCGCCGAGCATTTAAGTGGCCTGTTTTATATTCAGGAAGCCAGTTCAATGTTGCCTGTTGCCGCCTTGTTTGCTGACGGTAATGCACCACAGCGGGTGATGGATGTCGCTGCCGCGCCTG  >  NC_000913/1920137‑1920607

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