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%	169.5 / NA	48	E104G (GAA→GGA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/25);  total (23/25)

CGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCA  >  NC_000913/4619787‑4620039                                                                                                                               |                                                                                                                               caaaaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGgtgt                                                                                                                    >  2:171804/3‑139 (MQ=255)        acacGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGtt                                                                                                             <  1:1192280/139‑1 (MQ=255)             aaGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGtt                                                                                                        <  1:266466/139‑1 (MQ=255)                    tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                 >  2:873139/1‑139 (MQ=255)                    tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                 >  1:1176894/1‑139 (MQ=255)                    tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                 <  2:580397/139‑1 (MQ=255)                           ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGtt                                                                                                             <  1:659788/120‑1 (MQ=255)                           ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGtt                                                                                                             >  2:659788/1‑120 (MQ=255)                           ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                     <  1:778543/128‑1 (MQ=255)                           ttCATTCCGGCGGGAATGGACGGCAACGCAGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                     >  2:778543/1‑128 (MQ=255)                               ttCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTtctc                                                                                      <  2:719349/139‑1 (MQ=255)                                    gcggGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATc                                                                                 <  1:171804/139‑1 (MQ=255)                                                  cAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                 >  2:625636/1‑109 (MQ=255)                                                  cAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                 <  1:625636/109‑1 (MQ=255)                                                         gAAGTTATCGGAGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAag                                                            <  1:873139/139‑1 (MQ=255)                                                                  ggCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTCTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTgg                                                   >  1:717614/1‑139 (MQ=255)                                                                    cgcgTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATGCCCCGTCTGTTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAagag                                                          >  2:883961/1‑130 (MQ=255)                                                                    cgcgTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATGCCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAagag                                                          <  1:883961/130‑1 (MQ=255)                                                                               gggCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                      <  2:1176894/139‑1 (MQ=255)                                                                                  cgcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAag                                                            >  2:1165182/1‑114 (MQ=255)                                                                                  cgcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAag                                                            <  1:1165182/114‑1 (MQ=255)                                                                                  cgcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTATTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcg                                   >  1:839622/1‑139 (MQ=255)                                                                                    cACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcgcg                                 >  2:1049507/1‑139 (MQ=255)                                                                                       gAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTa                              >  1:1216739/1‑139 (MQ=255)                                                                                               aTCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCgg                                                                                                               <  2:644764/50‑1 (MQ=255)                                                                                               aTCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCgg                                                                                                               >  1:644764/1‑50 (MQ=255)                                                                                                  cGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGt                   <  2:1216739/139‑1 (MQ=255)                                                                                                   ggTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGgtg                                                                                                                     >  1:63400/1‑40 (MQ=255)                                                                                                   ggTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGgtg                                                                                                                     <  2:63400/40‑1 (MQ=255)                                                                                                    gTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACg                                                                              <  1:425358/78‑1 (MQ=255)                                                                                                    gTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACg                                                                              >  2:425358/1‑78 (MQ=255)                                                                                                     tAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTAc                <  1:1082490/139‑1 (MQ=255)                                                                                                            aCCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctg                                                       <  1:649875/93‑1 (MQ=255)                                                                                                            aCCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctg                                                       >  2:649875/1‑93 (MQ=255)                                                                                                              cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc       >  1:33045/1‑139 (MQ=255)                                                                                                              cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc       >  1:156566/1‑139 (MQ=255)                                                                                                               ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATc                                                                                 <  2:937972/64‑1 (MQ=255)                                                                                                               ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATc                                                                                 >  1:937972/1‑64 (MQ=255)                                                                                                                 gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGc    >  1:582886/1‑139 (MQ=255)                                                                                                                  tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                     >  1:969274/1‑41 (MQ=255)                                                                                                                  tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                     <  2:969274/41‑1 (MQ=255)                                                                                                                   cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCggg                    <  2:440392/121‑1 (MQ=255)                                                                                                                   cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCggg                    >  1:440392/1‑121 (MQ=255)                                                                                                                   cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGccg  <  2:839622/139‑2 (MQ=255)                                                                                                                   cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa  <  1:1055079/139‑1 (MQ=255)                                                                                                                   cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa  <  2:7872/139‑1 (MQ=255)                                                                                                                     ggTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa  >  2:830535/1‑137 (MQ=255)                                                                                                                     ggTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa  <  1:830535/137‑1 (MQ=255)                                                                                                                            ggggAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGa                                                                                              >  2:857683/1‑38 (MQ=255)                                                                                                                               |                                                                                                                               CGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCA  >  NC_000913/4619787‑4620039

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