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_0009132,130,8530CT100.0% 98.3 / NA 31S9N (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 (15/16);  total (15/16)

GCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGT  >  NC_000913/2130732‑2130973
                                                                                                                         |                                                                                                                        
gCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGaa                                                                                                         <  2:746199/139‑1 (MQ=255)
       aaaaaCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACa                                                                                                  <  1:14196/139‑1 (MQ=255)
                agTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACtcctcc                                                                                         >  1:594742/1‑139 (MQ=255)
                  tAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAg                                                                                       >  2:473674/1‑139 (MQ=255)
                          cAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgtt                                                                               >  2:629979/1‑139 (MQ=255)
                             gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACa                                                                                                  <  2:556486/117‑1 (MQ=255)
                             gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACa                                                                                                  >  1:556486/1‑117 (MQ=255)
                             gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTc                                                                            >  2:390543/1‑139 (MQ=255)
                                      tGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                          <  1:344703/100‑1 (MQ=255)
                                      tGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                          >  2:344703/1‑100 (MQ=255)
                                      tGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                          <  1:229793/100‑1 (MQ=255)
                                      tGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                          >  2:229793/1‑100 (MQ=255)
                                        aCGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATa                                                                                                >  1:588930/1‑108 (MQ=255)
                                        aCGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATa                                                                                                <  2:588930/108‑1 (MQ=255)
                                               tGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGgcggcgcggc                                                          <  1:629979/139‑1 (MQ=255)
                                               tGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGgcggcgcggc                                                          <  1:390543/139‑1 (MQ=255)
                                                     ggAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                          <  2:189132/85‑1 (MQ=255)
                                                     ggAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                          >  1:189132/1‑85 (MQ=255)
                                                      gAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTg                                                   >  1:5546/1‑139 (MQ=255)
                                                        aaCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                          <  1:505393/82‑1 (MQ=255)
                                                        aaCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                          >  2:505393/1‑82 (MQ=255)
                                                              cGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACt                                           >  1:199232/1‑139 (MQ=255)
                                                                     cTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTcagccag                                    <  2:199232/139‑1 (MQ=255)
                                                                       gggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                  <  2:5546/139‑1 (MQ=255)
                                                                                                aTAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGc                                                 >  2:145546/1‑99 (MQ=255)
                                                                                                aTAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGc                                                 <  1:145546/99‑1 (MQ=255)
                                                                                                  aaaTCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                                                                    >  2:77218/1‑62 (MQ=255)
                                                                                                  aaaTCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                                                                    <  1:77218/62‑1 (MQ=255)
                                                                                                       aGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGt  <  1:473674/139‑1 (MQ=255)
                                                                                                          ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCACCt                                                                                  >  2:298145/1‑56 (MQ=255)
                                                                                                          ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCACCt                                                                                  <  1:298145/56‑1 (MQ=255)
                                                                                                                         |                                                                                                                        
GCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGT  >  NC_000913/2130732‑2130973

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