Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,130,853	C→T	S9N (AGC→AAC)  *406* (TAG→TAA)	wcaD ← wcaC ←	putative colanic acid polymerase putative glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,130,853	0	C	T	100.0%	156.7 / NA	49	S9N (AGC→AAC)  *406* (TAG→TAA)	wcaD wcaC	putative colanic acid polymerase putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/27);  total (22/27)

AACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTGC  >  NC_000913/2130723‑2130977                                                                                                                                   |                                                                                                                             aaCGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGAtt                                                                                                                      <  2:607620/139‑1 (MQ=255) aaCGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGAtt                                                                                                                      <  1:542020/139‑1 (MQ=255)            gTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACTGGTAGTTACAGATTCTGATAGAAGt                                                                                                           <  2:613782/139‑1 (MQ=255)                        gcagTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATActcct                                                                                               >  1:353372/1‑139 (MQ=255)                                aaGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGA‑TTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCt                                                                                      <  1:90899/139‑1 (MQ=255)                                  gACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtg                                                                                     <  2:353372/139‑1 (MQ=255)                                      ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt                                                                                 <  2:1173644/139‑1 (MQ=255)                                       gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAg                                                                                                            <  1:152202/111‑1 (MQ=255)                                       gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAg                                                                                                            >  2:152202/1‑111 (MQ=255)                                        tAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACCCCTCCAGc                                                                                          <  1:286320/128‑1 (MQ=255)                                        tAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGAAATACCCCTCCAGc                                                                                          >  2:286320/1‑128 (MQ=255)                                                 gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATAc                                                                                                   <  2:127079/110‑1 (MQ=255)                                                 gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATAc                                                                                                   >  1:127079/1‑110 (MQ=255)                                                 gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                                                                        >  1:1174658/1‑121 (MQ=255)                                                 gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                                                                        <  2:1174658/121‑1 (MQ=255)                                                  aCGATGGTGCGGGGGAAACTTTCCCCATGCTGCGCGATCTTGACGTTAACCAGCAGATCAATCAGCGGCAGCAGCAGGTAGTAAACGACTCAGACAGCAGTTGACCAACTCCTCCAGCATCTCTTGTCCACTGTAggcg                                                                     >  2:1226619/1‑139 (MQ=255)                                                  aCGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcg                                                                     >  2:1109266/1‑139 (MQ=255)                                                       ggTGATGGGGGAACTTCCGGAAAGCTGGGCGATTTTGATGTTAATACGCAGATAAATCAGGGACAGCAACAGGTAGTTACAGCTTCTGATAGAAGTTGGCATACTCCTCGAGCATATGTTGTGTACTGTAggcggatcg                                                                >  2:1397549/1‑135 (MQ=255)                                                        gTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGGAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGTTAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGgcggcgcgg                                                               >  1:977210/1‑139 (MQ=255)                                                              gggAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGTGCGGCTGCGt                                                         >  1:1245005/1‑139 (MQ=255)                                                                 aaaCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGTGCGGCTGCGTTgg                                                      <  2:1245005/139‑1 (MQ=255)                                                                  aaCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                                                                        >  2:1455716/1‑104 (MQ=255)                                                                  aaCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                                                                        <  1:1455716/104‑1 (MQ=255)                                                                   aCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                                                                    >  1:945901/1‑123 (MQ=255)                                                                   aCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                                                                    <  2:945901/123‑1 (MQ=255)                                                                   aCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCt                                                    >  2:1449567/1‑139 (MQ=255)                                                                      ttCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGaa                                                 <  1:1449567/139‑1 (MQ=255)                                                                       tCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAAc                                                <  2:1463747/139‑1 (MQ=255)                                                                               cTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTcagccag                                        <  1:1109266/139‑1 (MQ=255)                                                                                 gggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                      <  2:728526/139‑1 (MQ=255)                                                                                  ggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACAt                                                                                                     >  2:1381645/1‑75 (MQ=255)                                                                                  ggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACAt                                                                                                     <  1:1381645/75‑1 (MQ=255)                                                                                            aCGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGATGGCTGAACTCAGCCAGCGTGGTACCAAat                           >  2:472736/1‑139 (MQ=255)                                                                                                       cagATAAATCAGCGGCAGCAACAGGTAGTTACAGAt                                                                                                                       >  2:95754/1‑36 (MQ=255)                                                                                                       cagATAAATCAGCGGCAGCAACAGGTAGTTACAGAt                                                                                                                       <  1:95754/36‑1 (MQ=37)                                                                                                       cagATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGc                                                                                          <  1:530627/65‑1 (MQ=255)                                                                                                       cagATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGc                                                                                          >  2:530627/1‑65 (MQ=255)                                                                                                          aTAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgtt                                                                                   <  1:95001/69‑1 (MQ=255)                                                                                                          aTAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgtt                                                                                   >  2:95001/1‑69 (MQ=255)                                                                                                                 aGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACa                                                                                                      <  1:18086/43‑1 (MQ=255)                                                                                                                 aGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACa                                                                                                      >  2:18086/1‑43 (MQ=255)                                                                                                                 aGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTc                                                                                <  2:695870/65‑1 (MQ=255)                                                                                                                 aGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTc                                                                                >  1:695870/1‑65 (MQ=255)                                                                                                                   cGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGgcggcgcgg                                                               <  2:1059964/80‑1 (MQ=255)                                                                                                                   cGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGgcggcgcgg                                                               >  1:1059964/1‑80 (MQ=255)                                                                                                                   cGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGc                                                     <  1:1374442/90‑1 (MQ=255)                                                                                                                   cGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGc                                                     >  2:1374442/1‑90 (MQ=255)                                                                                                                     gcagcaACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTGc  <  1:436891/139‑1 (MQ=255)                                                                                                                         caACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGc                                                           >  1:1462579/1‑78 (MQ=255)                                                                                                                         caACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGc                                                           <  2:1462579/78‑1 (MQ=255)                                                                                                                                   |                                                                                                                             AACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTGC  >  NC_000913/2130723‑2130977

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