Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,796,782	G→A	T268T (ACG→ACA)	waaC →	ADP‑heptose:LPS heptosyltransferase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,796,782	0	G	A	100.0%	93.6 / NA	27	T268T (ACG→ACA)	waaC	ADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/14);  total (13/14)

TGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTT  >  NC_000913/3796496‑3796889                                                                                                                                                                                                                                                                                               |                                                                                                            tGCCTACAGATGCTGTCGAATATGCCGTATTTCTTCATGCGACGACCCTTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGAGTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGAt                                                                                                 <  2:306323/299‑1 (MQ=255) tGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTACTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGAt                                                                                                 <  2:50924/299‑1 (MQ=255)  gcccaCAGATGCTGGNCAATATGCCGTATTTCTTNATGCGACCACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTNAACAGCGGCACTGGAt                                                                                                 <  1:585857/294‑1 (MQ=255)  gCCTACAGATGCTGGCNAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGNTACGGGGTTAAGCCATTTNACAGCGGCACTGGATa                                                                                                >  1:582993/1‑299 (MQ=255)  gCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTC‑‑GAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGAc                                                                                             >  1:42519/1‑300 (MQ=255)    cTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGAc                                                                                             >  1:173178/1‑300 (MQ=255)     tACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGGGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGAcc                                                                                            <  1:348600/300‑1 (MQ=255)     tACAGATGCTGGCGAATATGCCGTATGTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGAcc                                                                                            <  1:220446/300‑1 (MQ=255)          aTGCTGGCGAATATGCCGTAGTTCTTCCTGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGTCGAGAATTGATTGTTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGa                                                                                              <  2:40385/293‑1 (MQ=255)                        gCCGTATTTCTTCATGCGACGACCCGGGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGTGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGa                                                                         <  2:290429/300‑1 (MQ=255)                        gCCGTAGTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGTAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATCCATGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGTCTAAACGACTGGCGTAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGTGGCTAAATTTTTATTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGa                                                                         <  2:290267/300‑1 (MQ=255)                         ccGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGAc                                                                        >  1:415561/1‑300 (MQ=255)                           gTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGAcct                                                                      >  1:686458/1‑299 (MQ=255)                           gTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCAAGGTTTATGGAcct                                                                      >  1:686470/1‑299 (MQ=255)                                         gacgacCCTTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGAt                                                         <  1:162775/299‑1 (MQ=255)                                         gacgacCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACNGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGAt                                                         <  1:135315/299‑1 (MQ=255)                                             acCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCTGTGGGGCGCGCCGCATGAGGGAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTACTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGACACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAAtt                                                    <  2:51849/300‑1 (MQ=255)                                                      gatAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTAt                                           >  1:330467/1‑300 (MQ=255)                                                      gatAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTAt                                           >  2:97857/1‑300 (MQ=255)                                                      gatAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCAGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCAGGGATTAGTTGGTGGGTAt                                           >  2:65401/1‑300 (MQ=255)                                                           acacTGGCCGTAAGAACACTGGCGAGAATTGATTGTTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGTCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGgaa                                      <  2:586002/300‑1 (MQ=255)                                                           acacTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCAGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCAAGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGgaa                                      >  2:444380/1‑300 (MQ=255)                                                                    ggAAGAACACTGTCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGAGTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGAt                             <  2:5656/300‑1 (MQ=255)                                                                     gaagaaCACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGAt                             >  2:694002/1‑299 (MQ=255)                                                                               ggCGAGGATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCAATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTTTGTa                     >  2:358408/1‑297 (MQ=255)                                                                                       ttgattgTTTTACTGGCTGATTCATGAAGACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGTGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTTTGTAGGGCTccgggg          <  2:201325/300‑5 (MQ=255)                                                                                                ttACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTANTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCNTGCTGGGCGGGGCTAAATTTGTAGTGTCAGTGGATACGGGGTTAAGCCATTTAACAGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCTACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTTTGTAGGGCTccggggaatgagtt  >  2:537383/1‑287 (MQ=255)                                                                                                                                                                                                                                                                                               |                                                                                                            TGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTT  >  NC_000913/3796496‑3796889

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