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%	97.9 / NA	31	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 (12/19);  total (12/19)

TTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTGCTTAAC  >  NC_000913/2130717‑2130982                                                                                                                                         |                                                                                                                                  ttttttAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTa                                                                                                                                 <  1:115093/139‑1 (MQ=255)            tGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACAATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGAt                                                                                                                      <  2:130815/139‑1 (MQ=255)               tGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                                   <  1:386171/139‑1 (MQ=255)               tGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                                   <  2:670236/139‑1 (MQ=255)                                agTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGAtt                                                                                                                           >  2:780820/1‑114 (MQ=255)                                agTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGAtt                                                                                                                           <  1:780820/114‑1 (MQ=255)                                   aaCAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGc                                                                                               <  2:947499/139‑1 (MQ=255)                                     cAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                                                                             <  1:892035/139‑1 (MQ=255)                                        gACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtg                                                                                          <  1:38280/139‑1 (MQ=255)                                            ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                                   >  2:297691/1‑110 (MQ=255)                                            ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                                   <  1:297691/110‑1 (MQ=255)                                            ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACa                                                                                                           >  1:152647/1‑118 (MQ=255)                                            ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACa                                                                                                           <  2:152647/118‑1 (MQ=255)                                            ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGc                                                                                               <  1:417726/130‑1 (MQ=255)                                            ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGc                                                                                               >  2:417726/1‑130 (MQ=255)                                            ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtg                                                                                          >  2:952686/1‑135 (MQ=255)                                            ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtg                                                                                          <  1:952686/135‑1 (MQ=255)                                                                       aaaCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTgg                                                           >  2:718905/1‑139 (MQ=255)                                                                             tCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACa                                                                                                           >  2:898668/1‑85 (MQ=255)                                                                             tCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACa                                                                                                           <  1:898668/85‑1 (MQ=255)                                                                             tCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                                                                         <  2:681519/119‑1 (MQ=255)                                                                             tCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                                                                         >  1:681519/1‑119 (MQ=255)                                                                                   aGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                                                                             >  1:732891/1‑93 (MQ=255)                                                                                   aGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                                                                             <  2:732891/93‑1 (MQ=255)                                                                                       gggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                           <  1:718905/139‑1 (MQ=255)                                                                                         gCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTg                                         <  1:162502/139‑1 (MQ=255)                                                                                                                  aaaTCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGAttt                >  2:204293/1‑139 (MQ=255)                                                                                                                  aaaTCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGAttt                >  1:130500/1‑139 (MQ=255)                                                                                                                         cGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTa                                                                              <  2:436654/70‑1 (MQ=255)                                                                                                                         cGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTa                                                                              >  1:436654/1‑70 (MQ=255)                                                                                                                                aaCAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTGCTTAAc  <  1:204293/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  TTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTGCTTAAC  >  NC_000913/2130717‑2130982

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