Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,115,706 A→G F252S (TTC→TCC)  cpsG ← colanic acid biosynthesis phosphomannomutase CpsG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,115,7060AG100.0% 101.7 / NA 29F252S (TTC→TCC) cpsGcolanic acid biosynthesis phosphomannomutase CpsG
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/16);  total (13/16)

GGAGAGACGTGGATCGTGGATGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGTTAGGAATACCGTTGGGGAAATTGCCGT  >  NZ_CP009273/2115591‑2115832
                                                                                                                   |                                                                                                                              
ggAGAGACGTGGATCGTGGATGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCa                                                                                                         <  2:17694/139‑1 (MQ=255)
      aCGTGGATCGTGGATGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCa                                                                                                   <  1:215187/139‑1 (MQ=255)
                tggatGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCatat                                                                                         >  2:428930/1‑139 (MQ=255)
                tggatGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCatat                                                                                         >  1:67101/1‑139 (MQ=255)
                           cgcgCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCgcgc                                                                                   <  2:32228/134‑1 (MQ=255)
                           cgcgCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCgcgc                                                                                   <  2:226142/134‑1 (MQ=255)
                           cgcgCCGGGATGTTTTTCGAGGAATGCTTCTGCCAACAGGCGGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCgcgc                                                                                   <  2:457310/134‑1 (MQ=255)
                                cGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTtgat                                                                         <  1:559833/139‑1 (MQ=255)
                                 gggATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCgcgc                                                                                   >  1:393901/1‑128 (MQ=255)
                                 gggATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCgcgc                                                                                   <  2:393901/128‑1 (MQ=255)
                                                           ccaccaGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCgcgc                                                                                   >  2:168938/1‑102 (MQ=255)
                                                           ccAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCgcgc                                                                                   <  1:168938/102‑1 (MQ=255)
                                                           ccAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGc                                               <  2:420987/138‑1 (MQ=255)
                                                           ccAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCa                                              >  2:187981/1‑139 (MQ=255)
                                                                  gCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCgg                                        >  2:374303/1‑138 (MQ=255)
                                                                     gACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTg                                                               >  2:150762/1‑112 (MQ=255)
                                                                     gACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTg                                                               <  1:150762/112‑1 (MQ=255)
                                                                       cAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCagtag                                  <  1:428930/139‑1 (MQ=255)
                                                                       cAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCagtag                                  <  1:374303/139‑1 (MQ=255)
                                                                            tagtagCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGt                             <  2:316664/139‑1 (MQ=255)
                                                                              gtagCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGtt                            <  1:187981/138‑1 (MQ=255)
                                                                                    cTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGTTAGGAATa                     >  1:97649/1‑139 (MQ=255)
                                                                                    cTCAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGTTAGGAATa                     >  1:42476/1‑139 (MQ=255)
                                                                                      cAATAAACTGCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGTTAGGAATAcc                   >  2:333936/1‑139 (MQ=255)
                                                                                               gCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGTTAGGAATACCGTTGGGGa           >  1:455822/1‑138 (MQ=255)
                                                                                               gCCCTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGTTAGGAATACCGTTGGGGa           >  2:535426/1‑138 (MQ=255)
                                                                                                 ccTTTTTCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGTTAGGAATACCGTTGGGGAAAt        >  1:549446/1‑139 (MQ=255)
                                                                                                     tttCGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGTTAGGAATACCGTTGGGGAAATTg      <  2:549446/137‑1 (MQ=255)
                                                                                                        cGTCAAACAGGGAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGTTAGGAATACCGTTGGGGAAATTGCCGt  <  1:70988/138‑1 (MQ=255)
                                                                                                                   |                                                                                                                              
GGAGAGACGTGGATCGTGGATGATCTTCGCGCCGGGATTTTTTTCGAGGAATGCTTCTGCCAACAGGCCGACAATGTAGTAGCCCTCAATAAACTGCCCTTTTTCGTCAAACAGGAAACAGCGGTCAAAATCGCCATCAAAAGCAATGCCCATATCCGCGCCGTGTTTGATGACCGCATTGCGGGTGTCGTCGCGGCATTCCGGCAGTAGTGGGTTAGGAATACCGTTGGGGAAATTGCCGT  >  NZ_CP009273/2115591‑2115832

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