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%	286.3 / NA	76	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 (30/46);  total (30/46)

CGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAG  >  NZ_CP009273/1926927‑1927091                                                                                    |                                                                                  cgcgTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCag                                                                             <  1:2323952/90‑1 (MQ=255)     tgtTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtt                                                                         >  2:537716/1‑90 (MQ=255)     tgtTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtt                                                                         <  1:2216590/90‑1 (MQ=255)     tgtTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtt                                                                         <  2:2131681/90‑1 (MQ=255)            gTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtt                                                                         <  2:1024492/83‑1 (MQ=255)            gTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtt                                                                         >  1:1024492/1‑83 (MQ=255)               tttttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGcc                                                               <  1:2151035/90‑1 (MQ=255)                      cACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTAc                                                        >  2:1289734/1‑90 (MQ=255)                      cACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTAc                                                        >  2:1149147/1‑90 (MQ=255)                      cACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTAc                                                        <  1:207558/90‑1 (MQ=255)                              ccGGTCCAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgt                                                <  1:1452394/90‑1 (MQ=255)                                ggTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcg                                              <  2:2272410/90‑1 (MQ=255)                                ggTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcg                                              <  1:978317/90‑1 (MQ=255)                                 gTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTccatca                                                                              >  2:384054/1‑57 (MQ=255)                                 gTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTccatca                                                                              <  1:384054/57‑1 (MQ=255)                                  tACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcgtc                                            >  2:2619991/1‑90 (MQ=255)                                     aaCCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTa                                         >  2:2913912/1‑90 (MQ=255)                                      aCCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaa                                        <  2:2694049/90‑1 (MQ=255)                                      aCCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaa                                        <  2:419698/90‑1 (MQ=255)                                      aCCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaa                                        >  2:2312658/1‑90 (MQ=255)                                      aCCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaa                                        >  1:2034396/1‑90 (MQ=255)                                      aCCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaa                                        <  1:190337/90‑1 (MQ=255)                                        cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                      >  2:1578090/1‑90 (MQ=255)                                        cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                      <  2:704108/90‑1 (MQ=255)                                        cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                      <  1:1228790/90‑1 (MQ=255)                                        cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                      >  2:1785890/1‑90 (MQ=255)                                        cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaa                                        >  2:822933/1‑88 (MQ=255)                                        cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaa                                        <  1:822933/88‑1 (MQ=255)                                         gggCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaaa                                     >  1:61626/1‑90 (MQ=255)                                         gggCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaaa                                     <  1:1428020/90‑1 (MQ=255)                                         gggCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaaa                                     <  2:1185045/90‑1 (MQ=255)                                           gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTccatca                                                                              <  2:1488227/47‑1 (MQ=255)                                           gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTccatca                                                                              >  1:1488227/1‑47 (MQ=255)                                           gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                      <  1:2849305/87‑1 (MQ=255)                                           gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                      >  2:2849305/1‑87 (MQ=255)                                           gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAgg                                   >  2:2297119/1‑89 (MQ=255)                                           gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGt                                   >  2:2621285/1‑90 (MQ=255)                                            ccGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTg                                  <  1:2726502/90‑1 (MQ=255)                                              ggtggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGAtt                                                          <  2:2110664/64‑1 (MQ=255)                                              ggtggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGAtt                                                          >  1:2110664/1‑64 (MQ=255)                                              ggtggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGAt                                >  2:2148729/1‑90 (MQ=255)                                              ggtggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGAt                                >  2:40574/1‑90 (MQ=255)                                               gtggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATa                               <  2:80699/90‑1 (MQ=255)                                                tggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATAc                              >  2:2646224/1‑90 (MQ=255)                                                 ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGGGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                             <  2:187781/90‑1 (MQ=255)                                                 ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                             <  2:2792253/90‑1 (MQ=255)                                                 ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                             <  2:324053/90‑1 (MQ=255)                                                 ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                             <  2:1057425/90‑1 (MQ=255)                                                 ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                             <  2:231890/90‑1 (MQ=255)                                                 ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                             <  1:999301/90‑1 (MQ=255)                                                 ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                             <  1:730400/90‑1 (MQ=255)                                                 ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                             <  2:2406144/90‑1 (MQ=255)                                                    cAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAg                                    >  2:177950/1‑80 (MQ=255)                                                    cAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAg                                    <  1:177950/80‑1 (MQ=255)                                                      ggATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTg                        <  1:2621285/90‑1 (MQ=255)                                                           gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                   >  1:158404/1‑90 (MQ=255)                                                           gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                   <  2:158404/90‑1 (MQ=255)                                                           gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                   <  1:1863042/90‑1 (MQ=255)                                                           gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                   <  2:489740/90‑1 (MQ=255)                                                           gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                   <  2:2727833/90‑1 (MQ=255)                                                           gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                   <  1:40574/90‑1 (MQ=255)                                                           gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                   <  1:2619991/90‑1 (MQ=255)                                                                tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAATAAGTGATACAGGTTGCGCCCTGTTc              >  1:2939838/1‑90 (MQ=255)                                                                tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc              >  2:889898/1‑90 (MQ=255)                                                                tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc              <  2:2292523/90‑1 (MQ=255)                                                                tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc              <  2:2230891/90‑1 (MQ=255)                                                                tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc              >  2:159711/1‑90 (MQ=255)                                                                tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc              >  1:824727/1‑90 (MQ=255)                                                                tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc              >  1:2581692/1‑90 (MQ=255)                                                                 gCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCg             >  2:794572/1‑90 (MQ=255)                                                                      tGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCAcc        >  2:2519412/1‑90 (MQ=255)                                                                      tGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCAcc        >  2:752583/1‑90 (MQ=255)                                                                        ttttCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCgg      <  1:1108952/90‑1 (MQ=255)                                                                          ttCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGAc    <  2:585773/90‑1 (MQ=255)                                                                          ttCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGAc    <  2:1299334/90‑1 (MQ=255)                                                                           tCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACa   <  1:716977/90‑1 (MQ=255)                                                                           tCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACa   >  1:2456871/1‑90 (MQ=255)                                                                           tCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACa   <  1:1261762/90‑1 (MQ=255)                                                                            ccAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAg  <  1:823932/90‑1 (MQ=255)                                                                            ccAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAg  <  1:869284/90‑1 (MQ=255)                                                                                 tttcCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGAc    <  1:1663166/83‑1 (MQ=255)                                                                                    |                                                                                  CGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAG  >  NZ_CP009273/1926927‑1927091

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