Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,791,182 A→C Y292S (TAC→TCC)  waaL → O‑antigen ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,791,1820AC100.0% 127.5 / NA 37Y292S (TAC→TCC) waaLO‑antigen ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (29/8);  total (29/8)

GAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTT  >  NZ_CP009273/3791096‑3791263
                                                                                      |                                                                                 
gAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGa                                                                                >  1:985121/1‑90 (MQ=255)
  aTCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGaaa                                                                              >  1:372926/1‑90 (MQ=255)
     gCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTg                                                                           >  2:432296/1‑90 (MQ=255)
     gCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTg                                                                           >  1:459457/1‑90 (MQ=255)
                  aTTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAatat                                                                 >  2:876979/1‑87 (MQ=255)
                  aTTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAatat                                                                 <  1:876979/87‑1 (MQ=255)
                          aCTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGt                                                      >  1:230194/1‑90 (MQ=255)
                          aCTTAAACAGTTATACCAATGCTAATAGTGTTACCTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGt                                                      >  2:879780/1‑90 (MQ=255)
                                   gTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCttttt                                             <  2:505446/90‑1 (MQ=255)
                                           aaTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGa                                     >  1:178309/1‑90 (MQ=255)
                                           aaTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGa                                     >  2:139220/1‑90 (MQ=255)
                                            aTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGAt                                    >  1:126239/1‑90 (MQ=255)
                                             tGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATc                                   >  2:467991/1‑90 (MQ=255)
                                             tGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATc                                   >  2:834794/1‑90 (MQ=255)
                                                      tgtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTc                          >  2:364750/1‑90 (MQ=255)
                                                       gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa                         >  2:705625/1‑90 (MQ=255)
                                                       gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa                         >  2:710584/1‑90 (MQ=255)
                                                       gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa                         >  1:859730/1‑90 (MQ=255)
                                                       gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa                         >  1:688984/1‑90 (MQ=255)
                                                       gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa                         >  1:458677/1‑90 (MQ=255)
                                                       gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa                         >  1:379104/1‑90 (MQ=255)
                                                       gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa                         >  1:282791/1‑90 (MQ=255)
                                                        ttACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCAc                        >  2:936935/1‑90 (MQ=255)
                                                        ttACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCAc                        <  2:958259/90‑1 (MQ=255)
                                                         tACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTa                                                                      <  1:401140/43‑1 (MQ=255)
                                                         tACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTa                                                                      >  2:401140/1‑43 (MQ=255)
                                                          aCTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCAcgc                      >  2:740718/1‑90 (MQ=255)
                                                               ccTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTg                  <  2:174299/89‑1 (MQ=255)
                                                               ccTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTg                  >  1:174299/1‑89 (MQ=255)
                                                               ccTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGa                 >  2:471522/1‑90 (MQ=255)
                                                                  aGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAg              <  2:267912/90‑1 (MQ=255)
                                                                    gTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTa            <  2:871306/90‑1 (MQ=255)
                                                                    gTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTa            <  1:145148/90‑1 (MQ=255)
                                                                        aaGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGaa        >  2:127643/1‑90 (MQ=255)
                                                                             tGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGt   >  2:853044/1‑90 (MQ=255)
                                                                              ggCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGtt  >  2:566902/1‑90 (MQ=255)
                                                                              ggCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGtt  >  1:421976/1‑90 (MQ=255)
                                                                                      |                                                                                 
GAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTT  >  NZ_CP009273/3791096‑3791263

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