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%	248.5 / NA	68	E104G (GAA→GGA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/32);  total (36/32)

GGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGG  >  NC_000913/4619785‑4620048                                                                                                                                 |                                                                                                                                        ggCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGgt                                                                                                                               >  2:37510/1‑139 (MQ=255)  gCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGgtg                                                                                                                              >  2:1610891/1‑139 (MQ=255)  gCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGgtg                                                                                                                              >  2:686730/1‑139 (MQ=255)        gCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCgg                                                                                                                        <  2:48538/139‑1 (MQ=255)                  gTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                              <  1:1189915/139‑1 (MQ=255)                    tCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGt                                                                                                            <  2:478930/139‑1 (MQ=255)                      tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                          <  2:1382487/139‑1 (MQ=255)                      tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                          >  1:1382487/1‑139 (MQ=255)                            tttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGatat                                                                                                    >  1:1030352/1‑139 (MQ=255)                              tCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATAttt                                                                                                  <  2:513072/139‑1 (MQ=255)                                     ggcggGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGt                                                                                                            >  1:1434716/1‑122 (MQ=255)                                     ggcggGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGt                                                                                                            <  2:1434716/122‑1 (MQ=255)                                            aTGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGaaa                                                                                    <  1:37510/139‑1 (MQ=255)                                                                 aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctgc                                                               <  2:562429/139‑1 (MQ=255)                                                                 aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctgc                                                               >  1:752988/1‑139 (MQ=255)                                                                    ggCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGtt                                                                                                                      >  1:840862/1‑81 (MQ=255)                                                                    ggCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGtt                                                                                                                      <  2:840862/81‑1 (MQ=255)                                                                                    cgcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                               >  1:1347691/1‑136 (MQ=255)                                                                                    cgcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                               <  2:1347691/136‑1 (MQ=255)                                                                                    cgcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcg                                            >  1:327864/1‑139 (MQ=255)                                                                                      cACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATc                                                                                          <  2:505700/91‑1 (MQ=255)                                                                                      cACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATc                                                                                          >  1:505700/1‑91 (MQ=255)                                                                                      cACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctg                                                                <  2:142031/117‑1 (MQ=255)                                                                                      cACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctg                                                                >  1:142031/1‑117 (MQ=255)                                                                                         gAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATa                                                         <  2:1620196/121‑1 (MQ=255)                                                                                         gAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATa                                                         >  1:1620196/1‑121 (MQ=255)                                                                                                 aTCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCg                               >  1:311468/1‑139 (MQ=255)                                                                                                  tCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTccc                                                                                                                          <  2:1556674/47‑1 (MQ=255)                                                                                                  tCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTccc                                                                                                                          >  1:1556674/1‑47 (MQ=255)                                                                                                       tAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                               >  2:820292/1‑117 (MQ=255)                                                                                                       tAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                               <  1:1068295/117‑1 (MQ=255)                                                                                                       tAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                               >  2:1068295/1‑117 (MQ=255)                                                                                                       tAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                               <  1:820292/117‑1 (MQ=255)                                                                                                             tACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGt                   <  2:176370/139‑1 (MQ=255)                                                                                                              aCCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATa                                                         >  2:453054/1‑100 (MQ=255)                                                                                                              aCCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATa                                                         <  1:453054/100‑1 (MQ=255)                                                                                                               ccccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGt                                                                                                                       <  2:1084029/37‑1 (MQ=255)                                                                                                               ccccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGt                                                                                                                       >  1:1084029/1‑37 (MQ=255)                                                                                                                cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGt                   >  1:1604796/1‑136 (MQ=255)                                                                                                                cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGt                   <  2:1604796/136‑1 (MQ=255)                                                                                                                 ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgggg                                                                                                        >  1:765612/1‑50 (MQ=255)                                                                                                                 ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgggg                                                                                                        <  2:765612/50‑1 (MQ=255)                                                                                                                 ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctg                                                                <  2:1340565/90‑1 (MQ=255)                                                                                                                 ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctg                                                                >  1:1340565/1‑90 (MQ=255)                                                                                                                   gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAagag                                                                   <  2:543557/85‑1 (MQ=255)                                                                                                                   gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAagag                                                                   >  1:543557/1‑85 (MQ=255)                                                                                                                   gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATATCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGc             <  1:1096068/139‑1 (MQ=255)                                                                                                                    tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGa                                                                                                       <  1:1279747/48‑1 (MQ=255)                                                                                                                    tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGa                                                                                                       >  2:1279747/1‑48 (MQ=255)                                                                                                                    tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCa                                                                                         >  1:1151149/1‑62 (MQ=255)                                                                                                                    tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCa                                                                                         <  2:1151149/62‑1 (MQ=255)                                                                                                                    tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCAGGTTACCTCGGTAACTGcc            >  2:108579/1‑139 (MQ=255)                                                                                                                     cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAg                                                                                                             <  2:792563/41‑1 (MQ=255)                                                                                                                     cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAg                                                                                                             >  1:792563/1‑41 (MQ=255)                                                                                                                     cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCtt                                                                             <  1:985901/73‑1 (MQ=255)                                                                                                                     cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCtt                                                                             >  2:985901/1‑73 (MQ=255)                                                                                                                     cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa           <  1:1610891/139‑1 (MQ=255)                                                                                                                     cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa           <  2:311468/139‑1 (MQ=255)                                                                                                                     cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa           >  1:1116285/1‑139 (MQ=255)                                                                                                                     cTGGTCACTGGGGAATTGCCAGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGcga           >  2:975380/1‑137 (MQ=255)                                                                                                                     cTGGTCACTCGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa           >  2:1376448/1‑139 (MQ=255)                                                                                                                       ggTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCAct         >  2:31742/1‑139 (MQ=255)                                                                                                                        gTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgggg                                                                                                        <  2:208144/43‑1 (MQ=255)                                                                                                                        gTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgggg                                                                                                        >  1:208144/1‑43 (MQ=255)                                                                                                                        gTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGt                            <  2:17814/119‑1 (MQ=255)                                                                                                                        gTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGt                            >  1:17814/1‑119 (MQ=255)                                                                                                                         tCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCActct       >  2:1057543/1‑139 (MQ=255)                                                                                                                             tGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAATCGCCACTCTTCCg   >  1:1510056/1‑139 (MQ=255)                                                                                                                              ggggAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCgg  <  1:31742/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                        GGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGG  >  NC_000913/4619785‑4620048

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