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%	59.2 / NA	18	F525S (TTC→TCC)	feoB	Fe(2+) transporter permease subunit FeoB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/7);  total (11/7)

TTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTC  >  NZ_CP009273/3535255‑3535415                                                                                     |                                                                             ttAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGa                                                                         >  2:227686/1‑90 (MQ=255)  tAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAg                                                                        <  2:361242/90‑1 (MQ=255)         cTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTc                                                                 >  2:355331/1‑90 (MQ=255)           gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCaa                                                               >  2:105755/1‑90 (MQ=255)                        tGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTg                                                  >  1:199658/1‑90 (MQ=255)                               gTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGa                                           >  1:259919/1‑90 (MQ=255)                                           tCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATa                               <  2:199658/90‑1 (MQ=255)                                           tCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATa                               <  2:34531/90‑1 (MQ=255)                                                  gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTTGATAACATCaa                        >  1:108538/1‑90 (MQ=255)                                                        tGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTc                  <  1:180069/90‑1 (MQ=255)                                                             aaGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGgcgc             <  2:195877/90‑1 (MQ=255)                                                               gTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTg           >  1:225201/1‑90 (MQ=255)                                                                 gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc         >  2:92713/1‑90 (MQ=255)                                                                 gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc         >  2:302461/1‑90 (MQ=255)                                                                   tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt       >  2:283074/1‑90 (MQ=255)                                                                   tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt       >  2:223261/1‑90 (MQ=255)                                                                       catcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgt   <  2:247700/90‑1 (MQ=255)                                                                        atcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgtc  <  1:226298/90‑1 (MQ=255)                                                                                     |                                                                             TTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTC  >  NZ_CP009273/3535255‑3535415

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