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%	200.1 / 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 (26/27);  total (26/27)

CGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCG  >  NC_000913/3873294‑3873583                                                                                                                                                  |                                                                                                                                                 aggggCAGACGGGCATGCCGTAGCCCACCGAACGGCGGATCACTTCACGATTGATATTGGGCGTAACGAGGAGCTGACGGCCCATCCTGGAGAGCGCAGCGACCTGGTCAGGTTTCAGGACCGTACCTGCGCCAATCAACGCCTTATc                                                                                                                                                >  1:54317/2‑148 (MQ=255)  ggggCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCg                                                                                                                                            <  2:49912/151‑1 (MQ=255)   gggCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGt                                                                                                                                           <  2:395212/151‑1 (MQ=255)      cAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACg                                                                                                                                        <  2:59615/151‑1 (MQ=255)      cAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACg                                                                                                                                        >  1:301623/1‑151 (MQ=255)                 gCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATg                                                                                                                             <  2:497789/151‑1 (MQ=255)                  ccGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATTAGCTGACATCCCATCCTGGCGAGCGCATCGACCTCTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATgg                                                                                                                            <  1:53621/151‑1 (MQ=255)                     tAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCgg                                                                                                                         >  2:458005/1‑151 (MQ=255)                      aGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCggg                                                                                                                        <  1:344027/151‑1 (MQ=255)                      aGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCgg                                                                                                                         <  2:453097/150‑1 (MQ=255)                      aGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCgg                                                                                                                         <  1:545513/150‑1 (MQ=255)                          cACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATg                                                                                                                    <  2:119360/151‑1 (MQ=255)                            ccGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCt                                                                                                                  <  2:19558/151‑1 (MQ=255)                              gCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCttt                                                                                                                >  2:543187/1‑151 (MQ=255)                               cACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTg                                                                                                               <  1:473877/151‑1 (MQ=255)                                      gATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAATGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTcccc                                                                                                         >  1:1818/1‑149 (MQ=255)                                         cacttcactATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATtgt                                                                                                     <  2:207549/151‑1 (MQ=255)                                         cacttcactATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATtg                                                                                                      <  1:705785/150‑1 (MQ=255)                                          acttcactATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATtgtg                                                                                                    <  2:560367/151‑1 (MQ=255)                                                  aTGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAAtt                                                                                             >  2:468325/1‑150 (MQ=255)                                                    ggATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTc                                                                                            >  2:637369/1‑149 (MQ=255)                                                      atatTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCg                                                                                        <  2:219499/151‑1 (MQ=255)                                                        atTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCggg                                                                                      >  2:238300/1‑151 (MQ=255)                                                        atTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCggg                                                                                      <  2:166123/151‑1 (MQ=255)                                                                 aaCGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCaaga                                                                             <  2:59771/151‑3 (MQ=255)                                                                       gAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTc                                                                        >  1:424172/1‑150 (MQ=255)                                                                       gAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTc                                                                        >  2:112982/1‑150 (MQ=255)                                                                         gCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGa                                                                      >  1:378719/1‑150 (MQ=255)                                                                            gACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGGGGGGAATTCAGCGGGATTTCAACCGCGTCGAAccc                                                                  >  1:69092/1‑151 (MQ=255)                                                                                    aTCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAAGGCTTTGCCCCCATTGGGGGGAATTCAGGGGGATTTAAACCGCGGCGAACCCGGCg                                                              >  1:19327/1‑147 (MQ=255)                                                                                        tGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCAc                                                      >  2:173982/1‑151 (MQ=255)                                                                                                gcATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCa                                               >  2:594910/1‑150 (MQ=255)                                                                                                  aTCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACa                                            <  1:256855/151‑1 (MQ=255)                                                                                                    cGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACAt                                           <  1:627465/150‑1 (MQ=255)                                                                                                     gACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATgc                                         <  2:508168/151‑1 (MQ=255)                                                                                                         tGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGcc                                      >  2:309174/1‑150 (MQ=255)                                                                                                          gTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGcc                                      >  1:211613/1‑149 (MQ=255)                                                                                                          gTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCa                                     <  1:543135/150‑1 (MQ=255)                                                                                                          gTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGGCCAACCCGGCGTCAATCACCGCGCCAACATGCGcc                                      >  1:46926/1‑149 (MQ=255)                                                                                                              aGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGcc                                <  1:746856/151‑1 (MQ=255)                                                                                                              aGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGcc                                >  2:507781/1‑151 (MQ=255)                                                                                                              aGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGcc                                >  1:709515/1‑151 (MQ=255)                                                                                                               ggTTTCAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGCATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCt                               >  2:557322/1‑151 (MQ=255)                                                                                                                    cAGTACCGTACCTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCtcggc                           <  1:719585/150‑3 (MQ=255)                                                                                                                      gtaccgtaccTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCggg                        <  2:386282/151‑1 (MQ=255)                                                                                                                        accgtaccTGCGCCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGCCGGGc                       >  1:415596/1‑150 (MQ=255)                                                                                                                                   gcCAATCGACGTCTTATCGCCGTACGCATCAACGATGGCGGGAGTGCTTTGCGCCCATTGTGGGCACTTCAGCGGGATTTCACCCGCGTCGAACCCGGCGTCAGGCACCGCGCCAACATGCGGCAGCGCCTCGTCGGGCGTACTACcgcg            <  1:723040/150‑1 (MQ=255)                                                                                                                                   gcCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCgcgc           >  1:360350/1‑151 (MQ=255)                                                                                                                                   gcCAATCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCgcgc           >  1:214475/1‑151 (MQ=255)                                                                                                                                      aaTCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCaaa        <  1:336536/151‑1 (MQ=255)                                                                                                                                       aTCAACGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCaaa        >  1:20134/1‑150 (MQ=255)                                                                                                                                            cGCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCg  <  2:427115/151‑1 (MQ=255)                                                                                                                                             gCCTTATCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCg  >  1:710268/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCG  >  NC_000913/3873294‑3873583

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