Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,897,230	A→G	D309G (GAT→GGT)	manX →	PTS mannose transporter subunit IIAB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,897,230	0	A	G	100.0%	137.1 / NA	39	D309G (GAT→GGT)	manX	PTS mannose transporter subunit IIAB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/21);  total (18/21)

GTTTCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTC  >  NZ_CP009273/1897145‑1897316                                                                                      |                                                                                       gTTTCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCg                                                                                    >  1:465830/1‑90 (MQ=255)      ggTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACt                                                                               >  1:746484/1‑90 (MQ=255)      ggTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACt                                                                               <  1:214533/90‑1 (MQ=255)      ggTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACt                                                                               >  1:451973/1‑90 (MQ=255)       gTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTg                                                                              <  1:48817/90‑1 (MQ=255)           atgaAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGaaaa                                                                          >  2:241535/1‑90 (MQ=255)           atgaAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGaaaa                                                                          >  1:111780/1‑90 (MQ=255)           atgaAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGaaaa                                                                          >  2:560961/1‑90 (MQ=255)               aaaaGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat                                                                      <  1:558537/90‑1 (MQ=255)               aaaaGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat                                                                      <  2:442673/90‑1 (MQ=255)                              gTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCaa                                                       <  1:159945/90‑1 (MQ=255)                                    gAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGa                                                 <  1:541734/90‑1 (MQ=255)                                    gAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGa                                                 <  2:465830/90‑1 (MQ=255)                                      aaatGAATGAGAGAGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATa                                               <  1:442892/87‑1 (MQ=255)                                      aaCTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATa                                               <  1:743021/90‑1 (MQ=255)                                      aaCTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATa                                               <  2:707944/90‑1 (MQ=255)                                      aaCTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATa                                               <  2:332057/90‑1 (MQ=255)                                        cTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCaa                                                       >  1:409715/1‑80 (MQ=255)                                        cTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCaa                                                       <  2:409715/80‑1 (MQ=255)                                         tGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATc                                                   >  2:97535/1‑83 (MQ=255)                                         tGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATc                                                   <  1:97535/83‑1 (MQ=255)                                          gAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTa                                           <  1:480919/90‑1 (MQ=255)                                                gcgcgGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTAt                                     <  2:451973/90‑1 (MQ=255)                                                    ggTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATtgtg                                 >  2:532196/1‑90 (MQ=255)                                                          gAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGAtt                           <  2:712959/90‑1 (MQ=255)                                                               ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCaaa                                                      >  1:122192/1‑58 (MQ=255)                                                               ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCaaa                                                      <  2:122192/58‑1 (MQ=255)                                                               ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc                      >  2:118157/1‑90 (MQ=255)                                                               ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc                      >  2:300232/1‑90 (MQ=255)                                                               ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc                      >  2:602026/1‑90 (MQ=255)                                                                     ccGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATa                                               >  2:730046/1‑59 (MQ=255)                                                                     ccGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATa                                               <  1:730046/59‑1 (MQ=255)                                                                     ccGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGtt                >  1:688494/1‑90 (MQ=255)                                                                     ccGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGtt                >  1:729931/1‑90 (MQ=255)                                                                         aaGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTca            <  1:86453/90‑1 (MQ=255)                                                                         aaGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGGGTTGATTATCACTCAGTTTTca            >  2:679179/1‑90 (MQ=255)                                                                          aGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTcac           <  1:145394/90‑1 (MQ=255)                                                                              ttCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACtt       >  2:228240/1‑90 (MQ=255)                                                                                   ccGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTc  <  2:191984/90‑1 (MQ=255)                                                                                      |                                                                                       GTTTCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTC  >  NZ_CP009273/1897145‑1897316

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