Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,663	G→A	W354* (TGG→TAG)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,663	0	G	A	100.0%	110.3 / NA	32	W354* (TGG→TAG)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/17);  total (15/17)

TTCCTGGGGCCACCGTCGCGACGTCGC‑‑CGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATG  >  NC_000913/4620523‑4620802                                                                                                                                               |                                                                                                                                            ttCCTGGGGCCACCGTCGCGACGTCGC‑‑CGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCgg                                                                                                                                       <  5:18589/147‑1 (MQ=255)           cACCGTCGCGACGTCGC‑‑CGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCacac                                                                                                                           <  8:60072/149‑1 (MQ=255)            taaccggcgACGTCGCGACGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGAc                                                                                                                                >  7:83298/7‑145 (MQ=255)            taaccggcgACGTCGCGACGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGAGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGAc                                                                                                                                <  8:83298/139‑1 (MQ=255)                               ggTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGt                                                                                                         <  3:60739/149‑1 (MQ=255)                                             gcTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGcc                                                                                           <  6:200427/148‑1 (MQ=255)                                              cTGGAACTGTTCGACCGCCGTCTGCGGGAGCTGATGTCTCTGCTGCGCGTTGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTAc                                                                                                       <  8:142785/136‑1 (MQ=255)                                              cTGGAACTGTTCGACCGCCGTATGCCGGAGCTGATGTCTCTGCTGCGAGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACCCGCGTGAACACATTCCGGTAc                                                                                                       >  7:142785/1‑136 (MQ=255)                                                                      ccGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTgg                                                                   <  5:204820/148‑1 (MQ=255)                                                                      ccGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTg                                                                    <  1:225961/147‑1 (MQ=255)                                                                           gCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATc                                                             >  3:134914/1‑149 (MQ=255)                                                                            cTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGt                                                                                                   <  5:132793/110‑1 (MQ=255)                                                                            cTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGt                                                                                                   >  6:132793/1‑110 (MQ=255)                                                                              gATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGt                                                           >  5:68757/1‑148 (MQ=255)                                                                               aTGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTg                                                          >  6:167715/1‑148 (MQ=255)                                                                                          gcgcGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTcgcg                                               >  2:229748/1‑148 (MQ=255)                                                                                                tGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCg                                        <  4:134914/149‑1 (MQ=255)                                                                                                    gacATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGc                                      >  8:23594/1‑147 (MQ=255)                                                                                                           tGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGc                                                                                            <  4:52271/86‑1 (MQ=255)                                                                                                           tGATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGc                                                                                            >  3:52271/1‑86 (MQ=255)                                                                                                            gATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGc                                      >  8:80251/1‑139 (MQ=255)                                                                                                            gATCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGc                                      <  7:80251/139‑1 (MQ=255)                                                                                                             aTCCTCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTgg                            <  6:68757/148‑1 (MQ=255)                                                                                                                 tCACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCaaaa                       >  8:21072/1‑149 (MQ=255)                                                                                                                  cACCGCTGACCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCaaaa                       >  6:27809/1‑148 (MQ=255)                                                                                                                          aCCACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTAGAGGGCCCGAAAGTAAAACCGGGCTCACTGGGTCAGCGTGAAACCTTCGCGGATATCGGCCTGACTCTGGCAAAATATTTTg                <  8:157611/147‑1 (MQ=255)                                                                                                                            cACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGt                                                                                    <  5:207649/77‑1 (MQ=255)                                                                                                                            cACGGTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGt                                                                                    >  6:207649/1‑77 (MQ=255)                                                                                                                                gTTGCGATCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCt        <  1:229748/149‑1 (MQ=255)                                                                                                                                      aTCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATg  >  5:28407/1‑149 (MQ=255)                                                                                                                                       tCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCgg                                                                                                          <  6:213995/44‑1 (MQ=255)                                                                                                                                       tCCGACCTAGACCGGTACTGACCACACGCGTGAACACATTCCgg                                                                                                          >  5:213995/1‑44 (MQ=255)                                                                                                                                               |                                                                                                                                            TTCCTGGGGCCACCGTCGCGACGTCGC‑‑CGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATG  >  NC_000913/4620523‑4620802

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