Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,789,563	A→G	R173R (CGT→CGC)	waaH ←	UDP‑glucuronate:LPS(HepIII) glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,789,563	0	A	G	100.0%	65.9 / NA	21	R173R (CGT→CGC)	waaH	UDP‑glucuronate:LPS(HepIII) glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/4);  total (17/4)

CAGGCGGGTAATCTTAATATAGTGACGTTGATAATTAAGGTTTTTATTCCCTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCG  >  NC_000913/3789347‑3789806                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    cagGCGGGTAATCTTAATATAGTGACGTTGATAATTAAGGTTTTTATTCCCTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATACACCCCCATCCAGACAACGTGAGTcc                                                                                                                                                                                                                   >  1:76236/1‑251 (MQ=255)                       tGACGTTGATAATTAAGGTTTTTATTCCCTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAACGCCCCAt                                                                                                                                                                                             >  1:157897/1‑251 (MQ=255)                                                tCCCTTGTCTATGCATCCGACTCACTGACGTATTATGCAGATAATCTTTATATCATGATTGCTCGGTGTATCGCGCACGCAGCGCATTAAACATGCATTCTGTTGTCCAGACAATATCCTGATGATGAAATCCGGCAATAACTTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGCAGAAAGCCCCATCCGCCGCCAGTCCGGGCCGGATAAt                                                                                                                                                                    >  2:394528/1‑251 (MQ=255)                                                  ccTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATAc                                                                                                                                                                  >  1:292960/1‑251 (MQ=255)                                                                     tCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACACCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCg                                                                                                                                                >  2:353798/1‑250 (MQ=255)                                                                         tGACGTATTATGCCGATAATATTTATATAATGATTGCTCGGTATATCGCGCCCGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGCTGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACCATACCATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGACAGCCCCATCCGCAGCCAGTCCGGGCCGGTTACTACACCGGTTGAGCGAAGGCGATCgc                                                                                                                                           >  1:214358/1‑250 (MQ=255)                                                                                 tatGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGggatg                                                                                                                                   >  1:129354/1‑251 (MQ=255)                                                                                    gCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGggatggat                                                                                                                                >  2:368637/1‑251 (MQ=255)                                                                                         taataTTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACAGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGcc                                                                                                                            >  1:314761/1‑250 (MQ=255)                                                                                                  tataATGATTGCTCGGTATATCGAGCAAGCAGCGCCTTAAACCTGAATTCTGTTCTCCAGACAATATCCTGCTGATGTAATACGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCAAACGGCGCGAAGAAAGCCCCCACCACAGCCAGTCCGGGCAGGTTAATACACCGGATGAGCGCAACCGCTCGCTGTGGATGGCTTGCCAGGTTtcccc                                                                                                                  >  2:382521/1‑248 (MQ=255)                                                                                                    taATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCg                                                                                                                >  1:270491/1‑251 (MQ=255)                                                                                                                                       tAAACATGAAGTCTGTTGTCCAGACAATAGCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAAGAACATCGCGGCGATAAACCCCCATCCAGCCAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACAGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTgcgc                                                                              <  1:403628/250‑1 (MQ=255)                                                                                                                                             tGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCg                                                                       >  1:118638/1‑251 (MQ=255)                                                                                                                                                                     ccTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGtc                                               <  1:290308/251‑1 (MQ=255)                                                                                                                                                                           gatgTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGc                                         >  2:312700/1‑251 (MQ=255)                                                                                                                                                                                           aaaTTTAATGTTATTTTTCACAATAACAGCGCGGCGATAAACCCCCAGCCAGACAACGTGATGCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGTATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCTCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTGTCGTACAtggtgt                         <  2:287523/251‑2 (MQ=255)                                                                                                                                                                                                 aaTGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGa                   <  1:353798/251‑1 (MQ=255)                                                                                                                                                                                                     ttATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTc               >  2:321699/1‑251 (MQ=255)                                                                                                                                                                                                               cAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGc      >  2:330027/1‑250 (MQ=255)                                                                                                                                                                                                               cAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGATAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGc      >  1:143355/1‑250 (MQ=255)                                                                                                                                                                                                                   aaCATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCg  >  1:373896/1‑250 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    CAGGCGGGTAATCTTAATATAGTGACGTTGATAATTAAGGTTTTTATTCCCTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCG  >  NC_000913/3789347‑3789806

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