Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,117,670	T→C	A245A (GCA→GCG)	wcaK ←	putative colanic acid biosynthesis pyruvyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,117,670	0	T	C	100.0%	206.0 / NA	58	A245A (GCA→GCG)	wcaK	putative colanic acid biosynthesis pyruvyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/23);  total (35/23)

TTGTTATAGCTGTCAATGCCCGTACAGGTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCAC  >  NC_000913/2117536‑2117804                                                                                                                                       |                                                                                                                                       ttGTTATAGCTGTCAATGCCCGTACAGGTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAg                                                                                                                                   <  1:393386/140‑1 (MQ=255)  tGTTATAGCTGTCAATGCCCGTACAGGTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGt                                                                                                                                  <  2:923956/140‑1 (MQ=255)  tGTTATAGCTGTCAATGCCCGTACAGGTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGt                                                                                                                                  >  1:633729/1‑140 (MQ=255)            gTCAATGCCCGTACAGGTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAggg                                                                                                                        >  1:802063/1‑140 (MQ=255)            gTCAATGCCCGTACAGGTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAggg                                                                                                                        >  1:51031/1‑140 (MQ=255)             tCAATGCCCGTACAGGTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGt                                                                                                                       >  1:547360/1‑140 (MQ=255)                    ccGTACAGGTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGcc                                                                                                                <  2:929055/140‑1 (MQ=255)                    ccGTACAGGTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGcc                                                                                                                <  2:486262/140‑1 (MQ=255)                          acaggtggagCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGtt                                                                                                          <  2:559161/133‑1 (MQ=255)                           ggTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGttt                                                                                                         >  2:694388/1‑140 (MQ=255)                                 gagCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCtgtt                                                                                                   >  2:863043/1‑140 (MQ=255)                                  agCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCtgttg                                                                                                  >  1:184863/1‑140 (MQ=255)                                    cgcAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTtgtg                                                                                                >  1:838966/1‑140 (MQ=255)                                      cAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGtt                                                                                                          >  2:695995/1‑128 (MQ=255)                                      cAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGtt                                                                                                          <  1:695995/128‑1 (MQ=255)                                         tCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCa                                                                                           >  1:946387/1‑140 (MQ=255)                                          cACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCaa                                                                                          <  2:51031/140‑1 (MQ=255)                                                 tACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTccag                                                                                   <  2:802063/140‑1 (MQ=255)                                                 tACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGACACCGTTTTCTGTTGTGCGGCAACGTccag                                                                                   >  1:243763/1‑140 (MQ=255)                                                     cTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTccagccag                                                                               <  2:633729/140‑1 (MQ=255)                                                       tCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAgtg                                                                             >  2:586391/1‑140 (MQ=255)                                                       tCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAgtg                                                                             >  1:442719/1‑140 (MQ=255)                                                       tCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAgtg                                                                             >  1:379077/1‑140 (MQ=255)                                                       tCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGGAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAgtg                                                                             >  1:208039/1‑140 (MQ=255)                                                        cATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAgtgt                                                                            >  1:144540/1‑140 (MQ=255)                                                           cGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGa                                                                         >  2:870622/1‑140 (MQ=255)                                                           cGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGa                                                                         <  2:144540/140‑1 (MQ=255)                                                            gagaATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGaa                                                                        <  2:184863/140‑1 (MQ=255)                                                             agaATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAAc                                                                       >  2:220276/1‑140 (MQ=255)                                                              gaATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACg                                                                      >  2:688678/1‑140 (MQ=255)                                                                aTGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGc                                                                    >  2:673858/1‑140 (MQ=255)                                                                          ccaccCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCg                                                          <  1:870622/140‑1 (MQ=255)                                                                                  gCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCg                                                          <  1:569546/132‑1 (MQ=255)                                                                                  gCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCg                                                          >  2:569546/1‑132 (MQ=255)                                                                                        gCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCtgt                                            >  2:647506/1‑140 (MQ=255)                                                                                        gCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCtgt                                            >  1:251584/1‑140 (MQ=255)                                                                                            tttCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTgtgg                                        >  2:903134/1‑140 (MQ=255)                                                                                                  aCGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCg                                  >  2:443968/1‑140 (MQ=255)                                                                                                  aCGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCg                                  >  2:203124/1‑140 (MQ=255)                                                                                                        gttgAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGc                            <  1:688678/140‑1 (MQ=255)                                                                                                        gttgAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGc                            >  2:256306/1‑140 (MQ=255)                                                                                                          tgAGTGGTGCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCa                          <  1:647506/140‑1 (MQ=255)                                                                                                             gtggtgCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCAcgc                       <  2:774782/140‑1 (MQ=255)                                                                                                              tggtgCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCa                          <  1:124639/136‑1 (MQ=255)                                                                                                              tggtgCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCa                          >  2:124639/1‑136 (MQ=255)                                                                                                                 tgCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTa                   <  1:278169/140‑1 (MQ=255)                                                                                                                  gCCGAGACGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTAt                  >  2:764798/1‑140 (MQ=255)                                                                                                                       gaCGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACg             >  2:809933/1‑140 (MQ=255)                                                                                                                        aCGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCAcgcg                      >  2:592313/1‑130 (MQ=255)                                                                                                                        aCGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCAcgcg                      <  1:592313/130‑1 (MQ=255)                                                                                                                        aCGTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGc            <  1:863043/140‑1 (MQ=255)                                                                                                                          gTTTGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCa          <  2:838966/140‑1 (MQ=255)                                                                                                                           tttGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCAt         >  1:399292/1‑140 (MQ=255)                                                                                                                           tttGTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCAt         >  1:886067/1‑140 (MQ=255)                                                                                                                              gTCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGtt      <  1:220276/140‑1 (MQ=255)                                                                                                                               tCAAACGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTc     <  2:886067/140‑1 (MQ=255)                                                                                                                                  aaCGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCAc  >  1:90960/1‑140 (MQ=255)                                                                                                                                  aaCGGCGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCAc  >  1:153774/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                       TTGTTATAGCTGTCAATGCCCGTACAGGTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCAC  >  NC_000913/2117536‑2117804

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