Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,932,990	A→G	E362G (GAA→GGA)	pyk →	pyruvate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,932,990	0	A	G	100.0%	63.0 / NA	19	E362G (GAA→GGA)	pyk	pyruvate kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/9);  total (10/9)

AAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCT  >  NZ_CP009273/1932864‑1933110                                                                                                                               |                                                                                                                         aaaCTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGc                                                                                                              >  2:549666/1‑139 (MQ=255)        cGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCa                                                                                                       <  2:407350/139‑1 (MQ=255)             ggCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGggaga                                                                                                                       <  1:242132/118‑1 (MQ=255)             ggCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGggaga                                                                                                                       >  2:242132/1‑118 (MQ=255)                              aCCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAaccacc                                                                                 >  2:438480/1‑139 (MQ=255)                                          aTGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAac                                                                                     <  2:502635/123‑1 (MQ=255)                                          aTGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAac                                                                                     >  1:502635/1‑123 (MQ=255)                                                            ggTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGa                                                   <  1:419162/139‑1 (MQ=255)                                                                  gAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCa                                                                                                    >  1:577287/1‑84 (MQ=255)                                                                  gAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCa                                                                                                    <  2:577287/84‑1 (MQ=255)                                                                      aaaTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAGCCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCggg                                         <  1:186587/139‑1 (MQ=255)                                                                                          tCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGt                                                                       <  1:168321/89‑1 (MQ=255)                                                                                          tCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGt                                                                       >  2:168321/1‑89 (MQ=255)                                                                                                   cGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGt                                                                                                         <  2:285314/46‑1 (MQ=255)                                                                                                   cGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGt                                                                                                         >  1:285314/1‑46 (MQ=255)                                                                                                   cGTCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTAt              >  1:250617/1‑137 (MQ=255)                                                                                                     tCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGc          >  1:63791/1‑139 (MQ=255)                                                                                                     tCTGGACGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGc          >  1:477358/1‑139 (MQ=255)                                                                                                          aCGTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTgg     <  1:207913/139‑1 (MQ=255)                                                                                                            gTTCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTc   <  1:418825/139‑1 (MQ=255)                                                                                                             ttCAGTTCGACAATGTGGGAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCt  >  1:65062/1‑139 (MQ=255)                                                                                                                               |                                                                                                                         AAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCT  >  NZ_CP009273/1932864‑1933110

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