Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,913	A→G	E104G (GAA→GGA)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,913	0	A	G	100.0%	125.7 / NA	36	E104G (GAA→GGA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/16);  total (20/16)

GCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCT  >  NC_000913/4619786‑4620043                                                                                                                                |                                                                                                                                   gCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGgtg                                                                                                                         >  2:315202/1‑139 (MQ=255)    aaaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTcc                                                                                                                      <  2:393484/139‑1 (MQ=255)       gCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCgg                                                                                                                   <  2:967860/139‑1 (MQ=255)              aaGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGtt                                                                                                            <  1:553766/139‑1 (MQ=255)                   tCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTACCGGTTCTGTTTGAGt                                                                                                       <  1:1015834/139‑1 (MQ=255)                     tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                     >  2:339218/1‑139 (MQ=255)                            ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGtt                                                                                                                 <  2:1099691/120‑1 (MQ=255)                            ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGtt                                                                                                                 >  1:1099691/1‑120 (MQ=255)                                       gggAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAc                                                                                   <  1:315202/139‑1 (MQ=255)                                                 ggCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCt                                                                         <  1:217928/139‑1 (MQ=255)                                                               tATCGGCGCGTACGCATGGGCGCACGAAATGTCATGCGGTAAAGATCCCCCTTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctg                                                           <  1:685719/139‑1 (MQ=255)                                                                 tCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctgct                                                         >  2:1095509/1‑139 (MQ=255)                                                                               tGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCg                                           <  1:1095509/139‑1 (MQ=255)                                                                                   cgcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcg                                       >  1:28853/1‑139 (MQ=255)                                                                                      aCGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACgcgc                                    >  1:247484/1‑139 (MQ=255)                                                                                                  ccGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTg                             <  1:1078927/134‑1 (MQ=255)                                                                                                  ccGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTg                             >  2:1078927/1‑134 (MQ=255)                                                                                                    ggTAAATATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGt                                                                                                                  <  1:690637/47‑1 (MQ=255)                                                                                                    ggTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGt                                                                                                                  >  2:690637/1‑47 (MQ=255)                                                                                                      tAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCt                                   <  2:40885/124‑1 (MQ=255)                                                                                                      tAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCt                                   >  1:40885/1‑124 (MQ=255)                                                                                                            tACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGt              >  1:982209/1‑139 (MQ=255)                                                                                                             aCCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTa             >  2:460354/1‑139 (MQ=255)                                                                                                               cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc           >  2:248098/1‑139 (MQ=255)                                                                                                                  gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGc        >  1:92766/1‑139 (MQ=255)                                                                                                                   tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGcc       >  1:1089536/1‑139 (MQ=255)                                                                                                                   tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGcc       >  2:485721/1‑139 (MQ=255)                                                                                                                   tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGcc       <  2:28853/139‑1 (MQ=255)                                                                                                                    cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa      >  1:646806/1‑139 (MQ=255)                                                                                                                    cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa      >  2:270439/1‑139 (MQ=255)                                                                                                                    cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa      <  2:646806/139‑1 (MQ=255)                                                                                                                     tGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCAc     >  1:928811/1‑139 (MQ=255)                                                                                                                       gTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTg                                                                                                      >  1:883281/1‑40 (MQ=255)                                                                                                                       gTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTg                                                                                                      <  2:883281/40‑1 (MQ=255)                                                                                                                       gTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCActc   >  1:328305/1‑139 (MQ=255)                                                                                                                        tCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCActct  <  2:1089536/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                   GCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCT  >  NC_000913/4619786‑4620043

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