Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,130,344	(A)5→4	coding (535/1218 nt)	wcaD ←	putative colanic acid polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,130,340	0	A	.	100.0%	245.4 / NA	57	coding (539/1218 nt)	wcaD	putative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (27/30);  total (27/30)

AACAGGATAAAGGTCATAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAAT  >  NC_000913/2130208‑2130473                                                                                                                                     |                                                                                                                                      aaCAGGATAAAGGTCATAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGaaaatgcc                                                                                                                                >  1:620026/1‑136 (MQ=255)   cAGGATAAAGGTCATAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCgg                                                                                                                             >  1:73116/1‑140 (MQ=255)   cAGGATAAAGGTCATAACCCCTGAAAACGATCCGGATAATATTAT‑CCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGt                                                                                                                            <  2:928955/140‑1 (MQ=255)         aaaGGTCATAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGaa                                                                                                                       >  2:32101/1‑140 (MQ=255)          aaGGTCATAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGaaa                                                                                                                      <  1:15354/140‑1 (MQ=255)           aGGTCATAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAAt                                                                                                                     >  2:295122/1‑140 (MQ=255)                aTAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGaaa                                                                                                                      <  2:385272/134‑1 (MQ=255)                aTAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGaaa                                                                                                                      >  1:385272/1‑134 (MQ=255)                aTAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATa                                                                                                                    <  1:111119/136‑1 (MQ=255)                aTAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATa                                                                                                                    >  2:111119/1‑136 (MQ=255)                aTAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGt                                                                                                                >  1:1031632/1‑140 (MQ=255)                aTAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGt                                                                                                                >  2:978883/1‑140 (MQ=255)                aTAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGt                                                                                                                >  1:358234/1‑140 (MQ=255)                 tAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTg                                                                                                               <  1:939529/140‑1 (MQ=255)                     cccTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCATAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGt                                                                                                           >  2:959144/1‑140 (MQ=255)                         gAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTg                                                                                                       >  2:291476/1‑140 (MQ=255)                               gATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTgcg                                                                                                     >  2:938130/1‑136 (MQ=255)                               gATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTgcg                                                                                                     <  1:938130/136‑1 (MQ=255)                               gATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTg                                                                                                 >  2:172196/1‑140 (MQ=255)                                    ggATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTgccgcc                                                                                            <  1:403493/140‑1 (MQ=255)                                         atatTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAAtt                                                                                       <  2:654408/140‑1 (MQ=255)                                               tCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGa                                                                                 >  2:770274/1‑140 (MQ=255)                                                   tGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAAtt                                                                             <  1:172196/140‑1 (MQ=255)                                                    gCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAATTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTc                                                                            >  1:1048802/1‑140 (MQ=255)                                                    gCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTc                                                                            >  1:1023880/1‑140 (MQ=255)                                                    gCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTc                                                                            >  1:568390/1‑140 (MQ=255)                                                     cGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCa                                                                           >  2:549297/1‑140 (MQ=255)                                                     cGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCa                                                                           <  2:1023880/140‑1 (MQ=255)                                                      gagaATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAg                                                                          <  1:295122/140‑1 (MQ=255)                                                      gagaATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAg                                                                          <  1:959144/140‑1 (MQ=255)                                                        gaATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGca                                                                        <  1:467675/140‑1 (MQ=255)                                                        gaATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGca                                                                        <  2:568390/140‑1 (MQ=255)                                                         aaTCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGcac                                                                       <  1:32101/140‑1 (MQ=255)                                                                   tCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTg                                                             >  2:795843/1‑140 (MQ=255)                                                                    cTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACAt                                                                     <  1:589058/131‑1 (MQ=255)                                                                    cTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACAt                                                                     >  2:589058/1‑131 (MQ=255)                                                                               gTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGa                                                 <  1:770274/140‑1 (MQ=255)                                                                                 tttGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGata                                               >  1:272301/1‑140 (MQ=255)                                                                                          aaaCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAAt                                      <  1:795843/140‑1 (MQ=255)                                                                                           aaCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATc                                     >  1:291244/1‑140 (MQ=255)                                                                                             cTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGa                                   <  1:978883/140‑1 (MQ=255)                                                                                             cTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGa                                   <  2:272301/140‑1 (MQ=255)                                                                                             cTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGa                                   <  2:358234/140‑1 (MQ=255)                                                                                              tGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGAc                                  <  2:291244/140‑1 (MQ=255)                                                                                                tttGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGa                                   <  1:923460/137‑1 (MQ=255)                                                                                                tttGATGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGa                                   >  2:923460/1‑137 (MQ=255)                                                                                                    aTGCTGAGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCata                            <  1:290361/140‑1 (MQ=255)                                                                                                          aGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATActg                      >  2:1070525/1‑140 (MQ=255)                                                                                                          aGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATActg                      >  2:758202/1‑140 (MQ=255)                                                                                                          aGCCAAATTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATActg                      >  1:703128/1‑140 (MQ=255)                                                                                                                aTTGAGATTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCa                <  1:291476/140‑1 (MQ=255)                                                                                                                     gaTTAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGt           <  1:758202/140‑1 (MQ=255)                                                                                                                        tAATGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAg        <  2:51454/140‑1 (MQ=255)                                                                                                                          aTGCCAGAGCG‑AAAATGCCGGTTCGAAATAAATTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCAAAATCGACTCCATAATACTGCTGCCACCGGTGAGGa      <  2:1048802/140‑1 (MQ=255)                                                                                                                           tGCCAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGAt     >  2:451704/1‑140 (MQ=255)                                                                                                                              cAGAGCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAAt  <  2:703128/140‑1 (MQ=255)                                                                                                                                gagCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCa                >  1:176454/1‑124 (MQ=255)                                                                                                                                gagCG‑AAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCa                <  2:176454/124‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      AACAGGATAAAGGTCATAACCCCTGAAAACGATCCGGATAATATTATCCCTGCGAGAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAAT  >  NC_000913/2130208‑2130473

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