Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,670	T→C	G356G (GGT→GGC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,670	0	T	C	100.0%	166.7 / NA	45	G356G (GGT→GGC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/16);  total (29/16)

CTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCC  >  NC_000913/4620522‑4620817                                                                                                                                                     |                                                                                                                                                    ctTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCAc                                                                                                                                                   >  1:21293/1‑151 (MQ=255) ctTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCAc                                                                                                                                                   >  2:475359/1‑151 (MQ=255)   tCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACt                                                                                                                                                  <  1:238303/150‑1 (MQ=255)    ccTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGa                                                                                                                                                <  1:216487/151‑1 (MQ=255)        gggCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACcac                                                                                                                                            >  2:502587/1‑151 (MQ=255)                          cGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCg                                                                                                                          >  1:99248/1‑151 (MQ=255)                              ggTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTAc                                                                                                                      >  1:638184/1‑151 (MQ=255)                                      cgcgGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTatat                                                                                                              <  1:513182/151‑1 (MQ=255)                                              tGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGa                                                                                                       >  2:291269/1‑150 (MQ=255)                                               ggAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGaaa                                                                                                     <  1:109593/151‑1 (MQ=255)                                                       tCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAAc                                                                                              <  1:427928/150‑1 (MQ=255)                                                       tCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAAc                                                                                              <  1:18507/150‑1 (MQ=255)                                                             gccgTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTc                                                                                       >  1:203850/1‑151 (MQ=255)                                                             gccgTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTc                                                                                       >  1:644609/1‑151 (MQ=255)                                                               cgTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCAc                                                                                     >  1:426834/1‑151 (MQ=255)                                                                   tGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTggg                                                                                 >  1:320069/1‑151 (MQ=255)                                                                   tGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTggg                                                                                 >  2:279819/1‑151 (MQ=255)                                                                   tGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTggg                                                                                 >  2:143880/1‑151 (MQ=255)                                                                    gCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTg                                                                                   <  2:667963/148‑1 (MQ=255)                                                                     ccGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTc                                                                               >  2:603607/1‑151 (MQ=255)                                                                              aTGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGaaa                                                                      <  2:58187/151‑1 (MQ=255)                                                                                 tctcTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCt                                                                   >  2:479822/1‑151 (MQ=255)                                                                                    ctgctgCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTcg                                                                >  2:673323/1‑151 (MQ=255)                                                                                     tgctgcGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTcgc                                                               >  2:322364/1‑151 (MQ=255)                                                                                      gctgcGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTcg                                                                <  1:448146/149‑1 (MQ=255)                                                                                          cgcgATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGata                                                          <  1:508953/151‑1 (MQ=255)                                                                                          cgcgATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGat                                                           >  1:278589/1‑150 (MQ=255)                                                                                           gcgATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGatat                                                         >  1:139442/1‑151 (MQ=255)                                                                                           gcgATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGatat                                                         >  2:480167/1‑151 (MQ=255)                                                                                                gacgacATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGc                                                     <  2:278584/150‑1 (MQ=255)                                                                                                 acgacATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGcc                                                    >  2:274334/1‑150 (MQ=255)                                                                                                      aTCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGActc                                              <  2:202626/151‑1 (MQ=255)                                                                                                           gATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGc                                          >  2:484755/1‑150 (MQ=255)                                                                                                                 cACCGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAatat                                   >  1:594983/1‑151 (MQ=255)                                                                                                                   ccGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATAttt                                 >  1:662589/1‑151 (MQ=255)                                                                                                                    cGCTGACCACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATAtttt                                >  1:301635/1‑151 (MQ=255)                                                                                                                          ccACGGTTGCGATCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACt                          >  1:200440/1‑151 (MQ=255)                                                                                                                                      tCCGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGa               >  1:409831/1‑150 (MQ=255)                                                                                                                                       ccGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAat             >  1:342191/1‑151 (MQ=255)                                                                                                                                       ccGACCTGGACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAat             >  2:339635/1‑151 (MQ=255)                                                                                                                                              ggACCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCaa      <  2:626202/151‑1 (MQ=255)                                                                                                                                                aCCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCaa      >  1:239069/1‑149 (MQ=255)                                                                                                                                                aCCGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAg    <  2:178862/151‑1 (MQ=255)                                                                                                                                                  cGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGcc  <  1:38159/151‑1 (MQ=255)                                                                                                                                                  cGGCACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGcc  <  1:617633/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCC  >  NC_000913/4620522‑4620817

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