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%	146.1 / NA	46	A57T (GCA→ACA)	waaH	LPS(HepIII)‑glucuronic acid glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/26);  total (20/26)

TGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGCGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATTATATAAC  >  NC_000913/3789781‑3790043                                                                                                                                     |                                                                                                                                   tggGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATa                                                                                                                              >  2:55706/1‑139 (MQ=255)      tAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTg                                                                                                                         >  2:104312/1‑139 (MQ=255)           ttCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTAtt                                                                                                                    >  2:1081732/1‑139 (MQ=255)           ttCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTAtt                                                                                                                    >  2:37974/1‑139 (MQ=255)             cATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTc                                                                                                                  >  1:235804/1‑139 (MQ=255)             cATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTc                                                                                                                  >  1:598841/1‑139 (MQ=255)              aTCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCa                                                                                                                 >  1:867402/1‑139 (MQ=255)               tCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCaa                                                                                                                <  2:158019/139‑1 (MQ=255)                   cAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGa                                                                                                            >  1:1394/1‑139 (MQ=255)                    aGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGaa                                                                                                           <  1:353783/139‑1 (MQ=255)                      cATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAAtt                                                                                                         <  2:302621/139‑1 (MQ=255)                      cATCGACAAAAGCGACATATTTGCCCATTGCCACTTCAACCCCACGATTACGCGACACCGATGCACCAGCATTCGCCTGATGCAAAACACGAACGCGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAAtt                                                                                                         <  2:759568/139‑1 (MQ=255)                                aGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGaa                                                                                               >  1:172248/1‑139 (MQ=255)                                 gCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAAc                                                                                              <  1:55706/139‑1 (MQ=255)                                   gACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCa                                                                                            <  2:235804/139‑1 (MQ=255)                                   gACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCa                                                                                            <  1:812668/139‑1 (MQ=255)                                   gACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCa                                                                                            <  2:121386/139‑1 (MQ=255)                                                gTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAAtt                                                                                                         >  1:962971/1‑113 (MQ=255)                                                gTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAAtt                                                                                                         <  2:962971/113‑1 (MQ=255)                                                 ttGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCaa                                                                                                                <  1:391292/105‑1 (MQ=255)                                                 ttGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCaa                                                                                                                >  2:391292/1‑105 (MQ=255)                                                 ttGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTaataataat                                                                              >  2:21821/1‑139 (MQ=255)                                                  tGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATg                                                                             <  2:102623/139‑1 (MQ=255)                                                  tGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATg                                                                             >  1:102623/1‑139 (MQ=255)                                                            aTCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                   <  1:781403/139‑1 (MQ=255)                                                             tCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGc                                                                  >  1:438967/1‑139 (MQ=255)                                                                  cGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACa                                                                                                              <  1:620687/90‑1 (MQ=255)                                                                  cGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACa                                                                                                              >  2:620687/1‑90 (MQ=255)                                                                     ttACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                   <  1:915466/130‑1 (MQ=255)                                                                     ttACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                   >  2:915466/1‑130 (MQ=255)                                                                        cgcgCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                   >  1:724139/1‑127 (MQ=255)                                                                        cgcgCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                   <  1:714330/127‑1 (MQ=255)                                                                        cgcgCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCATTGAACCATCGTTAATAATAATGATTTCCagag                                                                   >  2:714330/1‑127 (MQ=255)                                                                        cgcgCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGCAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                   <  2:724139/127‑1 (MQ=255)                                                                           gcCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCg                                                    <  1:37974/139‑1 (MQ=255)                                                                                cGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAAtt                                               <  1:206423/139‑1 (MQ=255)                                                                                       ccAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTTTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGtt                                                                                       >  2:640972/1‑92 (MQ=255)                                                                                       ccAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGtt                                                                                       <  1:640972/92‑1 (MQ=255)                                                                                                               cgaacgTGCGGATAGTTTTCTGTGTAATACTTTGGTATTTCAACAGAATTATCCGCTGAACCATTGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGc                <  2:861784/139‑1 (MQ=255)                                                                                                               cgaacgTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGaa                                                                                               >  1:82701/1‑60 (MQ=255)                                                                                                               cgaacgTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGaa                                                                                               <  2:82701/60‑1 (MQ=255)                                                                                                               cgaacgTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGc                <  2:867402/139‑1 (MQ=255)                                                                                                                gaacgTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGcc               <  1:342824/139‑1 (MQ=255)                                                                                                                  acgTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCg             <  1:104312/139‑1 (MQ=255)                                                                                                                           tAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATtatata    <  2:34175/139‑1 (MQ=255)                                                                                                                             gTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATTATATAAc  <  2:598841/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                   TGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGCGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATTATATAAC  >  NC_000913/3789781‑3790043

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