Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,220,044	A→G	S236P (TCG→CCG)	yohF ←	SDR family oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,220,044	0	A	G	100.0%	122.6 / NA	36	S236P (TCG→CCG)	yohF	SDR family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/15);  total (21/15)

AATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAG  >  NZ_CP009273/2219961‑2220129                                                                                    |                                                                                      aaTGAAATGTCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                                 <  1:1034390/90‑1 (MQ=255) aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                                 <  2:984033/90‑1 (MQ=255) aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                                 <  2:275641/90‑1 (MQ=255)  aTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCg                                                                                >  2:417290/1‑90 (MQ=255)   tGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCgg                                                                               <  1:306765/90‑1 (MQ=255)    gAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCggt                                                                              <  1:109261/90‑1 (MQ=255)        tatCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGgtggtg                                                                          >  1:1160980/1‑90 (MQ=255)        tatCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGgtggtg                                                                          <  2:827677/90‑1 (MQ=255)                       cTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTcc                                                           >  1:1032700/1‑90 (MQ=255)                            cTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACa                                                      >  1:505147/1‑90 (MQ=255)                                   gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                                 >  2:98869/1‑56 (MQ=255)                                   gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                                 <  1:98869/56‑1 (MQ=255)                                       tGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGgtggtg                                                                          >  1:1134210/1‑59 (MQ=255)                                       tGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGgtggtg                                                                          <  2:1134210/59‑1 (MQ=255)                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCgg                                                                               <  1:1159472/53‑1 (MQ=255)                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCgg                                                                               <  1:168465/53‑1 (MQ=255)                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCgg                                                                               >  2:168465/1‑53 (MQ=255)                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCgg                                                                               >  2:1159472/1‑53 (MQ=255)                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  1:401892/1‑90 (MQ=255)                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  2:688759/1‑90 (MQ=255)                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  2:600304/1‑90 (MQ=255)                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  2:592091/1‑90 (MQ=255)                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  2:1086816/1‑90 (MQ=255)                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  1:219679/1‑90 (MQ=255)                                             gtgGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                                     >  2:1058640/1‑90 (MQ=255)                                                 aTTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGc                                 <  2:1090032/90‑1 (MQ=255)                                                  ttCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCa                                >  2:237691/1‑90 (MQ=255)                                                  ttCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCa                                <  2:416194/90‑1 (MQ=255)                                                   tCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCaa                               <  1:198917/90‑1 (MQ=255)                                                             aaaGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGt                     <  2:828622/90‑1 (MQ=255)                                                                   gccATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGcc               >  1:717564/1‑90 (MQ=255)                                                                   gccATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGcc               >  1:168091/1‑90 (MQ=255)                                                                   gccATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGcc               >  1:340451/1‑90 (MQ=255)                                                                    ccATCCACTATCGACGGCTGCCCGGGGGGGTAATTTGCGGCCTCCGAACAAAGCCACACCACCAAGCTGGCAATCGCATGCGTTGCGCCa              >  1:36684/1‑90 (MQ=255)                                                                               cAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCaa   <  2:717564/90‑1 (MQ=255)                                                                                aaCGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAg  >  2:313872/1‑90 (MQ=255)                                                                                    |                                                                                      AATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAG  >  NZ_CP009273/2219961‑2220129

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