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% 88.7 / NA 28S9N (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/16);  total (12/16)

AAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTT  >  NC_000913/2130739‑2130975
                                                                                                                  |                                                                                                                          
aaaaaCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACa                                                                                                    <  1:569442/139‑1 (MQ=255)
       gcagTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATActcct                                                                                             >  1:303168/1‑139 (MQ=255)
       gcagTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATActcct                                                                                             <  2:1121439/139‑1 (MQ=255)
         agTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACtcctcc                                                                                           <  1:378632/139‑1 (MQ=255)
              cAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                                                                      <  2:606442/139‑1 (MQ=255)
                aGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCt                                                                                    >  2:256716/1‑139 (MQ=255)
                aGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCt                                                                                    <  2:1073181/139‑1 (MQ=255)
                 gACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtg                                                                                   <  1:256716/139‑1 (MQ=255)
                     ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt                                                                               <  2:6856/139‑1 (MQ=255)
                     ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt                                                                               >  1:6856/1‑139 (MQ=255)
                     ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt                                                                               >  1:728112/1‑139 (MQ=255)
                                  cGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                                                                  <  1:889105/139‑1 (MQ=255)
                                        tGATGGGGAAACTTTCGCCAAGCTGGGCTATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGgcggcgcggc                                                            <  2:303168/139‑1 (MQ=255)
                                             gggAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGt                                                       >  1:369939/1‑139 (MQ=255)
                                                  aCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtg                                                                                   <  1:29082/106‑1 (MQ=255)
                                                  aCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtg                                                                                   >  2:29082/1‑106 (MQ=255)
                                                      tCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAAc                                              >  1:869625/1‑139 (MQ=255)
                                                                   cGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                            >  1:701132/1‑64 (MQ=255)
                                                                   cGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                            <  2:701132/64‑1 (MQ=255)
                                                                                aaCCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGc                    <  2:369939/139‑1 (MQ=255)
                                                                                              tCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCg      <  1:455137/139‑1 (MQ=255)
                                                                                                 gCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGtt   <  2:747299/139‑1 (MQ=255)
                                                                                                  cGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGttt  >  1:582141/1‑139 (MQ=255)
                                                                                                  cGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGttt  >  1:418490/1‑139 (MQ=255)
                                                                                                   ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATa                                                                                                  <  2:872100/42‑1 (MQ=255)
                                                                                                   ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATa                                                                                                  >  1:872100/1‑42 (MQ=255)
                                                                                                   ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTcagccag                                      <  2:788941/102‑1 (MQ=255)
                                                                                                   ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTcagccag                                      >  1:788941/1‑102 (MQ=255)
                                                                                                                  |                                                                                                                          
AAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTT  >  NC_000913/2130739‑2130975

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: