Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,132,499	A→G	L19L (CTT→CTC)	wcaB ←	putative colanic acid biosynthesis acyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,132,499	0	A	G	100.0%	218.8 / NA	61	L19L (CTT→CTC)	wcaB	putative colanic acid biosynthesis acyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/30);  total (31/30)

CGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTC  >  NC_000913/2132364‑2132633                                                                                                                                        |                                                                                                                                       cgcgGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGaa                                                                                                                                    <  2:392642/140‑1 (MQ=255) cgcgGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCCAACAAAAATGAGCAACACGATAGGCGAGaa                                                                                                                                    <  2:1164636/140‑1 (MQ=255)   cgGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAAcc                                                                                                                                  >  1:706019/1‑140 (MQ=255)   cgGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAAcc                                                                                                                                  >  2:1293165/1‑140 (MQ=255)   cgGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAAcc                                                                                                                                  >  2:365462/1‑140 (MQ=255)           tGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATgcagca                                                                                                                          >  2:1237868/1‑140 (MQ=255)                 tcgtAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCg                                                                                                                    <  1:541695/137‑1 (MQ=255)                  cATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGt                                                                                                                   >  1:1344194/1‑140 (MQ=255)                          aaaaaGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCa                                                                                                           >  2:1013100/1‑140 (MQ=255)                                    cGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGgcg                                                                                                 >  2:213266/1‑140 (MQ=255)                                     ggTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATATGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGgcgc                                                                                                <  2:1428278/140‑1 (MQ=255)                                         aTAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAg                                                                                            >  2:350061/1‑140 (MQ=255)                                              gCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCtt                                                                                       <  2:606761/140‑1 (MQ=255)                                              gCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCtt                                                                                       <  1:365462/140‑1 (MQ=255)                                              gCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCtt                                                                                       <  1:445601/140‑1 (MQ=255)                                               cGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTc                                                                                      >  2:246343/1‑140 (MQ=255)                                               cGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTc                                                                                      >  1:223278/1‑140 (MQ=255)                                                  tACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAAATCTTCCAg                                                                                   <  2:710093/140‑1 (MQ=255)                                                   aCAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGc                                                                                  <  2:223278/140‑1 (MQ=255)                                                   aCAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGc                                                                                  <  1:919425/140‑1 (MQ=255)                                                      gCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATc                                                                               >  1:1334470/1‑140 (MQ=255)                                                       cACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCt                                                                              >  1:1465419/1‑140 (MQ=255)                                                       cACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCt                                                                              >  1:675216/1‑140 (MQ=255)                                                         ccAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCTACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCttt                                                                            >  2:382715/1‑140 (MQ=255)                                                         ccAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCttt                                                                            >  1:292515/1‑140 (MQ=255)                                                          cagcagCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTa                                                                           >  2:1390351/1‑140 (MQ=255)                                                             cagcGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAgcg                                                                        <  1:213266/140‑1 (MQ=255)                                                             cagcGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAgcg                                                                        <  1:1477366/140‑1 (MQ=255)                                                               gcGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGcc                                                                      >  2:779890/1‑140 (MQ=255)                                                                 ggggCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGcccc                                                                    >  1:1250178/1‑140 (MQ=255)                                                                    gccgccCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGa                                                                 <  1:902638/140‑1 (MQ=255)                                                                     ccgccCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGaa                                                                <  1:1390351/140‑1 (MQ=255)                                                                               aTTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGc                                                      >  2:515831/1‑140 (MQ=255)                                                                                 tgttgAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCa                                                    <  1:1013100/140‑1 (MQ=255)                                                                                   ttgAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGa                                                  >  1:663858/1‑140 (MQ=255)                                                                                       ggACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGtt                                              <  1:692918/140‑1 (MQ=255)                                                                                       ggACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGtt                                              >  2:561591/1‑140 (MQ=255)                                                                                       ggACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGtt                                              >  2:824914/1‑140 (MQ=255)                                                                                         aCGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTg                                            <  2:972912/140‑1 (MQ=255)                                                                                         aCGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTg                                            >  1:897944/1‑140 (MQ=255)                                                                                           gTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGcc                                          >  2:760590/1‑140 (MQ=255)                                                                                             tCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGt                                        <  1:1293165/140‑1 (MQ=255)                                                                                             tCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGt                                        >  2:1068707/1‑140 (MQ=255)                                                                                                ttttGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTAc                                     <  1:1068707/140‑1 (MQ=255)                                                                                                ttttGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTAc                                     <  1:879758/140‑1 (MQ=255)                                                                                                   tGCGCCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCa                                  >  1:137882/1‑140 (MQ=255)                                                                                                    gcgcCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCAc                                 <  2:33787/140‑1 (MQ=255)                                                                                                     cgcCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCAcc                                >  2:1396021/1‑140 (MQ=255)                                                                                                     cgcCACACCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCAcc                                <  2:675216/140‑1 (MQ=255)                                                                                                         acacCGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGAca                            >  1:1281904/1‑140 (MQ=255)                                                                                                            ccGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGAcagca                         <  1:515831/140‑1 (MQ=255)                                                                                                             cGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGAcagcag                        >  1:1498196/1‑140 (MQ=255)                                                                                                             cGAACAAAAATGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGAcagcag                        <  2:1311310/140‑1 (MQ=255)                                                                                                                     aaTGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCa                >  1:1541118/1‑140 (MQ=255)                                                                                                                       tGAGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCaaa              <  2:1281904/140‑1 (MQ=255)                                                                                                                         aGCAACACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGt            <  1:561591/140‑1 (MQ=255)                                                                                                                            aaCACGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTgcg         <  2:1541118/140‑1 (MQ=255)                                                                                                                             acacGATAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTgcgc        <  1:1252663/140‑1 (MQ=255)                                                                                                                                  aTAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGt   <  1:350061/140‑1 (MQ=255)                                                                                                                                  aTAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGt   <  2:137882/140‑1 (MQ=255)                                                                                                                                   tAGGCGAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTc  >  2:915716/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       CGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTC  >  NC_000913/2132364‑2132633

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: ‑