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

CGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCC  >  NZ_CP009273/3535253‑3535419                                                                                       |                                                                                 cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCt                                                                               >  2:270716/1‑90 (MQ=255) cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCt                                                                               >  2:651394/1‑90 (MQ=255)   ttAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGa                                                                             >  1:738272/1‑90 (MQ=255)   ttAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGa                                                                             >  2:782622/1‑90 (MQ=255)    tAAAAGCCTGATTATCCAGACCTGGCAGGGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTTAGGATTTCCCTGAg                                                                            >  1:52549/1‑90 (MQ=255)    tAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAg                                                                            >  2:925177/1‑90 (MQ=255)    tAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAg                                                                            <  1:995632/90‑1 (MQ=255)    tAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAg                                                                            <  2:836329/90‑1 (MQ=255)           cTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTc                                                                     >  2:604440/1‑90 (MQ=255)            tGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCa                                                                    >  1:234545/1‑90 (MQ=255)              aTTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAAc                                                                  <  1:67610/90‑1 (MQ=255)               ttATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTTTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACa                                                                 <  1:352492/90‑1 (MQ=255)                 aTCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGc                                                               >  2:43829/1‑90 (MQ=255)                           ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                     >  2:813350/1‑90 (MQ=255)                           ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCAACGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                     >  1:35848/1‑90 (MQ=255)                               gCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCgg                                                 >  2:634272/1‑90 (MQ=255)                                cGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCggg                                                >  1:661794/1‑90 (MQ=255)                                 gTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGa                                               >  1:159459/1‑90 (MQ=255)                                 gTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGa                                               >  1:509432/1‑90 (MQ=255)                                  tCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGaa                                              <  1:600624/90‑1 (MQ=255)                                     gAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAt                                           <  1:106307/90‑1 (MQ=255)                                         ggCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcgtc                                       <  1:388415/90‑1 (MQ=255)                                          gCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAgcgc                                                                         >  2:177978/1‑55 (MQ=255)                                          gCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAgcgc                                                                         <  1:177978/55‑1 (MQ=255)                                           cTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGa                                     >  1:98371/1‑90 (MQ=255)                                           cTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGa                                     >  2:86884/1‑90 (MQ=255)                                                ttCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACa                                <  1:813350/90‑1 (MQ=255)                                                    gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAgc                                                                           <  1:885165/43‑1 (MQ=255)                                                    gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAgc                                                                           >  2:885165/1‑43 (MQ=255)                                                    gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCaa                            >  1:841028/1‑90 (MQ=255)                                                    gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCaa                            >  1:241714/1‑90 (MQ=255)                                                    gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCaa                            >  1:420614/1‑90 (MQ=255)                                                      gTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACg                          >  2:715478/1‑90 (MQ=255)                                                          tGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTc                      <  1:503044/90‑1 (MQ=255)                                                               aaGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGgcgc                 <  1:881590/90‑1 (MQ=255)                                                                   gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc             >  2:307743/1‑90 (MQ=255)                                                                     tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt           >  2:418552/1‑90 (MQ=255)                                                                     tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt           >  2:204542/1‑90 (MQ=255)                                                                          atcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCCTCGATAACATCAACGACTCGGCGCTGGcgaccgtc      >  1:596488/1‑90 (MQ=255)                                                                             gTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGc   >  2:852669/1‑90 (MQ=255)                                                                              tCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGcc  <  2:806805/90‑1 (MQ=255)                                                                                       |                                                                                 CGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCC  >  NZ_CP009273/3535253‑3535419

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