Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,621,161	C→T	G93G (GGC→GGT)	deoD →	purine nucleoside phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,621,161	0	C	T	100.0%	98.2 / NA	30	G93G (GGC→GGT)	deoD	purine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/12);  total (18/12)

CTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGT  >  NC_000913/4621037‑4621290                                                                                                                             |                                                                                                                                  cTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCt                                                                                                                           >  2:541045/1‑133 (MQ=255)  ttACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTg                                                                                                                          >  1:353181/1‑133 (MQ=255)    aCAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGcc                                                                                                                        <  2:555051/133‑1 (MQ=255)    aCAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGcc                                                                                                                        <  1:541045/133‑1 (MQ=255)          gCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGt                                                                                                                  >  1:25286/1‑133 (MQ=255)                   tttCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTgcg                                                                                                         <  2:353181/133‑1 (MQ=255)                              ggTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTa                                                                                              >  1:32239/1‑133 (MQ=255)                                              tCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCAATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTg                                                                              >  2:413481/1‑133 (MQ=255)                                                    ccTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGa                                                                        >  2:571977/1‑133 (MQ=255)                                                    ccTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGa                                                                        >  1:367072/1‑133 (MQ=255)                                                       gCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGACGTTATCGGTATGGGTGCCTGCACCGATTc                                                                     >  1:215105/1‑133 (MQ=255)                                                             tCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGt                                                               >  1:598542/1‑133 (MQ=255)                                                                      aaGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCAt                                                      <  2:598542/133‑1 (MQ=255)                                                                       aGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATc                                                     >  1:592622/1‑133 (MQ=255)                                                                       aGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGAGCCTGCACCGATTCTAAAGTTGACCGCATc                                                     >  1:161510/1‑133 (MQ=255)                                                                           cTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGtt                                                 >  2:404750/1‑133 (MQ=255)                                                                              aTCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTa                                              >  2:81022/1‑133 (MQ=255)                                                                              aTCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTa                                              <  1:404750/133‑1 (MQ=255)                                                                               tCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTaa                                             <  2:468855/133‑1 (MQ=255)                                                                                cACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGAGTGCCTGCACCGATTCTAAAGTTAAGCGCATCCGTTTTaaa                                            >  2:283586/1‑133 (MQ=255)                                                                                   caaTTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGAc                                         <  1:413481/131‑1 (MQ=255)                                                                                     aTTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCa                                       >  1:359530/1‑133 (MQ=255)                                                                                                       ttATCCGCGTAGGTTTCTGTGGTGTAGTACTGACGCACGTAGAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGTCGCTATCGc                     <  2:194283/133‑1 (MQ=255)                                                                                                        tATCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCt                    <  2:367072/133‑1 (MQ=255)                                                                                                          tCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGa                  >  2:378326/1‑133 (MQ=255)                                                                                                          tCCGCGTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGa                  <  1:56524/133‑1 (MQ=255)                                                                                                             gcgTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGTCGCTATCGCTGACtt               >  2:286697/1‑133 (MQ=255)                                                                                                              cgTGGGTTCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTc              >  2:55850/1‑133 (MQ=255)                                                                                                                    ttCCTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTACGGGTGCCTGCACCGATTCCAAAGTTAACAGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATg        <  2:359530/133‑1 (MQ=255)                                                                                                                       cTGTGGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTg     <  2:592622/133‑1 (MQ=255)                                                                                                                         gtggtgCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCg   <  1:613882/133‑1 (MQ=255)                                                                                                                          tgGTGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGTCGCTATCGCTGACTTCGACATGGTGCGt  <  2:102965/133‑1 (MQ=255)                                                                                                                             |                                                                                                                                  CTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGT  >  NC_000913/4621037‑4621290

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