Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,927,010	T→C	K2E (AAA→GAA)	kdgA ←	bifunctional 4‑hydroxy‑2‑oxoglutarate aldolase/2‑dehydro‑3‑deoxy‑phosphogluconate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,927,010	0	T	C	100.0%	116.1 / NA	34	K2E (AAA→GAA)	kdgA	bifunctional 4‑hydroxy‑2‑oxoglutarate aldolase/2‑dehydro‑3‑deoxy‑phosphogluconate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/19);  total (15/19)

ACCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACC  >  NZ_CP009273/1926925‑1927085                                                                                      |                                                                            aCCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTccatc                                                                         <  2:36731/90‑1 (MQ=255)    gcgTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCaga                                                                      >  2:728814/1‑90 (MQ=255)        gtTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGttt                                                                  >  2:505284/1‑90 (MQ=255)         ttCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtttt                                                                 >  1:409245/1‑90 (MQ=255)              gTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTtctctc                                                            >  1:754346/1‑90 (MQ=255)                 tttttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGcc                                                         <  2:564110/90‑1 (MQ=255)                 tttttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTG‑TTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCt                                                        <  2:334108/90‑1 (MQ=255)                                    tACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcgtc                                      >  2:121866/1‑90 (MQ=255)                                          cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtt                                                                   >  2:432496/1‑55 (MQ=255)                                          cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtt                                                                   <  1:432496/55‑1 (MQ=255)                                          cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                <  2:863983/90‑1 (MQ=255)                                          cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                <  2:730224/90‑1 (MQ=255)                                          cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                <  2:186357/90‑1 (MQ=255)                                           gggCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaaa                               >  2:1149984/1‑90 (MQ=255)                                            ggCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAg                              <  1:644383/90‑1 (MQ=255)                                             gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGt                             >  1:919655/1‑90 (MQ=255)                                                 gtggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATa                         <  2:343684/90‑1 (MQ=255)                                                 gtggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATa                         <  2:356254/90‑1 (MQ=255)                                                   ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                       <  1:543489/90‑1 (MQ=255)                                                        ggATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTg                  <  1:1138166/90‑1 (MQ=255)                                                           ttgattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTgcgc               <  2:1070783/90‑1 (MQ=255)                                                             gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  2:551014/90‑1 (MQ=255)                                                             gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  2:180758/90‑1 (MQ=255)                                                             gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  2:754346/90‑1 (MQ=255)                                                             gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  2:828861/90‑1 (MQ=255)                                                             gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  1:943933/90‑1 (MQ=255)                                                             gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             >  1:481266/1‑90 (MQ=255)                                                             gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  1:27300/90‑1 (MQ=255)                                                                  tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTTGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc        >  1:259287/1‑90 (MQ=255)                                                                  tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc        >  1:1028933/1‑90 (MQ=255)                                                                  tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc        >  1:406233/1‑90 (MQ=255)                                                                  tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc        <  1:270821/90‑1 (MQ=255)                                                                  tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc        >  2:97304/1‑90 (MQ=255)                                                                   gCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCg       >  2:51665/1‑90 (MQ=255)                                                                        tGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCAcc  >  2:569620/1‑90 (MQ=255)                                                                                      |                                                                            ACCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACC  >  NZ_CP009273/1926925‑1927085

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: ‑