Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,114,502	C→T	*465* (TAG→TAA)	wcaM ←	colanic acid biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,114,502	0	C	T	100.0%	221.9 / NA	68	*465* (TAG→TAA)	wcaM	colanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/33);  total (35/33)

TGTTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGT  >  NC_000913/2114367‑2114631                                                                                                                                        |                                                                                                                                  tGTTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTAc                                                                                                                                >  2:1235460/1‑139 (MQ=255)               cGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                                                  <  1:1235460/139‑1 (MQ=255)               cGCCAGCATCAATTACCCTGAAGTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                                                  <  1:787547/139‑1 (MQ=255)                               ccTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtc                                                                                                  >  1:653077/1‑139 (MQ=255)                                cTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtct                                                                                                 >  1:608796/1‑139 (MQ=255)                                 tGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtctt                                                                                                >  2:222579/1‑139 (MQ=255)                                   aaTTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtcttct                                                                                              >  1:50365/1‑139 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                     <  2:26989/127‑1 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                     >  1:26989/1‑127 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                     <  2:766217/127‑1 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                     >  1:766217/1‑127 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCcaca                                                                                         >  1:1385199/1‑139 (MQ=255)                                         ttACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCCCAt                                                                                        >  1:393432/1‑139 (MQ=255)                                         ttACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        <  2:393432/139‑1 (MQ=255)                                            cTGAATTACTTGTGATGTTACGGCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTCTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                                     <  2:280764/139‑1 (MQ=255)                                             tGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAAc                                                                                    >  1:1484210/1‑139 (MQ=255)                                               aaTTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACaa                                                                                  >  2:532164/1‑139 (MQ=255)                                                   aCTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                              <  1:222579/139‑1 (MQ=255)                                                   aCTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                              <  2:50365/139‑1 (MQ=255)                                                       gtgATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGa                                                                          >  1:141176/1‑139 (MQ=255)                                                        tgATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGAt                                                                         <  1:532164/139‑1 (MQ=255)                                                          aTGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTa                                                                                                                                 >  1:2398/1‑81 (MQ=255)                                                          aTGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTa                                                                                                                                 <  2:2398/81‑1 (MQ=255)                                                              tACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTaa                                                                                                          >  1:1110448/1‑100 (MQ=255)                                                              tACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTaa                                                                                                          <  2:1110448/100‑1 (MQ=255)                                                               aCGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCt                                                                  >  1:864296/1‑139 (MQ=255)                                                                     cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGat                                                            >  1:1135344/1‑139 (MQ=255)                                                                       tCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             >  2:908998/1‑104 (MQ=255)                                                                       tCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             <  1:908998/104‑1 (MQ=255)                                                                          ttGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCg                                                                                                                   >  2:228075/1‑79 (MQ=255)                                                                          ttGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCg                                                                                                                   <  1:228075/79‑1 (MQ=255)                                                                                 ttGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCt                                                <  2:1484210/139‑1 (MQ=255)                                                                                  tGCAGTAGCATTGTTCCTAAGTATGACCCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTa                                                                                                           <  1:127358/79‑1 (MQ=255)                                                                                  tGCAGTAGCATTGTTCCTAAGTATGACCCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTa                                                                                                           >  2:127358/1‑79 (MQ=255)                                                                                    cAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                                     >  1:1042638/1‑99 (MQ=255)                                                                                    cAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                                     <  2:1042638/99‑1 (MQ=255)                                                                                         gCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                                           >  2:621684/1‑56 (MQ=255)                                                                                         gCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                                           <  1:621684/56‑1 (MQ=255)                                                                                            ttGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCttttt                                     >  2:738926/1‑139 (MQ=255)                                                                                              gTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                                           <  2:393728/51‑1 (MQ=255)                                                                                              gTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                                           >  1:393728/1‑51 (MQ=255)                                                                                                 ccTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcattc                                >  1:845817/1‑139 (MQ=255)                                                                                                 ccTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcattc                                <  2:653077/139‑1 (MQ=255)                                                                                                       tATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        <  1:1348344/77‑1 (MQ=255)                                                                                                       tATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        >  2:1348344/1‑77 (MQ=255)                                                                                                        aTGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCAt                         <  1:1312105/139‑1 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        <  2:895284/75‑1 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        <  2:263497/75‑1 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        >  1:895284/1‑75 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        >  1:263497/1‑75 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATg                        <  2:1385199/139‑1 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATg                        >  1:155266/1‑139 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGCGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATg                        >  2:714932/1‑139 (MQ=255)                                                                                                          gACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGa                       <  1:738926/139‑1 (MQ=255)                                                                                                          gACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGa                       <  2:864296/139‑1 (MQ=255)                                                                                                               cATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                              <  2:635867/79‑1 (MQ=255)                                                                                                               cATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                              >  1:635867/1‑79 (MQ=255)                                                                                                                aTTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTcc                                                                                            <  1:885586/64‑1 (MQ=255)                                                                                                                aTTTTTCCAGGAATGGTCGCAAATTTACTACCTCAGTTCCGGCAATCTAAAGTTAATCTTCTcc                                                                                            >  2:885586/1‑64 (MQ=255)                                                                                                                 tttttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTg                <  2:155266/139‑1 (MQ=255)                                                                                                                 tttttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTg                >  1:1494377/1‑139 (MQ=255)                                                                                                                   tttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCa              <  1:714932/139‑1 (MQ=255)                                                                                                                      ccAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAaga           <  2:1494377/139‑1 (MQ=255)                                                                                                                       cAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                     >  1:42912/1‑48 (MQ=255)                                                                                                                       cAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                     <  2:42912/48‑1 (MQ=255)                                                                                                                         ggAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACa                                                                                   <  1:945766/64‑1 (MQ=255)                                                                                                                         ggAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACa                                                                                   >  2:945766/1‑64 (MQ=255)                                                                                                                               gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAgt  >  1:840061/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                  TGTTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGT  >  NC_000913/2114367‑2114631

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