Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,717	T→G	A163A (GCT→GCG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,717	0	T	G	100.0%	76.9 / NA	26	A163A (GCT→GCG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/13);  total (13/13)

TGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATC  >  NC_000913/4618474‑4618958                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      tgATAAACATTCCACTGGTGGCGTCGGCGATGTGCCGTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATaaa                                                                                                                                                                                                                                                 <  1:106495/250‑1 (MQ=255)   aTAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGt                                                                                                                                                                                                                                              <  1:206594/251‑1 (MQ=255)                                            aTGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCAt                                                                                                                                                                                                   <  1:268154/251‑1 (MQ=255)                                                          ggTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑cc                                                                                                                                                                                     <  2:170274/250‑1 (MQ=255)                                                          ggTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACGGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑cc                                                                                                                                                                                     <  1:126544/250‑1 (MQ=255)                                                                                         atCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGCATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGCAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGATCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGGGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGACGAAGAAAcacgcgggagg                                                                                                                                                     >  1:163296/1‑241 (MQ=255)                                                                                                     ggCCTCGGTCATACTGACGGTACGCTCGACAAACTGGAATCCCTCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCCCCGGCGGATAAACGTTTTTACGCTGACC‑GG‑GATATTAACCCAACCGGGGACCCCCTCCCGCTGATCAACCACC‑CTATTCTGGCGAAGAAACCGGCGGAAGGCCTGGAcgggcc                                                                                                                                         >  2:93857/1‑247 (MQ=255)                                                                                                                  cTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGCCGTg                                                                                                                            >  1:253679/1‑251 (MQ=255)                                                                                                                             cccGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCACCTCACTGGCCCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTCCCGCACCCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAg                                                                                                                <  1:4776/249‑1 (MQ=255)                                                                                                                                     aCTGGAATCCATCCCTGGCATAGGCATTTTCACGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTTTACGC‑GAACCGT‑GAAATTACCGCAACCGGGGACTACATCACGCTGATCA‑CCGCCTCTAATCTGGCGAAGAAACTGGCGGAAGGTCTGGACGAGCTGGGGCTGGACGGG‑GAAGTGGGTAGCGGagcgg                                                                                                        >  2:238092/1‑250 (MQ=255)                                                                                                                                          aaTCCATCCCTGGCTTCGACATTTTCCCGGATGACACCCGTTCCCGCGAAATTATTAAAGCCGTCGGCGTGGCGATTATCGGTCAGCCCAGCTCACTGGCTCCGGCGGAAAAACGTTTCTACGC‑GACCCGT‑GATATTCCCGCAACCGTGGACCCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATg                                                                                                   <  1:194403/251‑1 (MQ=255)                                                                                                                                                    tGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTCAGC‑GAACCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGGTCA‑CCGCCTCTCTTCTGGCGAAGAAACTGGCGGGAGGTCTGGACGCGCTGGTGATGGACGGG‑AAAGTGGGAAGCGGCGCGTTTATGGCGCCCTACg                                                                                         >  1:106300/1‑251 (MQ=255)                                                                                                                                                        ttCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGCTCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTGTATgacggccgccgaacc                                                                                     >  2:224464/1‑237 (MQ=255)                                                                                                                                                           gactttttCCCGGATGACAACCTTTCCCGCGAAATTATTAAAGCGGGCGTCGTGGCGATTACCGGTAAGCCCACCTCACTGGCCCCGGCGGATAATCTTTTCTCCGC‑GACCCTT‑GATATTACCGCACCCGTGGACTCCATCCCGCTTATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAActctc                                                                                  <  1:235125/247‑1 (MQ=255)                                                                                                                                                                          gACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAAAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGAGGGTAGCGGCGCGTTTATGCCGACCTACGAACGCTCTGAAGCCCTTGCCgg                                                                   >  1:200242/1‑250 (MQ=255)                                                                                                                                                                                 gTTTCCGCGAACTTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGCCCAGCTAACTGGCTCCGGAGGATAAACGTTTCTACGC‑GACCCGT‑GATATTAACGAAACCGTGGACTCAATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGAG‑AAAGGAGGTAGCGGCGCGTTTAGGCCGACCTACGAACTCTCTAGAGGCCTTGCCGATGCGAtt                                                            >  2:41727/1‑251 (MQ=255)                                                                                                                                                                                   ttCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTAGCGTTCAGACCAGCCCACGGGCTCCGGCGGATAACCTTTTCTCCGC‑G‑CCCGTTGATATTACCGCAACCGTGGACCCCACCCCGCTGACCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGAttg                                                           <  1:133157/250‑1 (MQ=255)                                                                                                                                                                                           aaTTATTAAAGCCGTGGTGGTGGCGATTATCGGTCAGACCAGCTCGCTGGCTCCGGCGTATAAACTTTTCGACGC‑GCCCGGA‑GATATGACCCCCACCGTGGCCTCCATCCCTCTGACCG‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACGCTCTGAAGCCCTTGCCGAAGCGCTTGTTGGCGTgg                                                  <  2:56203/251‑1 (MQ=255)                                                                                                                                                                                           aaTTATTAAACACGTCGGCGTGGCGATTATCGGTCAGACCAGATCACTGGCTCCGGCGGATAAACGTTTCTAAGC‑GACCCGT‑GATATTCCCGCAACCGTGGACTCCATACCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGACGGTCGGGACGCGCTGGGGATGGACGTGCAAAG‑GGGAAGCGGCGCGTTTATGCCGACCTACGAACTCTATGAAGCCCTTGCCGCAGCGATGGTTGGCGTgg                                                  >  2:104550/1‑251 (MQ=255)                                                                                                                                                                                               attaAAGACGGCGGCGTGGCGATTATCGGTCAGACCAGCCCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTaa                                              <  1:51356/251‑1 (MQ=255)                                                                                                                                                                                                  aaaGCCGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATCCCCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACg                                            <  1:191348/250‑1 (MQ=255)                                                                                                                                                                                                  aaaGACGTCGGCGTGGCGCTTATCGGTCAGACCCGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGGAGCGGCGCGTTTATGCCGACCTTCGAACTCTCTGAAGCCCTTGCCGAAGCGGTTGTTGGCGTGGCTAACgg                                           >  1:208980/1‑251 (MQ=255)                                                                                                                                                                                                  aaaGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGG‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGGCGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGGCCTACGAACTCTCTGCAGCCCTTGCCGCAGCGCTTGTTGGCGTGGCTAACgg                                           >  2:29872/1‑251 (MQ=255)                                                                                                                                                                                                                         cGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGCAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAACCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCgcgc                    >  1:265862/1‑251 (MQ=255)                                                                                                                                                                                                                                      tCACTGGCTCCGGCGGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGCTTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACAt       >  1:167958/1‑251 (MQ=255)                                                                                                                                                                                                                                           ggCTCCGGCGGAAAAACGTTTCACCGC‑GACCTGG‑GATATTCCCGCAACCGTGGACCCTCCCCCGCATATCA‑CCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATc  <  1:224464/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      TGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGC‑GACCCGT‑GATATTACCGCAACCGTGGACTCCATCCCGCTGATCA‑CCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTG‑AAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATC  >  NC_000913/4618474‑4618958

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