Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,112,456	Δ1 bp	coding (606/1086 nt)	rfbB ←	dTDP‑glucose 4,6‑dehydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,112,456	0	A	.	100.0%	87.1 / NA	19	coding (606/1086 nt)	rfbB	dTDP‑glucose 4,6‑dehydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (9/10);  total (9/10)

GAGCACTACATCGATGTTTTTCTTTTCGTTGTGCCCACCAATGTTATAAGTTTCACCCGCTTTACCTTCGGTTACGACGGTATATAACGCACGCGCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAAATGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTGTCGTCTCAGTAAATAAGGGTAATTCTTCTGTATTATTTACTTCATCTGGATGAGGCAAATCACCATAGACTTCGTCAGTAGAAATATGATGAAAACGGAAGCTATTTTTCTTGTCGCTA  >  NC_000913/2112249‑2112705                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          gAGCACTACATCGATGTTTTTCTTTTCGTTGTGCCCACCAATGTTATAAGTTTCACCCGCTTTACCTTCGGTTACGACGGTATATAACGCACGCGCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATATTTGTTCGAGCAAGTAGTCACAATTGTCgg                                                                                                                                                                                                               <  2:126752/251‑1 (MQ=255)   gCACTACATCGATGTTTTTCTTTTCGTTGTGCCCACCAATGTTATAAGTTTCACCCGCTTTACCTTCGGTTACGACGGTATATAACGCACGCGCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGt                                                                                                                                                                                                              >  1:646/1‑250 (MQ=255)             gATGTTTTTCTTTTCGTTGTGCCCACCAATGTTATAAGTTTCACCCGCTTTACCTTCGGTTACGACGGTATATAACGCACGCGCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGt                                                                                                                                                                                                   <  2:246734/251‑1 (MQ=255)                                                   tttCACCCGCTTTACCTTCGGTTACGACGGTATATAACGCACGCGCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCtt                                                                                                                                                             >  1:15027/1‑251 (MQ=255)                                                         ccGCTTTACCTTCGGTTACGACGGTATATAACGCACGCGCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCAAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATg                                                                                                                                                       >  1:325200/1‑251 (MQ=255)                                                                        ttACGACGGTATATAACGCACGCGCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGgctgc                                                                                                                                        >  1:370732/1‑251 (MQ=255)                                                                                  atataACGCACGCGCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTa                                                                                                                              <  1:300772/251‑1 (MQ=255)                                                                                         gcacgcGCATGATCTTCAACATACCACCAGTCGCGGCTCTGATCTCCTTTGCAATAAATAGGTCATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGAATAAATGATCGCTGGATGCTTTGGATGCGGCACAAGGGCTGCTTGGCGCGTCAGCTgtc                                                                                                                       >  2:367202/1‑251 (MQ=255)                                                                                         gcacgcGCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTgtc                                                                                                                       <  1:106757/251‑1 (MQ=255)                                                                                         gcacgcGAATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGCTGCGGAATAAGGGCTGCTTGGCGCGTAAGCTgtc                                                                                                                       >  1:278030/1‑251 (MQ=255)                                                                                           acgcTCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCAAAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTGTCCGTGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAATCTctcgt                                                                                                                     <  2:136837/251‑1 (MQ=255)                                                                                           acgcGCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTgtcgt                                                                                                                     <  2:136840/251‑1 (MQ=255)                                                                                                cATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTGTCGtctc                                                                                                                  >  1:27815/1‑249 (MQ=255)                                                                                                 aTGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACt                                                                                                                                                                                <  2:229486/186‑1 (MQ=255)                                                                                                 aTGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACt                                                                                                                                                                                >  1:229486/1‑186 (MQ=255)                                                                                                                                      ccTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTGTCGTCTCAGTAAATAAGGGTAATTCTTCTGTATTATTTACTTCATc                                                                           <  1:28270/250‑1 (MQ=255)                                                                                                                                                     ggTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTGTCGTCTCAGTAAATAAGGGTAATTCTTCTGTATTATTTACTTCATCTGGATGAGGCAAATCa                                                           <  2:646/251‑1 (MQ=255)                                                                                                                                                                     ccAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTGTCGTCTCAGTAAATAAGGGTAATTCTTCTGTATTATTTACTTCATCTGGATGAGGCAAATCACCATAGACTTCGTCa                                            <  1:155292/250‑1 (MQ=255)                                                                                                                                                                                 gAATAACCAGTGGAATAAGCTTTTCCGGGAA‑TGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTGTCGTCTCAGTAAATAAGGGTAATTCTTCTGTATTATTTACTTCATCTGGATGAGGCAAATCACCATAGACTTCGTCAGTAGAAATatgat                               >  2:54008/1‑251 (MQ=255)                                                                                                                                                                                                            gggaaTGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTGTCGTCTCAGTAAATAAGGGTAATTCTTCTGTATTATTTACTTCATCTGGATGAGGCAAATCACCATAGACTTCGTCAGTAGAAATATGATGAAAACGGAAGCTATTTTTCTTGTCg     <  1:89278/248‑1 (MQ=255)                                                                                                                                                                                                             ggaaTGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTGTCGTCTCAGTAAATAAGGGTAATTCTTCTGTATTATTTACTTCATCTGGATGAGGCAAATCACCATAGACTTCGTCAGTAGAAATATGATGAAAACGGAAGCTATTTTTCTTGTCGc    <  1:111003/249‑1 (MQ=255)                                                                                                                                                                                                               aaTGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTGTCGTCTCAGTAAATAAGGGTAATTCTTCTGTATTATTTACTTCATCTGGATGAGGCAAATCACCATAGACTTCGTCAGTAGAAATATGATGAAAACGGAAGCTATTTTTCTTGTCGCTa  <  1:238799/251‑1 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          GAGCACTACATCGATGTTTTTCTTTTCGTTGTGCCCACCAATGTTATAAGTTTCACCCGCTTTACCTTCGGTTACGACGGTATATAACGCACGCGCATGATCTTCAACATACAACCAGTCGCGGATCTGATCTCCTTTGCCATAAATAGGTAATGCCTTACCTTCCAGTGCATTAAGAATAACCAGTGGAATAAGCTTTTCCGGGAAATGATAAGGACCATAGTTGTTCGAGCAATTAGTCACAATTGTCGGTAAACCATATGTACGTTTCCACGCGCGGACTAAATGATCGCTGGATGCTTTGGATGCGGAATAAGGGCTGCTTGGCGCGTAAGCTGTCGTCTCAGTAAATAAGGGTAATTCTTCTGTATTATTTACTTCATCTGGATGAGGCAAATCACCATAGACTTCGTCAGTAGAAATATGATGAAAACGGAAGCTATTTTTCTTGTCGCTA  >  NC_000913/2112249‑2112705

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