Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	956,500	A→G	S45P (TCT→CCT)	ycaO ←	ribosomal protein S12 methylthiotransferase accessory factor YcaO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	956,500	0	A	G	100.0%	76.7 / NA	22	S45P (TCT→CCT)	ycaO	ribosomal protein S12 methylthiotransferase accessory factor YcaO
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/6);  total (16/6)

TTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTTCTTGGTTGCGCCTTTACCGTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATA  >  NC_000913/956389‑956618                                                                                                                |                                                                                                                       ttgAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTTCTTGGTTGCGCCTTTACCGTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGAtt                                                                                                   >  1:563532/1‑133 (MQ=255)      agaCGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTTCTTGGTTGCGCCTTTACCGTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGcc                                                                                              >  2:462695/1‑133 (MQ=255)                      cACCGAGTGCAGAAGCCAGCGCCGCTTTCTTGGTTGCGCCTTTACCGTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAAt                                                                              >  1:530867/1‑133 (MQ=255)                                       aGCGCCGCTTTCTTGGTTGCGCCTTTACCGTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTg                                                             >  1:104118/1‑133 (MQ=255)                                       aGCGCCGCTTTCTTGGTTGCGCCTTTACCGTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTg                                                             >  1:399223/1‑133 (MQ=255)                                       aGCGCCGCTTTCTTGGTTGCGCCTTTACCGTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCATCAGTCGGAAAGCAGAGGTCTg                                                             >  1:258422/1‑133 (MQ=255)                                                           gcCTTTACCGTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAagcga                                         <  2:104118/133‑1 (MQ=255)                                                            ccTTTACCGTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAagcgag                                        >  2:321293/1‑133 (MQ=255)                                                               ttACCGTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAgcgagcga                                     <  2:563532/133‑1 (MQ=255)                                                                tactgTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGAt                                    >  2:332901/5‑133 (MQ=255)                                                                       ggTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGa                                 <  1:416902/129‑1 (MQ=255)                                                                       ggTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGa                                 >  2:416902/1‑129 (MQ=255)                                                                            aaCACAGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCag                        >  1:102496/1‑133 (MQ=255)                                                                                 aGTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTTCTGGAAGGGAGCGATGGAATCTTCCAGAGCgg                   >  2:257364/1‑133 (MQ=255)                                                                                  gTGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGGAGCTAGCGATGGGATCTTCCAGAGCGGc                  >  2:282495/1‑133 (MQ=255)                                                                                   tGCGCACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCGGGAAGCGTGCGATGGAATCTTCCAGAGCGGCa                 >  2:228594/1‑133 (MQ=255)                                                                                       cACTCTTTGTTGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCtt             >  1:564295/1‑133 (MQ=255)                                                                                       cACTCTTTGTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCtt             >  1:224715/1‑133 (MQ=255)                                                                                              tGTCGCGAATATGTAAAGGCCAGACGTTTGGCACGGGATTCTGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCggg      <  2:258422/133‑1 (MQ=255)                                                                                               gTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGcc         <  1:521638/129‑1 (MQ=255)                                                                                               gTCGCGAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGcc         >  2:521638/1‑129 (MQ=255)                                                                                                  gcgAATATGTACAGGCCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATa  <  2:399223/133‑1 (MQ=255)                                                                                                                |                                                                                                                       TTGAGAGACGCTCGAAATATTCACCGAGTGCAGAAGCCAGCGCCGCTTTCTTGGTTGCGCCTTTACCGTTGGTAAAACACAGTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATA  >  NC_000913/956389‑956618

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