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

CGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTG  >  NZ_CP009273/3535253‑3535424                                                                                       |                                                                                      cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCt                                                                                    >  2:431191/1‑90 (MQ=255)   ttAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGa                                                                                  >  2:518473/1‑90 (MQ=255)    tAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAg                                                                                 <  1:194257/90‑1 (MQ=255)     aaaaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAgc                                                                                >  1:23942/1‑90 (MQ=255)           cTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTc                                                                          >  1:184493/1‑90 (MQ=255)             gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTTGTTCTGCGTGCTGGTAAAGTGTTCATTATCGTCAGCATTTCCCTGCGCGCTTTCaa                                                                        >  1:758384/1‑90 (MQ=255)             gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCaa                                                                        >  2:532794/1‑90 (MQ=255)             gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCaa                                                                        >  1:195334/1‑90 (MQ=255)                    cAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTtc                                                                 >  1:745720/1‑90 (MQ=255)                           ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                          >  1:245474/1‑90 (MQ=255)                           ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                          >  2:388896/1‑90 (MQ=255)                                 gTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGa                                                    >  2:536105/1‑90 (MQ=255)                                  tCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGaa                                                   <  2:459724/90‑1 (MQ=255)                                            ttCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGAt                                         >  1:298324/1‑90 (MQ=255)                                             tCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATa                                        >  2:658215/1‑90 (MQ=255)                                              cGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATaa                                       <  2:429911/90‑1 (MQ=255)                                                                aGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCt                     >  2:542750/1‑90 (MQ=255)                                                                   gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc                  >  1:490671/1‑90 (MQ=255)                                                                     tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAAc                                                                       <  2:419391/35‑1 (MQ=37)                                                                     tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAAc                                                                       >  1:419391/1‑35 (MQ=37)                                                                     tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt                >  2:357669/1‑90 (MQ=255)                                                                     tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt                >  1:395227/1‑90 (MQ=255)                                                                          atcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgtc           >  1:103491/1‑90 (MQ=255)                                                                           tcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCa          >  1:524760/1‑90 (MQ=255)                                                                           tcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCa          >  1:273034/1‑90 (MQ=255)                                                                           tcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCa          >  2:616916/1‑90 (MQ=255)                                                                               cAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCg      >  2:132051/1‑90 (MQ=255)                                                                               cAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCg      <  1:534575/90‑1 (MQ=255)                                                                               cAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCg      >  2:654901/1‑90 (MQ=255)                                                                                 gCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCggg    <  2:680528/90‑1 (MQ=255)                                                                                   aTTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTg  <  1:431191/90‑1 (MQ=255)                                                                                       |                                                                                      CGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTG  >  NZ_CP009273/3535253‑3535424

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