Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	859,587	G→A	A164T (GCA→ACA)	fsa →	fructose‑6‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	859,587	0	G	A	100.0%	122.3 / NA	31	A164T (GCA→ACA)	fsa	fructose‑6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/21);  total (10/21)

ATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACA  >  NZ_CP009273/859504‑859672                                                                                    |                                                                                      aTGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGa                                                                                 <  1:63776/90‑1 (MQ=255)     tCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGttt                                                                             >  1:202759/1‑90 (MQ=255)     tCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGttt                                                                             >  1:717871/1‑90 (MQ=255)             gTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAAccc                                                                     >  2:291772/1‑90 (MQ=255)                  ggCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGt                                                                >  2:654070/1‑90 (MQ=255)                    cATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCa                                                              <  1:19822/90‑1 (MQ=255)                         aGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGgcgc                                                         <  1:686770/90‑1 (MQ=255)                             tgtgACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGt                                                                               <  2:180799/64‑1 (MQ=255)                             tgtgACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGt                                                                               >  1:180799/1‑64 (MQ=255)                                gACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTg                                                  >  1:350367/1‑90 (MQ=255)                                          cacCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCa                                        <  2:643584/90‑1 (MQ=255)                                               gTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtg                                   >  1:78639/1‑90 (MQ=255)                                               gTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtg                                   >  2:50279/1‑90 (MQ=255)                                               gTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtg                                   <  1:291772/90‑1 (MQ=255)                                                 tATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtgtg                                 <  1:825940/90‑1 (MQ=255)                                                 tATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGAtgtg                                 <  2:56885/90‑1 (MQ=255)                                                   ttGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGaa                               >  2:170640/1‑90 (MQ=255)                                                   ttGAAAATGCATGCGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGaa                               <  1:50279/90‑1 (MQ=255)                                                         atgcatgcGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAAtt                         >  1:523154/1‑90 (MQ=255)                                                             atgcGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTActc                     <  2:626713/90‑1 (MQ=255)                                                             atgcGCCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTActc                     <  2:548862/90‑1 (MQ=255)                                                               gcgcCGCAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTg                   <  1:368315/90‑1 (MQ=255)                                                                     cAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCTACTg             <  1:326771/90‑1 (MQ=255)                                                                     cAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTg             <  1:750219/90‑1 (MQ=255)                                                                     cAGGCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTg             <  1:655363/90‑1 (MQ=255)                                                                        gCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGAt          <  2:406344/90‑1 (MQ=255)                                                                        gCGAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGAt          <  2:202578/90‑1 (MQ=255)                                                                          gAAAGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGAtgt        <  2:350367/90‑1 (MQ=255)                                                                             aGTGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGc     <  2:717872/90‑1 (MQ=255)                                                                               tGCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGcac   <  1:654069/90‑1 (MQ=255)                                                                                gCTGACAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGcaca  <  2:672274/90‑1 (MQ=255)                                                                                    |                                                                                      ATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAACCCCGCGTCAGGCGCTGGACTGCTTACTGGCAGGATGTGAATCAATTACTCTGCCACTGGATGTGGCACA  >  NZ_CP009273/859504‑859672

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