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%	119.8 / NA	35	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 (16/19);  total (16/19)

ACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAA  >  NC_000913/4621340‑4621596                                                                                                                                    |                                                                                                                              aCCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGAc                                                                                                                        <  2:796377/139‑1 (MQ=255)        ctcCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATc                                                                                                                    >  2:351133/1‑136 (MQ=255)        ctcCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATc                                                                                                                    <  1:351133/136‑1 (MQ=255)         tcCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCa                                                                                                                <  1:109704/139‑1 (MQ=255)           cGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCAcc                                                                                                              >  1:805732/1‑139 (MQ=255)                aCCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTAt                                                                                                          <  2:269631/138‑1 (MQ=255)                aCCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTAt                                                                                                          >  1:269631/1‑138 (MQ=255)                                 gACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTc                                                                                        <  2:805732/139‑1 (MQ=255)                                    ggCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACg                                                                                     >  2:837673/1‑139 (MQ=255)                                           tGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGAc                                                                              <  2:58888/139‑1 (MQ=255)                                                       tGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGa                                                                  >  2:352431/1‑139 (MQ=255)                                                               tACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGa                                                          <  1:352431/139‑1 (MQ=255)                                                                 cGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGa                                                          >  1:36281/1‑137 (MQ=255)                                                                 cGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGa                                                          <  2:36281/137‑1 (MQ=255)                                                                 cGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACt                                                        >  1:818911/1‑139 (MQ=255)                                                                                          gCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCg                               <  1:217918/139‑1 (MQ=255)                                                                                                          gCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCTAAATCGCACTGGAAt                      <  1:610175/132‑1 (MQ=255)                                                                                                          gCGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAAt                      >  2:610175/1‑132 (MQ=255)                                                                                                           cGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCt                                                                    <  2:51767/85‑1 (MQ=255)                                                                                                           cGTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCt                                                                    >  1:51767/1‑85 (MQ=255)                                                                                                           cGTCGCTGCAGAATTTGGCGCGAAAACCCTGAACATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTctg              >  2:786199/1‑139 (MQ=255)                                                                                                            gTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCa                                                                                                                >  2:14499/1‑40 (MQ=255)                                                                                                            gTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCa                                                                                                                <  1:14499/40‑1 (MQ=255)                                                                                                            gTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTAcc                                                     <  2:183769/99‑1 (MQ=255)                                                                                                            gTCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTAcc                                                     >  1:183769/1‑99 (MQ=255)                                                                                                             tCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCg                                                                      <  2:352652/81‑1 (MQ=255)                                                                                                             tCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCg                                                                      >  1:352652/1‑81 (MQ=255)                                                                                                             tCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTctg              <  2:533205/137‑1 (MQ=255)                                                                                                             tCGCTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTctg              >  1:533205/1‑137 (MQ=255)                                                                                                                cTGCAGAATTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTggg         <  2:818911/139‑1 (MQ=255)                                                                                                                       aTTTGGCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATaaa  <  2:446312/139‑1 (MQ=255)                                                                                                                           ggCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAg                                                                                   >  2:490851/1‑54 (MQ=255)                                                                                                                           ggCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAg                                                                                   <  1:490851/54‑1 (MQ=255)                                                                                                                           ggCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATa    >  1:912042/1‑133 (MQ=255)                                                                                                                           ggCGCGAAAACCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATa    <  2:912042/133‑1 (MQ=255)                                                                                                                                    |                                                                                                                              ACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAA  >  NC_000913/4621340‑4621596

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