Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,621,303	G→A	A141T (GCA→ACA)	deoD →	purine nucleoside phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,621,303	0	G	A	100.0%	152.4 / NA	45	A141T (GCA→ACA)	deoD	purine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/25);  total (20/25)

ACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGG  >  NC_000913/4621175‑4621432                                                                                                                                 |                                                                                                                                  aCGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCa                                                                                                                        <  2:599567/140‑1 (MQ=255) aCGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCa                                                                                                                        <  2:696540/140‑1 (MQ=255) aCGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCa                                                                                                                        <  1:216320/140‑1 (MQ=255)  cGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCAc                                                                                                                       >  2:880972/1‑140 (MQ=255)  cGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCAc                                                                                                                       >  1:565061/1‑140 (MQ=255)  cGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCAc                                                                                                                       >  1:488341/1‑140 (MQ=255)  cGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCAc                                                                                                                       >  1:295925/1‑140 (MQ=255)  cGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCAc                                                                                                                       >  2:841748/1‑140 (MQ=255)  cGTAAAACTGCGCGACGTCGTTATCCGGATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCAc                                                                                                                       >  2:123376/1‑140 (MQ=255)   gTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACt                                                                                                                      <  2:722037/140‑1 (MQ=255)   gTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACt                                                                                                                      <  2:565061/140‑1 (MQ=255)                 cgtcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTc                                                                                                        >  2:640538/1‑140 (MQ=255)                 cgtcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTc                                                                                                        >  2:215576/1‑140 (MQ=255)                  gtcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTcg                                                                                                       <  1:439985/140‑1 (MQ=255)                  gtcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTcg                                                                                                       <  1:640538/140‑1 (MQ=255)                   tcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTcgc                                                                                                      >  2:829710/1‑140 (MQ=255)                         tCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGt                                                                                                <  1:880972/140‑1 (MQ=255)                            gTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCGGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAAc                                                                                             <  2:988439/140‑1 (MQ=255)                            gTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAAc                                                                                             <  2:85693/140‑1 (MQ=255)                             tATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAAcc                                                                                            >  2:201105/1‑140 (MQ=255)                                 ggTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGtt                                                                                        <  1:841748/140‑1 (MQ=255)                                 ggTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGtt                                                                                        <  2:1156581/140‑1 (MQ=255)                                 ggTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGtt                                                                                        >  1:684667/1‑140 (MQ=255)                                    gCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTtctc                                                                                     >  1:1053346/1‑140 (MQ=255)                                      cTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCg                                                                                   <  2:295925/140‑1 (MQ=255)                                           ccGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGAc                                                                              <  1:345724/140‑1 (MQ=255)                                           ccGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGAc                                                                              <  1:1143681/140‑1 (MQ=255)                                           ccGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGAc                                                                              <  2:962095/140‑1 (MQ=255)                                                  cAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCt                                                                       <  1:802695/140‑1 (MQ=255)                                                      gTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTActc                                                                   >  1:754146/1‑140 (MQ=255)                                                           ccGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCgg                                                              <  2:561294/140‑1 (MQ=255)                                                                      ttAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATg                                                   >  1:1017399/1‑140 (MQ=255)                                                                                tGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGa                                         <  1:504254/140‑1 (MQ=255)                                                                                tGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCAGACGGCGAAATGTTCGACGTGa                                         >  2:322594/1‑140 (MQ=255)                                                                                      tGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGCaaa                                   >  1:853379/1‑140 (MQ=255)                                                                                            tATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACg                             >  2:660766/1‑140 (MQ=255)                                                                                                   gACTTCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATtct                      <  2:779213/140‑1 (MQ=255)                                                                                                      ttCGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCgg                   <  2:637815/140‑1 (MQ=255)                                                                                                        cGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCg                 >  1:39413/1‑140 (MQ=255)                                                                                                        cGACATGGTGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTGCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCg                 <  1:123376/140‑1 (MQ=255)                                                                                                                tGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATg         >  1:1114218/1‑140 (MQ=255)                                                                                                                tGCGTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATg         >  2:315142/1‑140 (MQ=255)                                                                                                                   gTAACGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGaa      <  2:684667/140‑1 (MQ=255)                                                                                                                       cGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCgg  <  2:464966/140‑1 (MQ=255)                                                                                                                       cGCAGTAGATACAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCgg  <  1:409093/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                  ACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGG  >  NC_000913/4621175‑4621432

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑