Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,917	T→C	I105I (ATT→ATC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,917	0	T	C	100.0%	73.2 / NA	24	I105I (ATT→ATC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/14);  total (10/14)

GTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCT  >  NC_000913/4619678‑4620162                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        ggcGGGGCCGACACCCTGGGCCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCTGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACTCATGGGCTCACGAAATGTCATCCGGTAAAGAGACCCCGTCTGGTCACTGGGAAATCGCCTGTGTccc                                                                                                                                                                                                                                              <  2:100666/249‑1 (MQ=255)   cGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCg                                                                                                                                                                                                                                             <  1:98255/250‑1 (MQ=255)      ggCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTc                                                                                                                                                                                                                                         <  1:113481/251‑1 (MQ=255)          gACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCGCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGtt                                                                                                                                                                                                                                     <  2:110301/251‑1 (MQ=255)                       atatCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGatat                                                                                                                                                                                                                        <  1:534/251‑1 (MQ=255)                                          cAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACa                                                                                                                                                                                                    <  1:265496/251‑1 (MQ=255)                                             aGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGGTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCt                                                                                                                                                                                                 >  1:2307/1‑251 (MQ=255)                                                    gCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGc                                                                                                                                                                                          >  1:233481/1‑250 (MQ=255)                                                    gCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTGCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCa                                                                                                                                                                                         <  1:89884/251‑1 (MQ=255)                                                      gtattaCGGTCGTAAAGGACCGCTCAATCTGCCAAATCTGACCCGTCTGTGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATGCGGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCGTGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAag                                                                                                                                                                                       <  1:33041/246‑1 (MQ=255)                                                                     aGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGGACGCATGGGGGCACGAAATGTCATTCGGTTAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGGGGGGGGAACTTTCTCCGAT‑ACAGAACACAGCTTTCCCG‑AAGAGCTGCTGGATAAAc                                                                                                                                                                        >  2:231067/1‑251 (MQ=255)                                                                                                         gCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGTAATGGGCGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTTTCCCGGTGCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAATCAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTa                                                                                                                                    <  1:138595/251‑1 (MQ=255)                                                                                                            ggCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTaa                                                                                                                                   <  2:67421/249‑1 (MQ=255)                                                                                                                        cGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTAGCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTGGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCg                                                                                                                     <  1:38084/251‑1 (MQ=255)                                                                                                                                         tCATTCCGGCGGGACTGGAAGGCAACGCTGAAGTTATCGGCGCGTCCGCATGGGCGCACGAAATGTCATCCGGTAAAGATAACCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCCCTCTTCCGGTACGGTCATTATGGAt                                                                                                    >  2:243515/1‑251 (MQ=255)                                                                                                                                                ggcggGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGAGAACGCAATGTCATCCGGTAAAGATACACCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGGGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGTGCTGCTGGCTAAACTGGTCGAACGCGCTAATCTGCCGGGTTGTCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGGTCAGCTgg                                                                                             >  2:180664/1‑251 (MQ=255)                                                                                                                                                                     tGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCg                                                                                                                                       >  1:32003/1‑188 (MQ=255)                                                                                                                                                                     tGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCg                                                                                                                                       <  2:32003/188‑1 (MQ=255)                                                                                                                                                                                          aTGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCt                                                   >  1:180728/1‑251 (MQ=255)                                                                                                                                                                                                     aaaTGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGAAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACTCGCTAATCTGCCGGGTGTCCGCGGTAACTGCCGCTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTcc                                        <  1:125120/251‑1 (MQ=255)                                                                                                                                                                                                         gTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCgtgt                                    >  1:12368/1‑251 (MQ=255)                                                                                                                                                                                                         gTCATCCGGGAACGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCCAGCCGATTTTCTATACCTCCGCTGCCTCCgtgt                                    >  2:191547/1‑251 (MQ=255)                                                                                                                                                                                                                    aaGATACCCCGTCTGGTAACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGCTGTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgcc                         >  1:88532/1‑251 (MQ=255)                                                                                                                                                                                                                                           gAAATCGCCGGTGTCCCGGTTCTTTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTTCCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGTCTCCGTGTTCCAGATTGCCTGCCACGAAGGAAcctttggcct  <  1:231067/251‑10 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        GTCGGGGCTGACACCCTGGGTCATATCGCAGAAGCTTGTGCCAAAGGCGAAGCTGATAACGGTCGTAAAGGCCCGCTCAATCTGCCAAATCTGACCCGTCTGGGGCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAC‑GAAAACAGCTT‑CCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCT  >  NC_000913/4619678‑4620162

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: ‑