Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,841,897	G→A	*136* (TAG→TAA)  T80I (ACT→ATT)	nudG → ynjH ←	CTP pyrophosphohydrolase, also hydrolyzes 2‑hydroxy‑dATP, 8‑hydroxy‑dGTP, 5‑hydroxy‑CTP, dCTP and 5‑methyl‑dCTP DUF1496 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,841,897	0	G	A	100.0%	192.5 / NA	60	*136* (TAG→TAA)  T80I (ACT→ATT)	nudG ynjH	CTP pyrophosphohydrolase, also hydrolyzes 2‑hydroxy‑dATP, 8‑hydroxy‑dGTP, 5‑hydroxy‑CTP, dCTP and 5‑methyl‑dCTP DUF1496 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/26);  total (34/26)

GGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGA  >  NC_000913/1841771‑1842027                                                                                                                               |                                                                                                                                   ggCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                               <  1:703808/132‑1 (MQ=255) ggCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                               >  2:703808/1‑132 (MQ=255)     caTGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGtc                                                                                                                    >  2:1511509/1‑139 (MQ=255)      aTGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGtca                                                                                                                   >  2:1244822/1‑139 (MQ=255)       tGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGtcat                                                                                                                  >  2:900865/1‑139 (MQ=255)          aCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCAc                                                                                                               >  2:1261331/1‑139 (MQ=255)              cAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTg                                                                                                           >  1:628712/1‑139 (MQ=255)               aGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                          <  2:85496/139‑1 (MQ=255)                ggCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGc                                                                                                         >  2:449992/1‑139 (MQ=255)                   gcTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                               >  2:155084/1‑114 (MQ=255)                   gcTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                               <  1:155084/114‑1 (MQ=255)                    ctggtctggtGCTCACCTGAAGAGGCGGTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                               <  1:1431511/113‑1 (MQ=255)                    ctggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                               >  2:1431511/1‑113 (MQ=255)                     tggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgc                                                                                                    <  1:1511509/139‑1 (MQ=255)                       gtctggtGCTCACCTGAAGAGGCGCTGCGATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                                  <  1:1554520/139‑1 (MQ=255)                       gtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGTATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                          <  1:1601980/131‑1 (MQ=255)                       gtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                          >  2:1601980/1‑131 (MQ=255)                       gtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                                  <  2:864860/139‑1 (MQ=255)                                 cACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCt                                                                                                                                >  2:1112004/1‑99 (MQ=255)                                 cACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCt                                                                                                                                <  1:1112004/99‑1 (MQ=255)                                 cACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCt                                                                                        >  2:1613873/1‑139 (MQ=255)                                 cACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCt                                                                                        <  1:900865/139‑1 (MQ=255)                                       aaGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCt                                                                                  >  1:4543/1‑139 (MQ=255)                                            gcgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                          <  1:1393359/110‑1 (MQ=255)                                            gcgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                          >  2:1393359/1‑110 (MQ=255)                                                  cAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCGTCGGCTTTAATTACCGc                                                                       <  2:1203657/139‑1 (MQ=255)                                                       tCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCtt                                                                  <  2:1343327/139‑1 (MQ=255)                                                             tggcccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGaat                                                            >  1:305080/4‑139 (MQ=255)                                                             ggCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGaat                                                            <  2:628712/139‑1 (MQ=255)                                                                 ccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAtt                                                        >  1:237873/1‑139 (MQ=255)                                                                   tgctgACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGttt                                                                                                                        <  1:488524/73‑1 (MQ=255)                                                                   tgctgACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGttt                                                                                                                        >  2:488524/1‑73 (MQ=255)                                                                    gctgACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTg                                                     <  1:981695/139‑1 (MQ=255)                                                                     ctgACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGa                                                    <  1:449992/139‑1 (MQ=255)                                                                          aTTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                                       >  2:1515737/1‑83 (MQ=255)                                                                          aTTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                                       <  1:1515737/83‑1 (MQ=255)                                                                            tCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGAt                                             <  2:1144941/139‑1 (MQ=255)                                                                                   ttaGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                                      <  2:1158393/139‑1 (MQ=255)                                                                                     agagGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                                       <  1:37940/72‑1 (MQ=255)                                                                                     agagGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                                       >  2:37940/1‑72 (MQ=255)                                                                                       agGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCAc                                  >  2:883149/1‑139 (MQ=255)                                                                                        ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCTGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGc                                                                       <  2:969141/101‑1 (MQ=255)                                                                                        ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGc                                                                       >  1:969141/1‑101 (MQ=255)                                                                                        ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACt                                 >  1:1200475/1‑139 (MQ=255)                                                                                        ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACt                                 >  1:245048/1‑139 (MQ=255)                                                                                        ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACt                                 >  1:215012/1‑139 (MQ=255)                                                                                          cGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGa                               >  2:332300/1‑139 (MQ=255)                                                                                            tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGt                             >  1:431958/1‑139 (MQ=255)                                                                                             ttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTa                            >  1:1589792/1‑139 (MQ=255)                                                                                              tATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCtt                                                                                       >  2:269803/1‑79 (MQ=255)                                                                                              tATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCtt                                                                                       <  1:269803/79‑1 (MQ=255)                                                                                                     ttACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGc                                                                                                             <  1:790957/50‑1 (MQ=255)                                                                                                     ttACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGc                                                                                                             >  2:790957/1‑50 (MQ=255)                                                                                                      tACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTgc                   >  2:558644/1‑139 (MQ=255)                                                                                                        cgcgCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACt                                 >  1:1304519/1‑123 (MQ=255)                                                                                                        cgcgCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACt                                 <  2:1304519/123‑1 (MQ=255)                                                                                                            ccgccAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGAc             >  2:386690/1‑139 (MQ=255)                                                                                                             cgccAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATc                                                  >  2:611965/1‑101 (MQ=255)                                                                                                             cgccAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATc                                                  <  1:611965/101‑1 (MQ=255)                                                                                                                       gCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctga  >  2:1012349/1‑139 (MQ=255)                                                                                                                               |                                                                                                                                   GGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGA  >  NC_000913/1841771‑1842027

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