Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,535,339	T→C	F525S (TTC→TCC)	feoB →	Fe(2+) transporter permease subunit FeoB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,535,339	0	T	C	100.0%	75.1 / NA	21	F525S (TTC→TCC)	feoB	Fe(2+) transporter permease subunit FeoB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/6);  total (15/6)

GATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAG  >  NZ_CP009273/3535265‑3535417                                                                           |                                                                               gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCaa                                                                 >  1:193199/1‑90 (MQ=255)  aTTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAAc                                                                <  1:259371/90‑1 (MQ=255)  aTTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAAc                                                                <  2:314247/90‑1 (MQ=255)          gACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTtctc                                                        >  2:357210/1‑90 (MQ=255)               ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAg                                                                          <  1:224281/67‑1 (MQ=255)               ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAg                                                                          >  2:224281/1‑67 (MQ=255)               ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                   >  1:375202/1‑90 (MQ=255)                     gTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGa                                             >  2:322347/1‑90 (MQ=255)                                ttCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGAt                                  >  2:314282/1‑90 (MQ=255)                                 tCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATa                                 >  2:99661/1‑90 (MQ=255)                                        gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCaa                          >  1:90157/1‑90 (MQ=255)                                        gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCaa                          >  1:344007/1‑90 (MQ=255)                                         cGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAAc                         >  2:180238/1‑90 (MQ=255)                                         cGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAAc                         >  2:578835/1‑90 (MQ=255)                                                       gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc           <  1:439279/90‑1 (MQ=255)                                                       gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc           >  1:288666/1‑90 (MQ=255)                                                       gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc           >  1:25566/1‑90 (MQ=255)                                                       gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc           >  2:460628/1‑90 (MQ=255)                                                         tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAAc                                                                <  2:227978/35‑1 (MQ=37)                                                         tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAAc                                                                >  1:227978/1‑35 (MQ=37)                                                                cGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAg  <  1:371313/90‑1 (MQ=255)                                                                           |                                                                               GATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAG  >  NZ_CP009273/3535265‑3535417

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