Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,118,164 G→A R81C (CGT→TGT)  wcaK ← colanic acid biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,118,1640GA100.0% 90.6 / NA 29R81C (CGT→TGT) wcaKcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/14);  total (15/14)

GATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCA  >  NC_000913/2118039‑2118289
                                                                                                                             |                                                                                                                             
gATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTa                                                                                                                  <  1:927516/139‑1 (MQ=255)
gATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTa                                                                                                                  <  1:431414/139‑1 (MQ=255)
          tcgtAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGa                                                                                                        >  1:1034139/1‑139 (MQ=255)
          tcgtAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGa                                                                                                        >  1:968963/1‑139 (MQ=255)
          tcgtAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGa                                                                                                        >  1:531299/1‑139 (MQ=255)
                ccTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAAcgc                                                                                                  <  2:968963/139‑1 (MQ=255)
                      agcagGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCccgccg                                                                                            <  2:194125/139‑1 (MQ=255)
                        cagGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCcgcg                                                                                          >  2:560479/1‑139 (MQ=255)
                             gcACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtgtt                                                                                     <  1:560479/139‑1 (MQ=255)
                                   aGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCccgccg                                                                                            <  1:255515/126‑1 (MQ=255)
                                   aGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCccgccg                                                                                            >  2:255515/1‑126 (MQ=255)
                                                       atggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGAc                                                                                                       <  1:736063/95‑1 (MQ=255)
                                                       atggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGAc                                                                                                       >  2:736063/1‑95 (MQ=255)
                                                              tatTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGcc                                                    <  1:444636/139‑1 (MQ=255)
                                                                                         cgcgTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGa                                                         >  1:742312/1‑107 (MQ=255)
                                                                                         cgcgTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGa                                                         <  2:742312/107‑1 (MQ=255)
                                                                                             tGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGgaagaa                     >  1:271348/1‑139 (MQ=255)
                                                                                                           gATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCccgccg                                                                                            <  1:226268/54‑1 (MQ=255)
                                                                                                           gATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCccgccg                                                                                            >  2:226268/1‑54 (MQ=255)
                                                                                                             tGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGc     >  2:783249/1‑139 (MQ=255)
                                                                                                                tGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGc                                                                     <  1:711723/72‑1 (MQ=255)
                                                                                                                tGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGc                                                                     >  2:711723/1‑72 (MQ=255)
                                                                                                                tGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCg       <  1:188454/134‑1 (MQ=255)
                                                                                                                tGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCg       >  2:188454/1‑134 (MQ=255)
                                                                                                                tGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtca  >  2:1009655/1‑139 (MQ=255)
                                                                                                                  gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGgaag                       <  2:513503/116‑1 (MQ=255)
                                                                                                                  gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGgaag                       >  1:513503/1‑116 (MQ=255)
                                                                                                                      cggcgACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAg                                                             >  2:375011/1‑74 (MQ=255)
                                                                                                                      cggcgACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAg                                                             <  1:375011/74‑1 (MQ=255)
                                                                                                                             |                                                                                                                             
GATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCA  >  NC_000913/2118039‑2118289

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