Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,117,071	T→C	F44S (TTC→TCC)	rpmE →	50S ribosomal protein L31

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,117,071	0	T	C	100.0%	69.2 / NA	22	F44S (TTC→TCC)	rpmE	50S ribosomal protein L31
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/11);  total (11/11)

GCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGA  >  NZ_CP009273/4116993‑4117155                                                                               |                                                                                     gCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAAc                                                                           <  1:314819/90‑1 (MQ=255) gCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAAc                                                                           <  1:394074/90‑1 (MQ=255)       aCGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTg                                                                     >  2:225231/1‑90 (MQ=255)           aaTGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGt                                                                 >  1:253152/1‑90 (MQ=255)            aTGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAAc                                                                           >  2:51215/1‑79 (MQ=255)            aTGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAAc                                                                           <  1:51215/79‑1 (MQ=255)                aaaTCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTa                                                            <  1:225231/90‑1 (MQ=255)                                 ggTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGc                                                    >  1:454211/1‑81 (MQ=255)                                 ggTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGc                                                    <  2:454211/81‑1 (MQ=255)                                       gACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTc                                     >  1:14602/1‑90 (MQ=255)                                           tGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaaca                                 >  1:498297/1‑90 (MQ=255)                                            gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                <  2:333930/90‑1 (MQ=255)                                                 tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt                           >  2:428247/1‑90 (MQ=255)                                                 tCGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt                           >  2:188807/1‑90 (MQ=255)                                                  cGACGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGttt                          >  1:9711/1‑90 (MQ=255)                                                     cGTGTGCAGCAAGTGCCACCCGTTCTCCCCTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCaa                       <  2:506338/90‑1 (MQ=255)                                                     cGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCaa                       <  1:188807/90‑1 (MQ=255)                                                     cGTGTGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCaa                       <  2:102136/90‑1 (MQ=255)                                                      gtgtGCAGCAAGTGCCACCCGTTCTCCACTGGCAAACACCGTGATGTTGCTACCGGGGGCCGCGTTGACCGCTTCAACAAGCGTTTCGAc                      >  2:395054/1‑90 (MQ=255)                                                               aaGTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgc             <  1:495969/90‑1 (MQ=255)                                                                 gTGCCACCCGTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcag           <  1:404270/90‑1 (MQ=255)                                                                          gTTCTCCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGa  >  2:7839/1‑90 (MQ=255)                                                                               |                                                                                     GCGGTAACGTAATGAAAATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGA  >  NZ_CP009273/4116993‑4117155

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