Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,621,458	A→G	E192E (GAA→GAG)	deoD →	purine nucleoside phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,621,458	0	A	G	100.0%	76.6 / NA	27	E192E (GAA→GAG)	deoD	purine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/16);  total (11/16)

GATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAA  >  NC_000913/4621219‑4621695                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              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                                                                                                                                                                                                        >  1:342153/1‑100 (MQ=255)                                                                                                                                                                                     aaaaTACGGCATCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGAc                                                                                                                                                                                                        <  2:342153/100‑1 (MQ=255)                                                                                                                                                                                        aTACGGCAGCCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCGGCAGAGTGTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACt                                              <  2:182377/251‑1 (MQ=255)                                                                                                                                                                                          aCGGCATCCTCGGCGTG‑GCAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCTCGAA‑GCCCTGACCATCTGCACCGTGTCTGCCCACATCCGCACTCAAGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATAGCACTGGAATCCGTTCTGCTGTGCGATAAAGAGTAATTGTGATTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGAcggcc                                           >  2:360575/1‑247 (MQ=255)                                                                                                                                                                                                 cctcGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGc                                      <  1:291392/249‑1 (MQ=255)                                                                                                                                                                                                          g‑gAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGTCCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTa                             <  1:236685/250‑1 (MQ=255)                                                                                                                                                                                                                                  acggcggtgCTGCAGAGATTGGCTCGAAAGCCCTGGCCATCTGCTCCGTTTCCGACCACATCCGCACTCACGAGCAGACCACTGTCGCTGAGCGTCAGAACAACTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTGCGCTGAAATGCGAACGTCGTTTTAAGCCGGAGCAGGAGACTGCTTCGGCTTTTTAGTAGCCATTCAtcaatagcaactg        <  2:247101/240‑13 (MQ=255)                                                                                                                                                                                                                                      cGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGTCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGACGGAGCAGAAGACTGCTCCGGATTTTTAGTATCTATTCAtcaatagcaactgatttat  >  2:212752/1‑232 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              GATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAA‑ATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTG‑GAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAA  >  NC_000913/4621219‑4621695

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