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%	209.1 / NA	65	*465* (TAG→TAA)	wcaM	colanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/34);  total (31/34)

TTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGT  >  NC_000913/2114369‑2114631                                                                                                                                      |                                                                                                                                  ttCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTc                                                                                                                              >  2:816763/1‑139 (MQ=255) ttCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTc                                                                                                                              >  1:1214248/1‑139 (MQ=255) ttCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTc                                                                                                                              >  1:167790/1‑139 (MQ=255)   cGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTccc                                                                                                                            <  2:167790/139‑1 (MQ=255)                   cATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCt                                                                                                            <  2:757703/139‑1 (MQ=255)                     tCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTaa                                                                                                          <  1:1248392/139‑1 (MQ=255)                                    tGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCca                                                                                           <  2:1214248/139‑1 (MQ=255)                                     gATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCcac                                                                                          <  2:998553/139‑1 (MQ=255)                                        tACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAtt                                                                                       <  1:816763/139‑1 (MQ=255)                                           tGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAAc                                                                                    >  2:1006983/1‑139 (MQ=255)                                             aaTTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACaa                                                                                  >  2:1258566/1‑139 (MQ=255)                                                 aCTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                              <  1:1195353/139‑1 (MQ=255)                                                 aCTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATCTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                              <  1:778648/139‑1 (MQ=255)                                                                   cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTccc                                                                                                                            <  1:336900/75‑1 (MQ=255)                                                                   cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTccc                                                                                                                            >  2:336900/1‑75 (MQ=255)                                                                    ttCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                     >  2:471385/1‑97 (MQ=255)                                                                    ttCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                     <  1:471385/97‑1 (MQ=255)                                                                                tGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCtt                                               >  1:552569/1‑139 (MQ=255)                                                                                  cAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                                           <  2:1170528/61‑1 (MQ=255)                                                                                  cAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                                           >  1:1170528/1‑61 (MQ=255)                                                                                  cAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTg                                             >  1:1078662/1‑139 (MQ=255)                                                                                    gTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTc                                           <  1:1258566/139‑1 (MQ=255)                                                                                       gCATTGTTCCT‑AGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        >  1:251725/1‑90 (MQ=255)                                                                                       gCATTGTTCCT‑AGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        <  2:251725/90‑1 (MQ=255)                                                                                               ccTAAGTGTGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTc                                                               <  2:230236/108‑1 (MQ=255)                                                                                               ccTAAGTGTGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTc                                                               >  1:230236/1‑108 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTccc                                                                                                                            >  1:466076/1‑39 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTccc                                                                                                                            <  2:466076/39‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAAt                                                                                                              <  1:302175/53‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAAt                                                                                                              >  2:302175/1‑53 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGt                                                                                                       <  1:774901/60‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGt                                                                                                       >  2:774901/1‑60 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                                     >  2:862189/1‑78 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                                     <  1:862189/78‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAAt                                                                     <  2:181062/94‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAAt                                                                     >  1:181062/1‑94 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACg                                                     <  1:159215/110‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACg                                                     >  2:159215/1‑110 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcattca                               >  1:1260556/1‑132 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATg                        <  2:552569/139‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATg                        >  1:527429/1‑139 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATATAAAGTTAATCTTCTCCACATTAGCAATATGGTGATTAATCCTGTCGATATCGACCGAGCTTTGTCCTTTttcattca                               <  2:1260556/132‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCACAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                                     <  2:965745/78‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCACAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                                     >  1:965745/1‑78 (MQ=255)                                                                                                        gACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGa                       <  2:1078662/139‑1 (MQ=255)                                                                                                         aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                                           <  1:828537/38‑1 (MQ=255)                                                                                                         aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                                           >  2:828537/1‑38 (MQ=255)                                                                                                         aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                                                  >  2:96045/1‑47 (MQ=255)                                                                                                         aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                                                  <  1:96045/47‑1 (MQ=255)                                                                                                            ccATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACAt                   <  2:920972/139‑1 (MQ=255)                                                                                                               tttttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                     >  2:468991/1‑54 (MQ=255)                                                                                                               tttttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                     <  1:468991/54‑1 (MQ=255)                                                                                                               tttttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCGCATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTg                >  2:858013/1‑139 (MQ=255)                                                                                                               tttttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        >  2:626621/1‑67 (MQ=255)                                                                                                               tttttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        <  1:626621/67‑1 (MQ=255)                                                                                                                 tttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCa              <  1:640132/139‑1 (MQ=255)                                                                                                                     cAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAagag          >  2:141198/1‑139 (MQ=255)                                                                                                                       ggAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGAc        <  1:141198/139‑1 (MQ=255)                                                                                                                            ggTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             >  1:568477/1‑49 (MQ=255)                                                                                                                            ggTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             <  2:568477/49‑1 (MQ=255)                                                                                                                             gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             >  2:1052498/1‑48 (MQ=255)                                                                                                                             gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             <  1:1052498/48‑1 (MQ=255)                                                                                                                             gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAgt  >  1:1026536/1‑139 (MQ=255)                                                                                                                             gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAAT‑CTGTCGATATCGACGGAGCTTTGTCCTTTttcattca                               <  1:890312/109‑1 (MQ=255)                                                                                                                             gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAAT‑CTGTCGATATCGACGGAGCTTTGTCCTTTttcattca                               >  2:890312/1‑109 (MQ=255)                                                                                                                                      |                                                                                                                                  TTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGT  >  NC_000913/2114369‑2114631

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