Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,501,175	T→C	I34V (ATC→GTC)	fecE ←	Fe(3+) dicitrate ABC transporter ATP‑binding protein FecE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,501,175	0	T	C	100.0%	97.2 / NA	27	I34V (ATC→GTC)	fecE	Fe(3+) dicitrate ABC transporter ATP‑binding protein FecE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/11);  total (16/11)

AACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCA  >  NZ_CP009273/4501054‑4501307                                                                                                                          |                                                                                                                                     aaCTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTcc                                                                                                                     >  2:830436/1‑139 (MQ=255)  aCTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTccc                                                                                                                    <  2:320267/139‑1 (MQ=255)        gcgATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTGAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTgg                                                                                                              >  2:154291/1‑139 (MQ=255)        gcgATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTgg                                                                                                              >  1:497880/1‑139 (MQ=255)                               ttATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTag                                                                                       >  2:210393/1‑138 (MQ=255)                                 aTCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAg                                                                                      >  1:120821/1‑138 (MQ=255)                                 aTCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAg                                                                                      >  1:156651/1‑138 (MQ=255)                                    gCCGAGAAATACGGTGCCAGACTGCGGCATTGAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTAc                                                                                   <  1:154291/138‑1 (MQ=255)                                    gCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTAc                                                                                   <  2:701880/138‑1 (MQ=255)                                                  tGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGt                                                                     >  1:309087/1‑138 (MQ=255)                                                  tGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGt                                                                     >  1:304332/1‑138 (MQ=255)                                                  tGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGt                                                                     >  2:591813/1‑138 (MQ=255)                                                       gACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGa                                                               >  2:41686/1‑139 (MQ=255)                                                            cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGt                                                           >  2:814733/1‑138 (MQ=255)                                                            cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGt                                                           >  2:770844/1‑138 (MQ=255)                                                            cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGt                                                           >  2:627328/1‑138 (MQ=255)                                                                     aaGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCa                                                 <  2:309087/139‑1 (MQ=255)                                                                           cgAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGtt                                                                                                                <  1:659023/70‑1 (MQ=255)                                                                           cgAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGtt                                                                                                                >  2:659023/1‑70 (MQ=255)                                                                            gAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTa                                          <  2:857662/139‑1 (MQ=255)                                                                            gAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTa                                          <  1:338677/139‑1 (MQ=255)                                                                              aaaaCAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTaaa                                        <  1:305083/139‑1 (MQ=255)                                                                                    gTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCAtt                                  >  1:115753/1‑139 (MQ=255)                                                                                         aCAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATc                              <  1:591813/138‑1 (MQ=255)                                                                                          cAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGc                            <  1:41686/139‑1 (MQ=255)                                                                                               tCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATtct                       >  2:513727/1‑139 (MQ=255)                                                                                                            aCCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAg          <  1:770844/139‑1 (MQ=255)                                                                                                                    gaccgaccAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCa  <  2:497880/139‑1 (MQ=255)                                                                                                                    gaccgaccAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCa  <  2:439597/139‑1 (MQ=255)                                                                                                                          |                                                                                                                                     AACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCA  >  NZ_CP009273/4501054‑4501307

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