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

ATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGCGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATTATATAACGGAA  >  NC_000913/3789777‑3790047                                                                                                                                         |                                                                                                                                       aTGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTgtgt                                                                                                                                      >  1:21844/1‑139 (MQ=255)    gtggGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAAt                                                                                                                                   >  1:435618/1‑139 (MQ=255)     tggGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATa                                                                                                                                  >  1:147254/1‑139 (MQ=255)                  aTCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCa                                                                                                                     >  1:357044/1‑139 (MQ=255)                   tCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACGAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCaa                                                                                                                    <  1:264373/139‑1 (MQ=255)                   tCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCaa                                                                                                                    <  1:1021930/139‑1 (MQ=255)                     gTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACa                                                                                                                  <  1:848099/139‑1 (MQ=255)                      tCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAg                                                                                                                 >  1:990870/1‑139 (MQ=255)                                    aGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGaa                                                                                                   <  2:163619/139‑1 (MQ=255)                                    aGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGaa                                                                                                   <  2:147254/139‑1 (MQ=255)                                    aGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGaa                                                                                                   <  1:996173/139‑1 (MQ=255)                                      cGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAAcc                                                                                                 <  1:341249/139‑1 (MQ=255)                                       gACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCa                                                                                                <  1:150919/139‑1 (MQ=255)                                        aCATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCAt                                                                                               >  2:720642/1‑139 (MQ=255)                                           tatTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGt                                                                                            <  1:887145/139‑1 (MQ=255)                                                     ttGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTaataataat                                                                                  >  2:356897/1‑139 (MQ=255)                                                      tGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATg                                                                                 >  1:389332/1‑139 (MQ=255)                                                      tGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATg                                                                                 >  2:148776/1‑139 (MQ=255)                                                      tGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGCTGAACCATCGTTAATAATAATg                                                                                 <  2:1025879/139‑1 (MQ=255)                                                       gCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCaa                                                                                                                    >  2:1034984/1‑103 (MQ=255)                                                       gCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCaa                                                                                                                    <  1:1034984/103‑1 (MQ=255)                                                                aTCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCa                                                                                                                     >  2:236418/1‑93 (MQ=255)                                                                aTCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCa                                                                                                                     <  1:236418/93‑1 (MQ=255)                                                                aTCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                       <  1:356897/139‑1 (MQ=255)                                                                aTCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                       <  2:357044/139‑1 (MQ=255)                                                                aTCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                       <  1:117631/139‑1 (MQ=255)                                                                aTCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                       <  1:258430/139‑1 (MQ=255)                                                                      cGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGtt                                                                                           <  1:869561/113‑1 (MQ=255)                                                                      cGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGtt                                                                                           >  2:869561/1‑113 (MQ=255)                                                                      cGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                       >  1:915729/1‑133 (MQ=255)                                                                      cGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                       <  2:915729/133‑1 (MQ=255)                                                                      cGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTcc                                                                 >  2:1015967/1‑139 (MQ=255)                                                                      cGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTcc                                                                 >  1:723950/1‑139 (MQ=255)                                                                              cgcCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGc                                                         <  1:219509/139‑1 (MQ=255)                                                                                  aCCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAAcc                                                                                                 <  1:1007435/95‑1 (MQ=255)                                                                                  aCCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAAcc                                                                                                 >  2:1007435/1‑95 (MQ=255)                                                                                   ccGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAAt                                                    >  1:571401/1‑139 (MQ=255)                                                                                      aTGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTaa                                                 <  2:21844/139‑1 (MQ=255)                                                                                           ccAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGAtt                                            >  2:803808/1‑139 (MQ=255)                                                                                           ccAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGAtt                                            >  2:10866/1‑139 (MQ=255)                                                                                           ccAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGAtt                                            >  1:402690/1‑139 (MQ=255)                                                                                              gCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCa                                         <  2:855324/139‑1 (MQ=255)                                                                                                 ttCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATAc                                      <  2:571401/139‑1 (MQ=255)                                                                                                 ttCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATAc                                      <  2:990870/139‑1 (MQ=255)                                                                                                         aTGCAACAAACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCgg                              >  2:598520/1‑139 (MQ=255)                                                                                                              acaaACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                       <  1:88139/93‑1 (MQ=255)                                                                                                              acaaACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCagag                                                                       >  2:88139/1‑93 (MQ=255)                                                                                                              acaaACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAA‑CATCGTTAATAATAATGATTTCCagag                                                                       <  2:293573/92‑1 (MQ=255)                                                                                                              acaaACGAACGTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAA‑CATCGTTAATAATAATGATTTCCagag                                                                       >  1:293573/1‑92 (MQ=255)                                                                                                                   cgaacgTGCTGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGc                    <  1:720642/139‑1 (MQ=255)                                                                                                                   cgaacgTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGc                    <  2:435618/139‑1 (MQ=255)                                                                                                                   cgaacgTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGc                    >  1:323595/1‑139 (MQ=255)                                                                                                                     aacgTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGccc                  >  2:336839/1‑139 (MQ=255)                                                                                                                      acgTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCg                 <  2:323595/139‑1 (MQ=255)                                                                                                                      acgTGCGGATAGTTTTCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCg                 <  1:10866/139‑1 (MQ=255)                                                                                                                                    ttCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATTATATAACGGa   >  1:592568/1‑139 (MQ=255)                                                                                                                                     tCTGTGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATTATATAACGGaa  <  1:803808/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       ATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGCGTAATACTTTGCTATTTCAACAGAATTATCCGTTGAACCATCGTTAATAATAATGATTTCCAGAGCAGTCCAGGTTTGCGTAATTAAAGATTCCATACAAGTGCGGAAATCATCGCCCGCATTATATAACGGAA  >  NC_000913/3789777‑3790047

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