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%	185.5 / NA	49	D58D (GAC→GAT)	dgoA	2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/34);  total (15/34)

AGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGAT  >  NC_000913/3873291‑3873585                                                                                                                                                     |                                                                                                                                                   aGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATc                                                                                                                                                  <  2:276479/151‑1 (MQ=255)  gCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCg                                                                                                                                                 <  1:752712/151‑1 (MQ=255)    cGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGcc                                                                                                                                               <  2:412692/151‑1 (MQ=255)      gggCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGt                                                                                                                                             <  1:182801/151‑1 (MQ=255)                    gCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGAt                                                                                                                                <  2:29430/150‑1 (MQ=255)                       gTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCg                                                                                                                            <  1:343231/151‑1 (MQ=255)                       gTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCg                                                                                                                            <  2:352472/151‑1 (MQ=255)                            ccACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAAt                                                                                                                       <  1:247946/151‑1 (MQ=255)                                cGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCtt                                                                                                                   >  1:575424/1‑151 (MQ=255)                                 gCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCttt                                                                                                                  <  1:609716/151‑1 (MQ=255)                                     ggcggATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTc                                                                                                              <  1:164554/151‑1 (MQ=255)                                         gATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCAt                                                                                                          <  2:23983/151‑1 (MQ=255)                                          aTCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCAtt                                                                                                         >  2:375407/1‑151 (MQ=255)                                                 cactATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGa                                                                                                  <  1:784127/151‑1 (MQ=255)                                                     aTGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAAtt                                                                                               <  2:436415/150‑1 (MQ=255)                                                      tGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTc                                                                                              >  1:237087/1‑150 (MQ=255)                                                       ggATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTc                                                                                              >  2:496601/1‑149 (MQ=255)                                                       ggATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAg                                                                                            <  1:292756/151‑1 (MQ=255)                                                                gCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTc                                                                                   <  1:712804/151‑1 (MQ=255)                                                                gCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTc                                                                                   <  2:472513/151‑1 (MQ=255)                                                                    aaCGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAAcc                                                                               <  1:481196/151‑1 (MQ=255)                                                                    aaCGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAAcc                                                                               >  1:42968/1‑151 (MQ=255)                                                                    aaCGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAAcc                                                                               <  1:639103/151‑1 (MQ=255)                                                                    aaCGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAAcc                                                                               >  1:282871/1‑151 (MQ=255)                                                                     aCGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAAc                                                                                <  1:477619/149‑1 (MQ=255)                                                                      cGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACcgc                                                                             <  2:363872/151‑1 (MQ=255)                                                                        aTGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGt                                                                           <  1:91896/151‑1 (MQ=255)                                                                           aGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGa                                                                        <  2:207839/151‑1 (MQ=255)                                                                               gACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAAccc                                                                    <  2:469024/151‑1 (MQ=255)                                                                                   gCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCg                                                                <  2:336490/151‑1 (MQ=255)                                                                                    cccATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCCCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGt                                                               >  2:494055/1‑151 (MQ=255)                                                                                      cATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCa                                                             <  2:467441/151‑1 (MQ=255)                                                                                      cATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCAGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCa                                                             <  1:584629/151‑1 (MQ=255)                                                                                         ccTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATc                                                          <  1:375382/151‑1 (MQ=255)                                                                                              cGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACcgc                                                     >  1:829196/1‑151 (MQ=255)                                                                                               gAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACcgcg                                                    <  2:569522/151‑1 (MQ=255)                                                                                                    cATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAAc                                               <  2:69887/151‑1 (MQ=255)                                                                                                        gACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACa                                              >  1:728543/1‑148 (MQ=255)                                                                                                           cTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGcc                                        >  1:193181/1‑151 (MQ=255)                                                                                                                      tCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCtcgtc                             >  1:544219/1‑151 (MQ=255)                                                                                                                      tCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCtcgt                              <  2:282896/150‑1 (MQ=255)                                                                                                                        aGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCgg                           <  2:805659/151‑1 (MQ=255)                                                                                                                           accgtaccTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGc                         <  2:725436/150‑1 (MQ=255)                                                                                                                                         aaTCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCaa           >  2:236416/1‑150 (MQ=255)                                                                                                                                            cAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACAGCGTCGAAACCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATg       >  1:363721/1‑151 (MQ=255)                                                                                                                                             aaCGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATgg      <  2:622019/151‑1 (MQ=255)                                                                                                                                             aaCGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATgg      >  2:186510/1‑151 (MQ=255)                                                                                                                                              aCGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGc     <  2:142777/151‑1 (MQ=255)                                                                                                                                                 ccTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGAt  >  1:850827/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGAT  >  NC_000913/3873291‑3873585

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