Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,046,051	A→G	V632A (GTA→GCA)	gfcD ←	putative O‑antigen capsule production periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,046,051	0	A	G	100.0%	187.2 / NA	56	V632A (GTA→GCA)	gfcD	putative O‑antigen capsule production periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/29);  total (27/29)

GTGACAGCATTTGTCCACCATCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATACGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAG  >  NC_000913/1045905‑1046198                                                                                                                                                   |                                                                                                                                                    gTGACAGCATTTGTCCACCATCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGtt                                                                                                                                                 >  1:505975/1‑151 (MQ=255)    aCAGCATTTGTCCACCATCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACt                                                                                                                                              <  2:779641/151‑1 (MQ=255)    aCAGCATTTGTCCACCATCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACt                                                                                                                                              <  1:252135/151‑1 (MQ=255)        cATTTGTCCACCATCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGAtt                                                                                                                                          >  1:160441/1‑151 (MQ=255)        cATTTGTCCACCATCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGAtt                                                                                                                                          <  1:691180/151‑1 (MQ=255)         aTTTGTCCACCATCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTg                                                                                                                                         >  2:1625/1‑151 (MQ=255)         aTTTGTCCACCATCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTg                                                                                                                                         <  2:353988/151‑1 (MQ=255)           ttGTCCACCATCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTgc                                                                                                                                        <  2:690334/150‑1 (MQ=255)                 accaTCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCAt                                                                                                                                 >  2:86692/1‑151 (MQ=255)                   caTCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATAc                                                                                                                               <  1:394869/151‑1 (MQ=255)                       aCGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCt                                                                                                                           <  2:549324/151‑1 (MQ=255)                        cgcgTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCt                                                                                                                           >  2:412140/1‑150 (MQ=255)                           gTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCg                                                                                                                       >  2:436916/1‑151 (MQ=255)                                 ggTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCAcgc                                                                                                                 >  2:41288/1‑151 (MQ=255)                                      ccACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATCCCCCTACCGCCACGCCGCTg                                                                                                            <  2:364427/151‑1 (MQ=255)                                       cACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGt                                                                                                           >  1:728694/1‑151 (MQ=255)                                                 ccGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGtt                                                                                                  <  1:304335/150‑1 (MQ=255)                                                       ggttgcggttAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATCCCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGAc                                                                                           <  1:794908/151‑1 (MQ=255)                                                               ttAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATc                                                                                    <  1:599389/150‑1 (MQ=255)                                                                     ccAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTgcgc                                                                             >  2:589186/1‑151 (MQ=255)                                                                       aaTGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCt                                                                           >  2:897447/1‑151 (MQ=255)                                                                           gTCATAAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCAGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTAt                                                                       >  2:514405/1‑151 (MQ=255)                                                                                 aaGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTCCTGATTGCCGCCCATCCCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTg                                                                 <  1:784768/151‑1 (MQ=255)                                                                                 aaGTCAAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATCCCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCtt                                                                    <  1:518392/148‑1 (MQ=255)                                                                                     cGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGcc                                                               >  1:87737/1‑149 (MQ=255)                                                                                      gAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGcc                                                               >  1:262969/1‑148 (MQ=255)                                                                                      gAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCa                                                              <  2:657476/149‑1 (MQ=255)                                                                                       aaTGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGat                                                           <  2:160418/151‑1 (MQ=255)                                                                                           gAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTg                                                       >  2:553970/1‑151 (MQ=255)                                                                                                  gaTATAAAACCCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAAc                                                 <  1:370016/150‑1 (MQ=255)                                                                                                     ataAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCt                                              >  2:533540/1‑150 (MQ=255)                                                                                                     ataAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCt                                              >  2:645138/1‑150 (MQ=255)                                                                                                       aaaaaCCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTaa                                            <  1:273188/150‑1 (MQ=255)                                                                                                        aaaaCCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGt                                          >  1:468121/1‑151 (MQ=255)                                                                                                         aaaCCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATCCCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGtt                                         <  2:763005/151‑1 (MQ=255)                                                                                                         aaaCCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGtt                                         >  2:186947/1‑151 (MQ=255)                                                                                                           aCCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTtc                                       >  1:637297/1‑151 (MQ=255)                                                                                                            ccTTTGCTAAAGCCGCCTTCGCCGTAGTCTTCTTTAGATGCGTTTCTGTTTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTtca                                      <  2:317656/151‑1 (MQ=255)                                                                                                               ttACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTtcatc                                    <  1:792521/150‑1 (MQ=255)                                                                                                               ttACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTtcat                                     <  1:415641/149‑1 (MQ=255)                                                                                                                tACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTtcatc                                    >  1:272707/1‑149 (MQ=255)                                                                                                                  cTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTtcatca                                   <  2:124014/148‑1 (MQ=255)                                                                                                                  cTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGt                                 >  2:32095/1‑150 (MQ=255)                                                                                                                   tAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGt                                 <  2:628087/149‑1 (MQ=255)                                                                                                                   tAAAGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTAc                               >  1:613310/1‑151 (MQ=255)                                                                                                                    aaaGCCGCCTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACg                              >  2:684478/1‑151 (MQ=255)                                                                                                                       gccgccTTCGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGcc                            <  1:1625/150‑1 (MQ=255)                                                                                                                               cGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAg                      <  2:211937/148‑1 (MQ=255)                                                                                                                               cGCCGTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGt                   <  1:569405/151‑1 (MQ=255)                                                                                                                                   gTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTcggc               <  2:229730/151‑1 (MQ=255)                                                                                                                                   gTAGTCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTcggc               >  2:627730/1‑151 (MQ=255)                                                                                                                                       tCATCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGt           >  2:840224/1‑151 (MQ=255)                                                                                                                                          tCTTTAGATGCGTTACTGATTGCCGCCCATCCCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGtt          <  1:154171/149‑1 (MQ=255)                                                                                                                                          tCTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTcc        >  2:584558/1‑151 (MQ=255)                                                                                                                                           cTTTAGATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTcc        >  1:481106/1‑150 (MQ=255)                                                                                                                                                gATGCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAg  <  1:654980/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GTGACAGCATTTGTCCACCATCACGCGTCAATGGTGTCCACGAAACCACCGCGCGGTTGCGGTTAGGTCCAATGGTCATCAAGTCGAATGGAATCGAGATATAAAAACCTTTACTAAAGCCGCCTTCGCCGTAGTCATCTTTAGATACGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAG  >  NC_000913/1045905‑1046198

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