Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,621,410	T→C	I176I (ATT→ATC)	deoD →	purine nucleoside phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,621,410	0	T	C	100.0%	93.8 / NA	32	I176I (ATT→ATC)	deoD	purine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/17);  total (15/17)

TCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCT  >  NC_000913/4621167‑4621654                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      tCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGc                                                                                                                                                                                                                                                 <  1:367105/250‑1 (MQ=255)    gCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTg‑g                                                                                                                                                                                                                                            <  1:194992/251‑1 (MQ=255)        cACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAAt                                                                                                                                                                                                                                        >  1:66409/1‑251 (MQ=255)           gTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATgg                                                                                                                                                                                                                                      >  1:118651/1‑250 (MQ=255)            tAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCATCGGCGTG‑GAAATGGa                                                                                                                                                                                                                                     >  1:299225/1‑250 (MQ=255)                       gACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTAt                                                                                                                                                                                                                         >  1:69285/1‑251 (MQ=255)                         cgtcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCt                                                                                                                                                                                                                       >  1:182377/1‑251 (MQ=255)                          gtcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTa                                                                                                                                                                                                                      <  1:260269/251‑1 (MQ=255)                                            gCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTTATGCTCGCGTGTGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGagtt                                                                                                                                                                                                    <  1:291941/251‑4 (MQ=255)                                                     gATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGa                                                                                                                                                                                            >  1:80628/1‑250 (MQ=255)                                                     gATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGGATCTACGGCGTCGCTGCAGAGTTTGGCGCGa                                                                                                                                                                                            >  1:153391/1‑250 (MQ=255)                                                            acgttaactgttgtcGTTTTATAGACCATGACTTTGCCGCTATCGCGGACTTTGCCATTGTGCGTAACGCTGTTGACGCGGATAAAGCACTGGTTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCTGCGTG‑GAAATGGAAGCGGCTGGTATCTACGTCGTCGCTGCAGAGTTTGGCGCTAAAGCCCTg                                                                                                                                                                                    <  2:281037/237‑1 (MQ=255)                                                                   ccgtaTCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAAACTGTTCTCCGCTGACCTGTTCTACTCTCCTGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGAAACTACGGCGTCGGTGCGGAGTTTGGCGCGAAAACCCttgccaatc                                                                                                                                                                             >  2:291392/5‑243 (MQ=255)                                                                          cGTTTTAAAGATCATGACTTTGCCGCTATCGCTGACTTTGACATGGGGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGt                                                                                                                                                                      <  1:211020/251‑1 (MQ=255)                                                                           gTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCCTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGtg                                                                                                                                                                     <  1:227135/251‑2 (MQ=255)                                                                                                       cGCTGACGTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAg                                                                                                                                          <  1:287269/250‑1 (MQ=255)                                                                                                            aCTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACt                                                                                                                                    <  1:22409/251‑1 (MQ=255)                                                                                                                             aacgcgggagacgcggCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTAAGCGTCAGa                                                                                                                    <  1:279701/235‑1 (MQ=255)                                                                                                                                                  gCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCa                                                                                               >  1:50447/1‑250 (MQ=255)                                                                                                                                                             tGATGCTCGCGTGGGAAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTgg                                                                                   >  1:298285/1‑251 (MQ=255)                                                                                                                                                                  cgcgggTGGGTACCCTGTTCACCGCTGACCTGTTCTCCTCTCCGGACGGCGAA‑TGGGTTGTCGGGATGGAAAAATCCGGCATCCTCGGCGTG‑GAAATGGAAGCGTCTGGTATTTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGTATca                                                                              <  2:69285/246‑2 (MQ=255)                                                                                                                                                                    cgcgTGGGTATCCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATcc                                                                              <  1:267410/249‑1 (MQ=255)                                                                                                                                                                            tAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTgg                                                                    >  1:63611/1‑251 (MQ=255)                                                                                                                                                                                       cgtcTGACCTGTCCTACTCTCCGGACGGCGAACATGG‑CGACGTGATGGCAAAATACGGCATCCGCGGAGGG‑GAAATGGAAGCGGCTGGTATCTACGGTGTCGCGGCATAGGTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGAACACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTTAGACCACCTTCAACGACGTGATCAAAAGCGCACTGGAATCCTTTCTGCTGGGCGATAAag                                                           <  2:118651/246‑1 (MQ=255)                                                                                                                                                                                             aCCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATCCGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGTCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCAc                                                                                                                                                <  2:82965/158‑1 (MQ=255)                                                                                                                                                                                             aCCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCAc                                                                                                                                                >  1:82965/1‑158 (MQ=255)                                                                                                                                                                                                                          ggtcaaCGTGATGGAAAAATACGGCATCCCCGGC‑TGTGAAATGGATGCGGCGGGTATCTACGGCGTCGCTGCAGAGGTTGGCGAGAAATCCCTGACCATCGGCACCGTGTCTGACCACATCCGCACTCACGAGCATACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAATTCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGccaacgccgtgt                        <  2:103004/245‑12 (MQ=255)                                                                                                                                                                                                                          gTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCCCCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCggacgtcgtgtt                       >  1:19323/1‑240 (MQ=255)                                                                                                                                                                                                                                         aaaaTACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGAc                                                                                                                                                               <  2:342153/100‑1 (MQ=255)                                                                                                                                                                                                                                         aaaaTACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGAc                                                                                                                                                               >  1:342153/1‑100 (MQ=255)                                                                                                                                                                                                                                            aTACGGCAGCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCGGCAGAGTGTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACt     <  2:182377/251‑1 (MQ=255)                                                                                                                                                                                                                                              aCGGCATCCTCGGCGTG‑GCAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCTCGAA‑GCCCTGACCATCTGCACCGTGTCTGCCCACATCCGCACTCAAGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATAGCACTGGAATCCGTTCTGCTGTGCGATAAAGAGTAATTGTGATTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGAcggcc  >  2:360575/1‑247 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      TCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCT  >  NC_000913/4621167‑4621654

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