Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,621,471	G→A	A197T (GCC→ACC)	deoD →	purine nucleoside phosphorylase 1, nicotinamide 1‑beta‑D‑riboside phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,621,471	0	G	A	100.0%	141.7 / NA	43	A197T (GCC→ACC)	deoD	purine nucleoside phosphorylase 1, nicotinamide 1‑beta‑D‑riboside phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/22);  total (21/22)

TAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTG  >  NC_000913/4621338‑4621605                                                                                                                                      |                                                                                                                                       tAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTg                                                                                                                                   <  2:521174/139‑1 (MQ=255) tAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTg                                                                                                                                   <  1:569233/139‑1 (MQ=255)           tcCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCa                                                                                                                         <  2:80745/139‑1 (MQ=255)               cTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGt                                                                                                                     >  2:863927/1‑139 (MQ=255)                            ctctCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCGCGGCGTGGAAATGGAAGCGGCTGGTATCTAAGGCGTCGCTGCAGAATTTGGCGCGAAAAGCCTGACCATCTGCACCGTATCTGACAACAta                                                                                                        >  2:660376/1‑138 (MQ=255)                              ctcCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCg                                                                                                      <  1:1030522/139‑1 (MQ=255)                                 cGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCa                                                                                                                         >  2:732524/1‑117 (MQ=255)                                 cGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCa                                                                                                                         <  1:732524/117‑1 (MQ=255)                                     cGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGATGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCAc                                                                                               <  2:714124/139‑1 (MQ=255)                                        cGAAATGTGCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAg                                                                                            <  2:672071/139‑1 (MQ=255)                                         gAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCt                                                                                                                 <  2:540912/117‑1 (MQ=255)                                         gAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCt                                                                                                                 >  1:540912/1‑117 (MQ=255)                                                       gATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCt                                                                             >  1:1051745/1‑139 (MQ=255)                                                                         tctcGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCg                                                                                                      >  1:540021/1‑96 (MQ=255)                                                                         tctcGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCg                                                                                                      <  2:540021/96‑1 (MQ=255)                                                                          ctcGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCa                                                          <  2:1051745/139‑1 (MQ=255)                                                                               cGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGAc                                                     <  1:863927/139‑1 (MQ=255)                                                                                 tGGATATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTg                                                                                                                           <  2:722234/67‑1 (MQ=255)                                                                                 tGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATATGGCGCGAAAACCCTGACCATCTg                                                                                                                           >  1:722234/1‑67 (MQ=255)                                                                                              ggCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCa                                      >  1:359627/1‑139 (MQ=255)                                                                                                            gCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCAc                                                                                                                        <  1:122400/43‑1 (MQ=255)                                                                                                            gCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCAc                                                                                                                        >  2:122400/1‑43 (MQ=255)                                                                                                             cGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGc                                                                         >  2:1054275/1‑89 (MQ=255)                                                                                                             cGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGc                                                                         <  1:1054275/89‑1 (MQ=255)                                                                                                             cGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAAt                               <  1:888566/131‑1 (MQ=255)                                                                                                             cGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAAt                               >  2:888566/1‑131 (MQ=255)                                                                                                             cGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCACACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTctg                       >  2:88847/1‑139 (MQ=255)                                                                                                              gTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGAc                                                                                                              <  2:202820/51‑1 (MQ=255)                                                                                                              gTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGAc                                                                                                              >  1:202820/1‑51 (MQ=255)                                                                                                              gTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTctgc                      >  1:374811/1‑139 (MQ=255)                                                                                                               tCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCa                                                          >  2:1104932/1‑102 (MQ=255)                                                                                                               tCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCa                                                          <  1:1104932/102‑1 (MQ=255)                                                                                                               tCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTctgct                     >  2:17321/1‑139 (MQ=255)                                                                                                               tCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCACACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTctgct                     <  1:88847/139‑1 (MQ=255)                                                                                                                 gCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTAt                                                                                                                   <  1:174986/43‑1 (MQ=255)                                                                                                                 gCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTAt                                                                                                                   >  2:174986/1‑43 (MQ=255)                                                                                                                 gCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTgg                   <  1:17321/139‑1 (MQ=255)                                                                                                                  cTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTggg                  >  1:949802/1‑139 (MQ=255)                                                                                                                  cTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTggg                  >  2:813963/1‑139 (MQ=255)                                                                                                                    gCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCg                <  1:660376/139‑1 (MQ=255)                                                                                                                       gAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATa             <  2:914567/139‑1 (MQ=255)                                                                                                                             ggCTCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTTAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGt       >  2:711430/1‑139 (MQ=255)                                                                                                                                  gAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATtg  >  2:580653/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                       TAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTG  >  NC_000913/4621338‑4621605

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