Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,495,584	A→G	K161E (AAA→GAA)	nanY →	2,7‑anhydro‑N‑acetylneuraminate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,495,584	0	A	G	100.0%	78.9 / NA	22	K161E (AAA→GAA)	nanY	2,7‑anhydro‑N‑acetylneuraminate hydratase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/15);  total (7/15)

ATCGGCGAAATATTATCATGTCATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGAAAAAGATGAAAGAACAATCTGGTGGACATCTATATCATCATATACATGAGTTAGATTGTGTTCAGCATTTACTTGGA  >  NZ_CP009273/4495515‑4495661                                                                      |                                                                              aTCGGCGAAATATTATCATGTCATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCt                                                           <  1:121220/90‑1 (MQ=255)      cGAAATATTATCATGTCATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCtggtgg                                                      <  1:352349/90‑1 (MQ=255)        aaaTATTATCATGTCATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGAc                                                    <  2:590301/90‑1 (MQ=255)         aaTATTATCATGTCATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACa                                                   <  2:234258/90‑1 (MQ=255)                 cATGTCATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATatc                                           <  2:189698/90‑1 (MQ=255)                   tGTCATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATatcat                                         <  1:153260/90‑1 (MQ=255)                    gTCATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATatcatc                                        >  2:34024/1‑90 (MQ=255)                    gTCATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATatc                                           >  1:303715/1‑87 (MQ=255)                    gTCATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATatc                                           <  2:303715/87‑1 (MQ=255)                      cATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACaa                                                              <  1:145493/66‑1 (MQ=255)                      cATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACaa                                                              >  2:145493/1‑66 (MQ=255)                       aTACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATATCATCata                                     <  1:113413/90‑1 (MQ=255)                       aTACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATATCATCata                                     >  2:101122/1‑90 (MQ=255)                            aaGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATATCATCATATACAt                                >  2:297468/1‑90 (MQ=255)                             agagAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATATCATCATATACATg                               <  2:83169/90‑1 (MQ=255)                                   aTGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATatcat                                         >  1:396647/1‑74 (MQ=255)                                   aTGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATatcat                                         <  2:396647/74‑1 (MQ=255)                                     ggCTGGGAAAACAAACAAGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATATCATCATATACATGAGTTAGAt                       <  1:34024/90‑1 (MQ=255)                                                   acaaGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGAc                                                    >  2:155352/1‑47 (MQ=255)                                                   acaaGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGAc                                                    <  1:155352/47‑1 (MQ=255)                                                   acaaGAGAGACTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATATCATCATATACATGAGTTAGATTGTGTTCAGCAttt         <  2:431834/90‑1 (MQ=255)                                                          agaCTTTCCTGGGAAAAGATGAAAGAACAATCTGGTGGACATCTATATCATCATATACATGAGTTAGATTGTGTTCAGCATTTACTTGga  <  2:194153/90‑1 (MQ=255)                                                                      |                                                                              ATCGGCGAAATATTATCATGTCATACTAAGAGAAATGGCTGGGAAAACAAACAAGAGAGACTTTCCTGGAAAAAGATGAAAGAACAATCTGGTGGACATCTATATCATCATATACATGAGTTAGATTGTGTTCAGCATTTACTTGGA  >  NZ_CP009273/4495515‑4495661

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