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

TTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAG  >  NC_000913/2114369‑2114627                                                                                                                                      |                                                                                                                              ttCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTc                                                                                                                          >  2:481570/1‑139 (MQ=255)                          tACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                 >  2:532398/1‑139 (MQ=255)                             ccTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtc                                                                                              >  2:171945/1‑139 (MQ=255)                              cTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtct                                                                                             >  1:422772/1‑139 (MQ=255)                                gAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtcttc                                                                                           <  1:532398/139‑1 (MQ=255)                                 aaTTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtcttct                                                                                          >  1:95890/1‑139 (MQ=255)                                      aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCcaca                                                                                     >  2:792435/1‑139 (MQ=255)                                       ttACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                    <  1:481570/139‑1 (MQ=255)                                        tACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTCCTCCACAtt                                                                                   <  1:2002/139‑1 (MQ=255)                                           tGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCca                                                                                       <  2:566893/132‑1 (MQ=255)                                           tGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAAc                                                                                >  2:242231/1‑139 (MQ=255)                                           tGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAACATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCca                                                                                       >  1:566893/1‑132 (MQ=255)                                               ttACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAata                                                                            >  1:503285/1‑139 (MQ=255)                                                 aCTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                          <  1:801051/139‑1 (MQ=255)                                                 aCTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                          <  2:503285/139‑1 (MQ=255)                                                 aCTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                          <  1:147315/139‑1 (MQ=255)                                                     gtgATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGa                                                                      <  2:490645/139‑1 (MQ=255)                                                       gATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGAtt                                                                    >  1:214458/1‑139 (MQ=255)                                                           ttACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                         <  2:368906/114‑1 (MQ=255)                                                           ttACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                         >  1:368906/1‑114 (MQ=255)                                                             aCGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATAT‑GTGATTAATCCTg                                                             >  2:336669/1‑139 (MQ=255)                                                                   cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTa                                                                                                                             <  1:703179/70‑1 (MQ=255)                                                                   cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTa                                                                                                                             >  2:703179/1‑70 (MQ=255)                                                                    ttCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtcttc                                                                                           <  2:796004/103‑1 (MQ=255)                                                                    ttCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtcttc                                                                                           >  1:796004/1‑103 (MQ=255)                                                                        ttGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGa                                                   >  2:213537/1‑139 (MQ=255)                                                                             gATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTcc                                                                                                                         >  2:364899/1‑64 (MQ=255)                                                                             gATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTcc                                                                                                                         <  1:364899/64‑1 (MQ=255)                                                                                  cAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                                 >  2:59014/1‑99 (MQ=255)                                                                                  cAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                                 <  1:59014/99‑1 (MQ=255)                                                                                   aGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCg                                                                                                               >  2:517429/1‑68 (MQ=255)                                                                                   aGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCg                                                                                                               <  1:517429/68‑1 (MQ=255)                                                                                     tAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCg                                                                                                               >  2:727725/1‑66 (MQ=255)                                                                                     tAGCATTGTTCCTAAGTATGACCCCATTTTTCCAGGAATGGTCGAAAATTTACTCCCTCAGTTCCg                                                                                                               <  1:727725/66‑1 (MQ=255)                                                                                       gCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAatat                                                                           <  2:402973/100‑1 (MQ=255)                                                                                       gCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAatat                                                                           >  1:402973/1‑100 (MQ=255)                                                                                          ttGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                          >  2:236571/1‑98 (MQ=255)                                                                                          ttGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                          <  1:236571/98‑1 (MQ=255)                                                                                            gTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttca                               >  2:57413/1‑139 (MQ=255)                                                                                             ttCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATAT‑GTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcatt                             <  1:336669/139‑1 (MQ=255)                                                                                                cTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcattca                           <  2:753965/139‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACa                                                                               <  1:620497/80‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACa                                                                               >  2:620497/1‑80 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcattc                            <  2:209485/131‑1 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcattc                            >  1:209485/1‑131 (MQ=255)                                                                                                       tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATg                    <  1:792435/139‑1 (MQ=255)                                                                                                         aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGaa                  >  1:283320/1‑139 (MQ=255)                                                                                                         aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGACGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGaa                  >  1:667216/1‑139 (MQ=255)                                                                                                            ccATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACAt               <  2:95890/139‑1 (MQ=255)                                                                                                             cATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcattca                           >  2:782340/1‑126 (MQ=255)                                                                                                             cATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcattca                           <  1:782340/126‑1 (MQ=255)                                                                                                               tttttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                    >  2:71206/1‑67 (MQ=255)                                                                                                               tttttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                    <  1:71206/67‑1 (MQ=255)                                                                                                                 tttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCa          <  1:213537/139‑1 (MQ=255)                                                                                                                         aaTGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGAcag  <  1:242231/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                              TTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAG  >  NC_000913/2114369‑2114627

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