Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,655,057	T→C	T29A (ACG→GCG)	rspA ←	bifunctional D‑altronate/D‑mannonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,655,057	0	T	C	100.0%	242.0 / NA	68	T29A (ACG→GCG)	rspA	bifunctional D‑altronate/D‑mannonate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/37);  total (31/37)

ACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGC  >  NC_000913/1654923‑1655193                                                                                                                                       |                                                                                                                                         aCCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATa                                                                                                                                      <  1:33090/139‑1 (MQ=255)   cTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATAcc                                                                                                                                    >  2:1139182/1‑139 (MQ=255)     ttATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGt                                                                                                                                  >  1:224724/1‑139 (MQ=255)      tataGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTc                                                                                                                                 >  1:951745/1‑139 (MQ=255)       ataGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTcc                                                                                                                                <  2:783070/139‑1 (MQ=255)        taGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCt                                                                                                                               <  1:595547/139‑1 (MQ=255)             aaaCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtg                                                                                                                          <  2:1194389/139‑1 (MQ=255)                                 tACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGaa                                                                                                      >  1:1180708/1‑139 (MQ=255)                                    ggTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAAtt                                                                                                   <  2:951745/139‑1 (MQ=255)                                    ggTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAAtt                                                                                                   <  2:380067/139‑1 (MQ=255)                                    ggTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAAtt                                                                                                   >  1:395504/1‑139 (MQ=255)                                      tGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTAc                                                                                                 >  2:42645/1‑139 (MQ=255)                                      tGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTAc                                                                                                 >  1:897265/1‑139 (MQ=255)                                      tGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTAc                                                                                                 >  1:784420/1‑139 (MQ=255)                                        cgcATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGc                                                                                               >  1:619781/1‑139 (MQ=255)                                         gcATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGcc                                                                                              <  2:784420/139‑1 (MQ=255)                                          cATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGccc                                                                                             <  1:1139182/139‑1 (MQ=255)                                                  ccAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAgg                                                                                     <  1:506040/139‑1 (MQ=255)                                                         gCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtggt                                                                                                                        <  2:897824/97‑1 (MQ=255)                                                         gCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtggt                                                                                                                        >  1:897824/1‑97 (MQ=255)                                                            gCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAAc                                                                           >  2:486148/1‑139 (MQ=255)                                                             cGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATAcc                                                                                                                                    <  1:1414414/81‑1 (MQ=255)                                                             cGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATAcc                                                                                                                                    >  2:1414414/1‑81 (MQ=255)                                                             cGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGAc                                                                                        <  2:667885/125‑1 (MQ=255)                                                             cGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGAc                                                                                        >  1:667885/1‑125 (MQ=255)                                                               gACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAg                                                                                                                            <  1:1252420/87‑1 (MQ=255)                                                               gACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAg                                                                                                                            >  2:1252420/1‑87 (MQ=255)                                                               gACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTc                                                                        <  2:1462913/139‑1 (MQ=255)                                                                aCAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAAtt                                                                                                   <  1:513915/111‑1 (MQ=255)                                                                aCAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCa                                                                       <  2:897265/139‑1 (MQ=255)                                                                aCAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTCAGGGTGCCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAAtt                                                                                                   >  2:513915/1‑111 (MQ=255)                                                                 cAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAg                                                                      <  1:1105062/139‑1 (MQ=255)                                                                  aaaGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACg                                                                                                <  2:420184/112‑1 (MQ=255)                                                                  aaaGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACg                                                                                                >  1:420184/1‑112 (MQ=255)                                                                         aTCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGccc                                                                                             >  2:36085/1‑108 (MQ=255)                                                                         aTCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGccc                                                                                             <  1:36085/108‑1 (MQ=255)                                                                         aTCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACg                                                              <  2:1180708/139‑1 (MQ=255)                                                                         aTCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACg                                                              <  2:395504/139‑1 (MQ=255)                                                                         aTCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACg                                                              <  2:619781/139‑1 (MQ=255)                                                                         aTCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACg                                                              <  2:1378263/139‑1 (MQ=255)                                                                              gCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCtt                                                         >  1:63809/1‑139 (MQ=255)                                                                                    aaGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGtt                                                   >  2:578151/1‑139 (MQ=255)                                                                                       agGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGa                                                <  2:315231/139‑1 (MQ=255)                                                                                            aCGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCt                                           <  1:42645/139‑1 (MQ=255)                                                                                              ggAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGcc                                                                                              <  2:789458/86‑1 (MQ=255)                                                                                              ggAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGcc                                                                                              >  1:789458/1‑86 (MQ=255)                                                                                              ggAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTg                                         >  2:895690/1‑139 (MQ=255)                                                                                               gAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGc                                        <  2:224724/139‑1 (MQ=255)                                                                                                        cGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtg                                                                                                                          <  2:971633/48‑1 (MQ=255)                                                                                                        cGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtg                                                                                                                          >  1:971633/1‑48 (MQ=255)                                                                                                                   ggTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATg                    <  1:486148/139‑1 (MQ=255)                                                                                                                    gTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGa                   >  1:1265877/1‑139 (MQ=255)                                                                                                                     tGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGAc                                                                                        >  1:900312/1‑69 (MQ=255)                                                                                                                     tGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGAc                                                                                        <  2:900312/69‑1 (MQ=255)                                                                                                                      ggCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAAt                                                                                                              >  2:1229900/1‑46 (MQ=255)                                                                                                                      ggCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAAt                                                                                                              <  1:1229900/46‑1 (MQ=255)                                                                                                                      ggCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCAt                     >  1:618565/1‑135 (MQ=255)                                                                                                                      ggCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCAt                     <  2:618565/135‑1 (MQ=255)                                                                                                                       gCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACtt                                                                         >  2:102543/1‑82 (MQ=255)                                                                                                                       gCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACtt                                                                         <  1:102543/82‑1 (MQ=255)                                                                                                                        cATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAg                                                                                      >  2:337031/1‑68 (MQ=255)                                                                                                                        cATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAg                                                                                      <  1:337031/68‑1 (MQ=255)                                                                                                                        cATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATc               >  1:1633637/1‑139 (MQ=255)                                                                                                                             ccAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGc          >  2:1624602/1‑139 (MQ=255)                                                                                                                                aGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAAc       >  2:396795/1‑139 (MQ=255)                                                                                                                                aGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAAc       <  2:1633637/139‑1 (MQ=255)                                                                                                                                 ggCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGCAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACt      <  2:993519/139‑1 (MQ=255)                                                                                                                                    ccGCAAAACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCT‑CGCAACTGAGc  <  2:708437/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         ACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGC  >  NC_000913/1654923‑1655193

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