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%	161.3 / NA	50	*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 (27/23);  total (27/23)

GGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACT  >  NC_000913/1841786‑1842030                                                                                                                |                                                                                                                                      ggCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCC‑GCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                             <  2:715863/137‑1 (MQ=255) ggCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCC‑GCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                             >  1:715863/1‑137 (MQ=255)   cgcTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                                          <  1:623774/139‑1 (MQ=255)     ctggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                                        <  2:1005389/139‑1 (MQ=255)      tggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgc                                                                                                       >  1:1110550/1‑139 (MQ=255)                  cACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGtt                                                                                                                            >  2:211303/1‑106 (MQ=255)                  cACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGtt                                                                                                                            <  1:211303/106‑1 (MQ=255)                        aaGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCt                                                                                     >  1:408123/1‑139 (MQ=255)                         agagGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCtt                                                                                    <  2:478571/139‑1 (MQ=255)                                   cAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGc                                                                          <  2:1110550/139‑1 (MQ=255)                                        tCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCtt                                                                     <  2:1251749/139‑1 (MQ=255)                                           gCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGtcat                                                                                                                     <  1:1169953/88‑1 (MQ=255)                                           gCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGtcat                                                                                                                     >  2:1169953/1‑88 (MQ=255)                                           gCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGtca                                                                                                                      <  1:1422764/87‑1 (MQ=255)                                           gCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGtca                                                                                                                      >  2:1422764/1‑87 (MQ=255)                                              ggCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGaat                                                               >  2:774103/1‑139 (MQ=255)                                                ccccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGaataa                                                             >  1:957510/1‑139 (MQ=255)                                                  ccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAtt                                                           >  1:269280/1‑139 (MQ=255)                                                  ccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAtt                                                           <  1:923377/139‑1 (MQ=255)                                                              ccATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATa                                               >  2:599570/1‑139 (MQ=255)                                                                    ttaGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                                         <  2:1318737/139‑1 (MQ=255)                                                                    ttaGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                                         <  2:408123/139‑1 (MQ=255)                                                                     taGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgcag                                        >  1:419515/1‑139 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                                  >  1:865711/1‑45 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                                  <  2:865711/45‑1 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGc                                                                          <  1:1289581/101‑1 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGc                                                                          >  2:1289581/1‑101 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCt                                                                      >  1:780869/1‑105 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCt                                                                      <  2:780869/105‑1 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCt                                                                      <  2:1227064/105‑1 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCt                                                                      >  1:1227064/1‑105 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACg                                           >  1:479971/1‑132 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACg                                           <  2:479971/132‑1 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACt                                    >  1:1206378/1‑139 (MQ=255)                                                                         ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACt                                    <  2:269280/139‑1 (MQ=255)                                                                             tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCgg                                                                                       >  1:1473980/1‑84 (MQ=255)                                                                             tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCgg                                                                                       <  2:1473980/84‑1 (MQ=255)                                                                              ttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTa                               <  2:957510/139‑1 (MQ=255)                                                                                      ttACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGc                                                                          >  2:1289639/1‑88 (MQ=255)                                                                                      ttACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGc                                                                          <  1:1289639/88‑1 (MQ=255)                                                                                       tACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                                         >  1:878783/1‑120 (MQ=255)                                                                                       tACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                                         <  2:878783/120‑1 (MQ=255)                                                                                       tACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTgc                      >  2:405463/1‑139 (MQ=255)                                                                                       tACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTgc                      >  1:812936/1‑139 (MQ=255)                                                                                       tACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTgc                      >  1:1297437/1‑139 (MQ=255)                                                                                             ccgccAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTgc                      <  2:394454/133‑1 (MQ=255)                                                                                             ccgccAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTgc                      >  1:394454/1‑133 (MQ=255)                                                                                                        gCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctga     >  1:522910/1‑139 (MQ=255)                                                                                                        gCGGAATAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGCATAATTTTGATCCTGATAAACGCAGCCCTGAGTACAGGGCTGCGCTGACTGCCCGCctga     >  2:1054230/1‑139 (MQ=255)                                                                                                           gATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctgaact  >  2:1341897/1‑139 (MQ=255)                                                                                                                |                                                                                                                                      GGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACT  >  NC_000913/1841786‑1842030

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