Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,789,913	C→T	A57T (GCA→ACA)	waaH ←	LPS(HepIII)‑glucuronic acid glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,789,913	0	C	T	100.0%	104.0 / NA	33	A57T (GCA→ACA)	waaH	LPS(HepIII)‑glucuronic acid glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/20);  total (13/20)

TGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGCGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATTATATAACGG  >  NC_000913/3789778‑3790045                                                                                                                                        |                                                                                                                                     tggtggGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTa                                                                                                                                   <  1:628063/139‑1 (MQ=255)        aTAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACttt                                                                                                                            >  2:414607/1‑139 (MQ=255)            aCTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTa                                                                                                                        >  1:764163/1‑139 (MQ=255)                 aTCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCa                                                                                                                   >  1:352233/1‑139 (MQ=255)                 aTCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCa                                                                                                                   <  1:552279/139‑1 (MQ=255)                                      gACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCa                                                                                              <  2:303353/139‑1 (MQ=255)                                      gACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCa                                                                                              <  1:414607/139‑1 (MQ=255)                                                  cGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTtaataata                                                                                  <  1:9713/139‑1 (MQ=255)                                                    ttGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTaataataat                                                                                >  2:690099/1‑139 (MQ=255)                                                      gCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATg                                                                               >  1:150456/1‑138 (MQ=255)                                                      gCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATg                                                                               <  2:150456/138‑1 (MQ=255)                                                               aTCCCATGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGTACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAatgatg                                                                               <  2:539443/129‑1 (MQ=255)                                                               aTCCCATGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATg                                                                               >  1:539443/1‑129 (MQ=255)                                                                tCCCACGATTACGCGCCACCGATGTGCCAGCATTCGCCTGATGCAACAAACCAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCa                                                                                              >  1:512817/1‑113 (MQ=255)                                                                tCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACCTTGCTATTTCAACAGAATTATCCGTTGAACCa                                                                                              <  2:512817/113‑1 (MQ=255)                                                                    aCGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACttt                                                                                                                            >  2:769055/1‑79 (MQ=255)                                                                    aCGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACttt                                                                                                                            <  1:769055/79‑1 (MQ=255)                                                                                cACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                     >  2:136112/1‑122 (MQ=255)                                                                                cACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                     <  1:136112/122‑1 (MQ=255)                                                                                     aTGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCa                                                                                              >  1:749045/1‑92 (MQ=255)                                                                                     aTGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCa                                                                                              <  2:749045/92‑1 (MQ=255)                                                                                      tGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTaaa                                              <  1:690099/139‑1 (MQ=255)                                                                                        cgcCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAAc                                                                                                                 <  1:434033/70‑1 (MQ=255)                                                                                        cgcCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAAc                                                                                                                 >  2:434033/1‑70 (MQ=255)                                                                                               aTTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATa                                     <  1:324874/139‑1 (MQ=255)                                                                                                        aTGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTaaa                                              <  1:722403/121‑1 (MQ=255)                                                                                                        aTGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTaaa                                              >  2:722403/1‑121 (MQ=255)                                                                                                             acaaACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAAcc                                                                                               <  2:195949/67‑1 (MQ=255)                                                                                                             acaaACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAAcc                                                                                               >  1:195949/1‑67 (MQ=255)                                                                                                                     acgTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTCCGGAAATCATCGCCCg               <  2:436007/139‑1 (MQ=255)                                                                                                                       gTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCa             <  2:147092/139‑1 (MQ=255)                                                                                                                                gTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATTATATAAc    <  1:178358/139‑1 (MQ=255)                                                                                                                                  tttCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATTATATAACgg  <  1:284534/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     TGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGCGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATTATATAACGG  >  NC_000913/3789778‑3790045

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