Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,910,133	C→T	V185I (GTC→ATC)	mazG ←	nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,910,133	0	C	T	100.0%	202.5 / NA	58	V185I (GTC→ATC)	mazG	nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/33);  total (25/33)

CGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCG  >  NC_000913/2910001‑2910267                                                                                                                                     |                                                                                                                                       cgcgAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGc                                                                                                                                  >  1:1551938/1‑139 (MQ=255) cgcgAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGc                                                                                                                                  <  2:869237/139‑1 (MQ=255)        acgacgCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCcgc                                                                                                                           >  1:1715715/1‑139 (MQ=255)           acgCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCt                                                                                                                        >  1:531535/1‑139 (MQ=255)           acgCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCt                                                                                                                        <  1:1328741/139‑1 (MQ=255)                       ttcgttcgCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCAcc                                                                                                            <  1:1013467/139‑1 (MQ=255)                               cTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAt                                                                                                    <  2:1715715/139‑1 (MQ=255)                               cTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAt                                                                                                    <  2:531535/139‑1 (MQ=255)                                tttttGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCGGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAtc                                                                                                   <  1:182804/139‑1 (MQ=255)                                tttttGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAtc                                                                                                   <  2:709223/139‑1 (MQ=255)                                      cAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCg                                                                                             <  2:1551938/139‑1 (MQ=255)                                                 cTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                                                                                  <  2:1283880/139‑1 (MQ=255)                                                 cTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                                                                                  <  1:169653/139‑1 (MQ=255)                                                       tCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTAc                                                                            <  2:1082591/139‑1 (MQ=255)                                                        cGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTa                                                                                           >  2:229152/1‑123 (MQ=255)                                                        cGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATCTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTa                                                                                           <  1:229152/123‑1 (MQ=255)                                                              tAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCa                                                                     >  2:1203979/1‑139 (MQ=255)                                                                  tGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGcg                                                                 <  2:1467744/139‑1 (MQ=255)                                                                     cGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCcgcg                                                                                                                          <  1:598448/79‑1 (MQ=255)                                                                     cGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCcgcg                                                                                                                          >  2:598448/1‑79 (MQ=255)                                                                     cGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAACCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGcgtcg                                                              >  2:1652664/1‑139 (MQ=255)                                                                       aGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCgtcgtc                                                            <  2:1394496/139‑1 (MQ=255)                                                                                  ccGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGcc                                                 >  1:1687950/1‑139 (MQ=255)                                                                                      ggCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCcgcg                                                                                                                          >  1:1507826/1‑62 (MQ=255)                                                                                      ggCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCcgcg                                                                                                                          <  2:1507826/62‑1 (MQ=255)                                                                                       gCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTg                                                                                                                               <  2:970107/56‑1 (MQ=255)                                                                                       gCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTg                                                                                                                               >  1:970107/1‑56 (MQ=255)                                                                                       gCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCtt                                                                                                                       <  1:176942/64‑1 (MQ=255)                                                                                       gCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCtt                                                                                                                       >  2:176942/1‑64 (MQ=255)                                                                                       gCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCAcc                                                                                                            >  1:1564111/1‑75 (MQ=255)                                                                                       gCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCAcc                                                                                                            <  2:1564111/75‑1 (MQ=255)                                                                                       gCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAtct                                                                                                  <  1:636845/85‑1 (MQ=255)                                                                                       gCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAtct                                                                                                  >  2:636845/1‑85 (MQ=255)                                                                                       gCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGt                                            >  1:1181702/1‑139 (MQ=255)                                                                                            cagcagGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAtc                                                                                                   >  2:327688/1‑79 (MQ=255)                                                                                            cagcagGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAtc                                                                                                   <  1:327688/79‑1 (MQ=255)                                                                                             agcagGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGgcgc                                      <  1:32683/139‑1 (MQ=255)                                                                                                  gTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGCTAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGATCAAGCGTCGTCCAATCGAAGcc                                                 <  2:1131610/123‑1 (MQ=255)                                                                                                  gTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTGGTAGACTTTATCGACTACCGGATCAAGCGTCGTCCAATCGAAGcc                                                 >  1:1131610/1‑123 (MQ=255)                                                                                                    cccccATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCcgcg                                                                                                                          >  1:142149/1‑48 (MQ=255)                                                                                                    cccccATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCcgcg                                                                                                                          <  2:142149/48‑1 (MQ=255)                                                                                                     ccccATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGcc                                                 <  1:1595238/120‑1 (MQ=255)                                                                                                     ccccATTTCCTCCTCCAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGcc                                                 >  2:1595238/1‑120 (MQ=255)                                                                                                               tcctccAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCtcgtc                                                                                                       >  1:60412/1‑56 (MQ=255)                                                                                                               tcctccAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCtcgtc                                                                                                       <  2:60412/56‑1 (MQ=255)                                                                                                                cctccAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGATCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTTTGGATTTtgtg                   <  1:1034875/139‑1 (MQ=255)                                                                                                                    cAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGcgca               <  1:1577128/139‑1 (MQ=255)                                                                                                                    cAGTTTAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGcgca               <  1:1165646/139‑1 (MQ=255)                                                                                                                        ttAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGcgcacgca           >  1:929285/1‑139 (MQ=255)                                                                                                                         tAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGcgca               >  2:1653547/1‑134 (MQ=255)                                                                                                                         tAGCCTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGcgca               <  1:1653547/134‑1 (MQ=255)                                                                                                                            ccTGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAATGTTTCTGGATTTTTTGCGCACGCATTaa       >  2:1070438/1‑139 (MQ=255)                                                                                                                              tGGTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAg     <  2:293782/139‑1 (MQ=255)                                                                                                                                gTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCg                                                      >  2:139653/1‑88 (MQ=255)                                                                                                                                gTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCg                                                      <  1:139653/88‑1 (MQ=255)                                                                                                                                gTCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGcc   >  1:1038748/1‑139 (MQ=255)                                                                                                                                 tCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGTGCACGCATTAAAGCCg  >  2:1299885/1‑139 (MQ=255)                                                                                                                                 tCGATAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCg  <  2:1038748/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       CGCGAAAACGACGCTCGAATTTTTCGTTCGCTTTTTGCAATGCGATTTCTGCTTTCGTCCCTAAATGGCGAGCCAGATTAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCG  >  NC_000913/2910001‑2910267

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