Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,913	A→G	E104G (GAA→GGA)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,913	0	A	G	100.0%	246.4 / NA	68	E104G (GAA→GGA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/33);  total (35/33)

GCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCG  >  NC_000913/4619786‑4620047                                                                                                                                 |                                                                                                                                       gCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGATTTGCCGgtg                                                                                                                             >  2:281558/1‑139 (MQ=255) gCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGgtg                                                                                                                             >  1:724501/1‑139 (MQ=255)          cacGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTc                                                                                                                    <  1:281558/139‑1 (MQ=255)           acGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCGGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCt                                                                                                                   <  2:834928/139‑1 (MQ=255)                 gTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                             >  1:1263865/1‑139 (MQ=255)                     tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGc                                                                                                                                  >  1:256862/1‑114 (MQ=255)                     tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGc                                                                                                                                  <  2:256862/114‑1 (MQ=255)                     tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                         <  1:1423869/139‑1 (MQ=255)                        cGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGa                                                                                                      <  1:1411627/139‑1 (MQ=255)                          gTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGata                                                                                                    >  1:1583020/1‑139 (MQ=255)                            ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTc                                                                                                                           >  1:1028574/1‑114 (MQ=255)                            ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTc                                                                                                                           <  2:1028574/114‑1 (MQ=255)                            ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATAtt                                                                                                  >  1:172901/1‑139 (MQ=255)                            ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATAtt                                                                                                  >  2:942619/1‑139 (MQ=255)                                    ggcggGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGAt                                                                                          <  2:1263865/139‑1 (MQ=255)                                               aCGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAg                                                                               >  2:284776/1‑139 (MQ=255)                                                   cAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTc                                                                           >  2:718357/1‑139 (MQ=255)                                                              ttATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgct                                                                <  1:284776/139‑1 (MQ=255)                                                              ttATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgct                                                                <  1:942619/139‑1 (MQ=255)                                                              ttATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgct                                                                <  2:1583020/139‑1 (MQ=255)                                                                aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctg                                                               >  2:1585471/1‑139 (MQ=255)                                                                aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCgg                                                                                                                       <  2:1209698/82‑1 (MQ=255)                                                                aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCgg                                                                                                                       >  1:1209698/1‑82 (MQ=255)                                                                     cgcgTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctg                                                               >  2:81002/1‑133 (MQ=255)                                                                     cgcgTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctg                                                               <  1:81002/133‑1 (MQ=255)                                                                      gcgTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATa                                                        <  2:724501/139‑1 (MQ=255)                                                                       cgTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATaa                                                       >  1:116740/1‑139 (MQ=255)                                                                        gTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATaaa                                                      >  2:1337455/1‑139 (MQ=255)                                                                           cGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgct                                                                <  1:96486/126‑1 (MQ=255)                                                                           cGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgct                                                                >  2:96486/1‑126 (MQ=255)                                                                               tGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                             <  1:1491493/77‑1 (MQ=255)                                                                               tGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                             >  2:1491493/1‑77 (MQ=255)                                                                                gggCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                              <  2:172901/139‑1 (MQ=255)                                                                                     cACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGatat                                                                                                   <  1:1360435/81‑1 (MQ=255)                                                                                     cACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGatat                                                                                                   >  2:1360435/1‑81 (MQ=255)                                                                                          aaTGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAAt                                    >  1:1454299/1‑139 (MQ=255)                                                                                          aaTGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAAt                                    <  2:1454299/139‑1 (MQ=255)                                                                                             gTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCtt                                                                            <  1:607946/96‑1 (MQ=255)                                                                                             gTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCtt                                                                            >  2:607946/1‑96 (MQ=255)                                                                                                aTCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTccc                                                                         >  1:1391036/1‑96 (MQ=255)                                                                                                aTCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTccc                                                                         <  2:1391036/96‑1 (MQ=255)                                                                                                aTCC‑GGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCg                              >  1:1337521/1‑139 (MQ=255)                                                                                                       tAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTAc                        >  1:1573618/1‑139 (MQ=255)                                                                                                             tACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGt                  <  2:865964/139‑1 (MQ=255)                                                                                                                cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAg                                                                                                            >  1:1387169/1‑46 (MQ=255)                                                                                                                cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAg                                                                                                            <  2:1387169/46‑1 (MQ=255)                                                                                                                cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCt                                  >  2:1265209/1‑120 (MQ=255)                                                                                                                cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCt                                  <  1:1265209/120‑1 (MQ=255)                                                                                                                cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc               >  1:16246/1‑139 (MQ=255)                                                                                                                cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc               >  2:1256725/1‑139 (MQ=255)                                                                                                                 ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACt              <  2:116740/139‑1 (MQ=255)                                                                                                                   gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCt                                                                             >  2:122443/1‑74 (MQ=255)                                                                                                                   gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCt                                                                             <  1:122443/74‑1 (MQ=255)                                                                                                                   gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGc            <  2:16246/139‑1 (MQ=255)                                                                                                                    tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGa                                                                                                      >  2:821227/1‑48 (MQ=255)                                                                                                                    tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGa                                                                                                      <  1:821227/48‑1 (MQ=255)                                                                                                                    tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAag                                                                    >  1:1004226/1‑82 (MQ=255)                                                                                                                    tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAag                                                                    <  2:1004226/82‑1 (MQ=255)                                                                                                                     cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGt                                                                                                           >  2:1123951/1‑42 (MQ=255)                                                                                                                     cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGt                                                                                                           <  1:1123951/42‑1 (MQ=255)                                                                                                                     cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa          <  2:1573618/139‑1 (MQ=255)                                                                                                                       ggTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTtctc                                                                                              >  2:198895/1‑53 (MQ=255)                                                                                                                       ggTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTtctc                                                                                              <  1:198895/53‑1 (MQ=255)                                                                                                                        gTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGt                                                                                                           >  2:555923/1‑39 (MQ=255)                                                                                                                        gTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGt                                                                                                           <  1:555923/39‑1 (MQ=255)                                                                                                                        gTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa          >  2:61615/1‑136 (MQ=255)                                                                                                                        gTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa          <  1:61615/136‑1 (MQ=255)                                                                                                                             tGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCg  >  2:1273772/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                       GCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCC‑GGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCG  >  NC_000913/4619786‑4620047

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