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%	168.0 / NA	52	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 (25/27);  total (25/27)

TTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTG  >  NC_000913/2130719‑2130976                                                                                                                                       |                                                                                                                            ttttAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACa                                                                                                                         <  2:1217340/139‑1 (MQ=255)         cTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGa                                                                                                                 >  1:387882/1‑139 (MQ=255)               cGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAg                                                                                                           <  1:1448902/139‑1 (MQ=255)                  tCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTg                                                                                                        >  1:256135/1‑139 (MQ=255)                       aaaCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATa                                                                                                   <  2:387882/139‑1 (MQ=255)                           agcagTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATActcc                                                                                               >  2:163982/1‑139 (MQ=255)                              agTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACtcctcc                                                                                            >  2:282769/1‑139 (MQ=255)                                     aGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCt                                                                                     <  1:163982/139‑1 (MQ=255)                                          ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTg                                                                                                        <  1:1201214/115‑1 (MQ=255)                                          ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTg                                                                                                        >  2:1201214/1‑115 (MQ=255)                                           gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATActcct                                                                                              >  2:323409/1‑124 (MQ=255)                                           gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATActcct                                                                                              <  1:323409/124‑1 (MQ=255)                                                           ggTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGgcggcgcg                                                               >  1:1358981/1‑139 (MQ=255)                                                                  gggAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATa                                                                                                   >  1:895152/1‑96 (MQ=255)                                                                  gggAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATa                                                                                                   <  2:895152/96‑1 (MQ=255)                                                                      aaCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATActcct                                                                                              <  1:335771/97‑1 (MQ=255)                                                                      aaCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATActcct                                                                                              >  2:335771/1‑97 (MQ=255)                                                                      aaCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGc                                                    >  1:1444759/1‑139 (MQ=255)                                                                      aaCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGc                                                    <  2:1444759/139‑1 (MQ=255)                                                                       aCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCt                                                   >  1:702734/1‑139 (MQ=255)                                                                           tCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcg                                                                    <  2:381500/118‑1 (MQ=255)                                                                           tCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCCGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcg                                                                    >  1:381500/1‑118 (MQ=255)                                                                             gCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTc                                             >  1:890909/1‑139 (MQ=255)                                                                              ccAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTca                                            >  2:918382/1‑139 (MQ=255)                                                                                     gggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                     <  2:1182398/139‑1 (MQ=255)                                                                                     gggCGAATTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACACAGCCAGCg                                     <  2:359231/139‑1 (MQ=255)                                                                                      ggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGt                                    <  2:1358981/139‑1 (MQ=255)                                                                                           ttttGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTAc                               <  2:256135/139‑1 (MQ=255)                                                                                                     aaCCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                             >  2:337255/1‑51 (MQ=255)                                                                                                     aaCCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                             <  1:337255/51‑1 (MQ=255)                                                                                                     aaCCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                                                                   >  2:34990/1‑93 (MQ=255)                                                                                                     aaCCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                                                                   <  1:34990/93‑1 (MQ=255)                                                                                                     aaCCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGc                                                    >  1:1142984/1‑108 (MQ=255)                                                                                                     aaCCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGc                                                    <  2:1142984/108‑1 (MQ=255)                                                                                                          gcagATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAg                                                                                                           >  2:108400/1‑48 (MQ=255)                                                                                                          gcagATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAg                                                                                                           <  1:108400/48‑1 (MQ=255)                                                                                                                aaaTCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCa                                                                             <  1:911988/72‑1 (MQ=255)                                                                                                                aaaTCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCa                                                                             >  2:911988/1‑72 (MQ=255)                                                                                                                   tCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCg       <  1:282769/139‑1 (MQ=255)                                                                                                                     aGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCg       >  2:1330970/1‑137 (MQ=255)                                                                                                                     aGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCg       <  1:1330970/137‑1 (MQ=255)                                                                                                                       cGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTcagc                                          >  2:1263370/1‑100 (MQ=255)                                                                                                                       cGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTcagc                                          <  1:1263370/100‑1 (MQ=255)                                                                                                                        ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                     >  2:518089/1‑104 (MQ=255)                                                                                                                        ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTg  <  2:702734/139‑1 (MQ=255)                                                                                                                        ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCAAGCg                                     <  1:518089/104‑1 (MQ=255)                                                                                                                          cagcaACAGGTAGTCACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTgcg                >  1:455902/1‑123 (MQ=255)                                                                                                                          cagcaACAGGTAGTCACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTgcg                <  2:455902/123‑1 (MQ=255)                                                                                                                            gcaACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCt                                                                                     >  1:922198/1‑52 (MQ=255)                                                                                                                            gcaACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCt                                                                                     <  2:922198/52‑1 (MQ=255)                                                                                                                            gcaACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTa                                >  2:743470/1‑105 (MQ=255)                                                                                                                            gcaACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTa                                <  1:743470/105‑1 (MQ=255)                                                                                                                                       |                                                                                                                            TTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTG  >  NC_000913/2130719‑2130976

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: ‑