Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,439	G→A	D58D (GAC→GAT)	dgoA ←	2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,439	0	G	A	100.0%	201.8 / NA	53	D58D (GAC→GAT)	dgoA	2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/22);  total (31/22)

AGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGAT  >  NC_000913/3873291‑3873585                                                                                                                                                     |                                                                                                                                                   aGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATc                                                                                                                                                  >  2:370413/1‑151 (MQ=255) aGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATc                                                                                                                                                  >  2:122013/1‑151 (MQ=255)        gCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTAc                                                                                                                                           >  2:230705/1‑151 (MQ=255)          aGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACg                                                                                                                                          <  1:613302/150‑1 (MQ=255)                         aGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCgg                                                                                                                           <  2:592159/150‑1 (MQ=255)                           cccACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGaa                                                                                                                        <  2:366166/151‑1 (MQ=255)                                 gCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCttt                                                                                                                  <  1:657843/151‑1 (MQ=255)                                  cACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCttt                                                                                                                  >  1:904551/1‑150 (MQ=255)                                   aCGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGc                                                                                                                <  2:126164/151‑1 (MQ=255)                                    cggcggATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCttt                                                                                                                  >  2:80315/1‑148 (MQ=255)                                      gcggATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTcc                                                                                                             >  2:560862/1‑151 (MQ=255)                                       cggATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTccc                                                                                                            <  1:564578/151‑1 (MQ=255)                                       cggATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTccc                                                                                                            >  2:529292/1‑151 (MQ=255)                                           tcacttcactATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCAtt                                                                                                         >  2:465736/1‑150 (MQ=255)                                                tcactATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTgggg                                                                                                   >  1:820948/1‑151 (MQ=255)                                                tcactATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTggg                                                                                                    <  2:365462/150‑1 (MQ=255)                                                 cactATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTgggg                                                                                                   <  2:710724/150‑1 (MQ=255)                                                         atatTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGGATTCAGCg                                                                                          >  2:472033/1‑151 (MQ=255)                                                         atatTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTACCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCg                                                                                          >  1:714501/1‑151 (MQ=255)                                                           atTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCgg                                                                                         <  2:159900/150‑1 (MQ=255)                                                           atTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCgg                                                                                         <  1:57384/150‑1 (MQ=255)                                                             tGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGAt                                                                                      >  1:741533/1‑151 (MQ=255)                                                               ggCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGAttt                                                                                    <  2:778622/151‑1 (MQ=255)                                                                gCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTc                                                                                   <  2:506685/151‑1 (MQ=255)                                                                       gATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACcgc                                                                             >  2:751406/1‑150 (MQ=255)                                                                             cTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAAc                                                                      >  1:152912/1‑151 (MQ=255)                                                                               gACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCCACCGCGTCGAAccc                                                                    >  2:27119/1‑151 (MQ=255)                                                                                aCAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAAccc                                                                    >  2:850957/1‑150 (MQ=255)                                                                                aCAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCg                                                                   >  1:661715/1‑151 (MQ=255)                                                                                  aGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCgg                                                                  >  2:280297/1‑150 (MQ=255)                                                                                      cATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACTCATCAACGATGGCGGAAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGt                                                               >  2:607521/1‑149 (MQ=255)                                                                                         ccTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATc                                                          <  2:319033/151‑1 (MQ=255)                                                                                           tGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCAc                                                        >  2:231869/1‑151 (MQ=255)                                                                                             gCGAGCGCATCGGCCTGTTCAGCTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCAGCAACGATGGAGTGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACAGCGTCGCACCCgg                                                                  <  2:662866/139‑1 (MQ=255)                                                                                               gAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACcgcg                                                    >  1:126327/1‑151 (MQ=255)                                                                                               gAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACcgcg                                                    <  2:288656/151‑1 (MQ=255)                                                                                                 gcgcATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGGATGCTTTGCTCCCGTTGTGGGGAATTCAGCGGGTTTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGcc                                                  >  1:498127/1‑151 (MQ=255)                                                                                                       cGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATg                                            >  1:19200/1‑151 (MQ=255)                                                                                                            tGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGcc                                        >  2:562981/1‑150 (MQ=255)                                                                                                             gTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGcc                                        >  2:479858/1‑149 (MQ=255)                                                                                                             gTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGcc                                        <  1:415061/149‑1 (MQ=255)                                                                                                                       cAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCtcgtcg                            >  2:682738/1‑151 (MQ=255)                                                                                                                       cAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCtcgtc                             <  1:682341/150‑1 (MQ=255)                                                                                                                        aGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCgg                           <  1:114476/151‑1 (MQ=255)                                                                                                                          taccgtaccTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCggg                          >  2:466440/1‑150 (MQ=255)                                                                                                                           accgtaccTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCg                        <  2:594882/151‑1 (MQ=255)                                                                                                                                     cgcCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACcgcg              >  1:189835/1‑151 (MQ=255)                                                                                                                                         aaTCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCaaa          <  2:565399/151‑1 (MQ=255)                                                                                                                                             aaCGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGGCGGGCGTAATACCGCGCAAAATgg      >  2:168298/1‑151 (MQ=255)                                                                                                                                             aaCGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCCGCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATgg      >  1:22540/1‑151 (MQ=255)                                                                                                                                               cGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCg    >  2:260984/1‑151 (MQ=255)                                                                                                                                               cGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCg    <  2:879179/151‑1 (MQ=255)                                                                                                                                                 ccTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGAt  <  2:59069/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGAT  >  NC_000913/3873291‑3873585

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