Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,798,234 A→T 9.8% intergenic (+27/+5) waaL → / ← waaU O‑antigen ligase/putative ADP‑heptose:LPS heptosyltransferase 4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,798,2340AT 9.8% 144.6 / 7.9 51intergenic (+27/+5)waaL/waaUO‑antigen ligase/putative ADP‑heptose:LPS heptosyltransferase 4
Reads supporting (aligned to +/- strand):  ref base A (17/29);  new base T (2/3);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

TGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATA  >  NC_000913/3798095‑3798376
                                                                                                                                           |                                                                                                                                                
tgGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaaa                                                                                                                                               >  1:695235/1‑143 (MQ=255)
tgGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTAaaaa                                                                                                                                               >  2:211844/1‑143 (MQ=255)
gggggAGTGGACTCAGTGATGCGATCATATGGGCACGGAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAAGAAACAAGTTTAAGAAGTGAGTTaaaa                                                                                                                                               <  1:961750/142‑1 (MQ=255)
 ggggATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                                                              >  1:403339/1‑143 (MQ=255)
    gaatGGACTGAGTGATGTGATCATATGGGAACGCAGCAGTCCAATTATCATTATATCCGCTATAGTCCTCATACACGTCATTAACAATCGTAACAAAACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                                                                           <  1:1013879/140‑1 (MQ=255)
    gATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                                                                           <  1:778566/143‑1 (MQ=255)
       tGGTCTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCTCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATACCCAAGTTTAAGAAGTGAGTTAAAACTCActt                                                                                                                                        <  2:640748/143‑1 (MQ=255)
       tGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActt                                                                                                                                        <  1:620926/143‑1 (MQ=255)
       gggACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActt                                                                                                                                        <  2:328673/142‑1 (MQ=255)
        ggACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                                                                       <  1:539937/143‑1 (MQ=255)
         gACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtct                                                                                                                                      >  2:449385/1‑143 (MQ=255)
           cTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCACTTCTTa                                                                                                                                    <  2:6339/143‑1 (MQ=255)
            tCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTAt                                                                                                                                   <  1:606973/143‑1 (MQ=255)
                tgatgtgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtat                                                                                                                               >  1:189090/1‑143 (MQ=255)
                    gtgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTACGAATAAACAAGTTTAAGAAGTGAGTGAAAACTCACTTCTTATCtataaga                                                                                                                           >  2:861718/1‑140 (MQ=255)
                       atCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATCCAACtt                                                                                                                        >  1:488565/1‑143 (MQ=255)
                       atCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACtt                                                                                                                        >  1:270288/1‑143 (MQ=255)
                          atatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAt                                                                                                                     <  1:695832/143‑1 (MQ=255)
                              gggCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATctct                                                                                                                 <  1:371068/143‑1 (MQ=255)
                              gggCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATctct                                                                                                                 <  1:897605/143‑1 (MQ=255)
                                   cGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACGATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACt                                                                                                            <  1:147753/143‑1 (MQ=255)
                                      agcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGaa                                                                                                         >  1:699307/1‑143 (MQ=255)
                                      agcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGaa                                                                                                         <  2:602141/143‑1 (MQ=255)
                                             cAATTATCATTATATCCGCGATAGTCCGCTTACTCGTCATTAATAATCGTAACAATACAAGTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTa                                                                                                  <  1:64801/143‑1 (MQ=255)
                                                ttatacttATATCCGCTATAGTACTCTTACTCGTCATTAATGATCGTAACAATACAATTAATTAAGAATAAACACGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAg                                                                                               <  2:241915/137‑1 (MQ=255)
                                                    cATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                           <  1:540391/143‑1 (MQ=255)
                                                           tCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                                                                    <  2:687307/143‑1 (MQ=255)
                                                               cTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGatat                                                                                <  1:672021/143‑1 (MQ=255)
                                                                   aGTCCTCTTACTCGTTATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCttt                                                                           <  2:392293/143‑1 (MQ=255)
                                                                             cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTatat                                                                 <  1:138420/143‑1 (MQ=255)
                                                                               cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGt                                                               <  1:108879/143‑1 (MQ=255)
                                                                               cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTAAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGt                                                               >  1:1004153/1‑143 (MQ=255)
                                                                                 tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTGACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTGTATGTTg                                                             >  2:990281/1‑143 (MQ=255)
                                                                                      aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt                                                        <  1:668426/143‑1 (MQ=255)
                                                                                             gTAACAAGACAATTAATTAAGAATCAACAATTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTCACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt                                                <  2:310422/143‑1 (MQ=255)
                                                                                              tAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTg                                               >  2:237567/1‑143 (MQ=255)
                                                                                               aaCAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAACCTCACTTCTTATCTATACAACTTAATCTCTTTTCTGAATTAGTTAAAGTTTCGGTTTCATTTCGTTGGGCCCAAAATGTTGGGGAGACTATTTGt                                              <  2:754409/143‑1 (MQ=255)
                                                                                                         ttagttaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTat                                    <  1:668770/143‑1 (MQ=255)
                                                                                                         ttaattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTat                                    <  1:115464/143‑1 (MQ=255)
                                                                                                         ttaattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTat                                    <  1:608749/143‑1 (MQ=255)
                                                                                                               aaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATgatgat                              >  1:639439/1‑143 (MQ=255)
                                                                                                                aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGAtt                             >  1:789255/1‑143 (MQ=255)
                                                                                                                   aTAAACAAGTTTAAGAAGTGAGTTAAAGCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCgg                          <  2:952021/143‑1 (MQ=255)
                                                                                                                       aCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGAc                      >  1:154420/1‑143 (MQ=255)
                                                                                                                          aGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATTGATATCTTTTACCGTATACGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCaa                   >  1:967510/1‑143 (MQ=255)
                                                                                                                           gTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCaaa                  <  1:402888/143‑1 (MQ=255)
                                                                                                                           gTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACCGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGCTGATTCGGGGACCaaa                  >  2:853425/1‑143 (MQ=255)
                                                                                                                                 gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAAt            >  1:806794/1‑143 (MQ=255)
                                                                                                                                 gaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAAt            >  1:102991/1‑143 (MQ=255)
                                                                                                                                 gaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAAt            <  1:716399/143‑1 (MQ=255)
                                                                                                                                  aagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTAGAGGTTGGGGATACTATTTGTATAGATTTATGGTGATTCGGGGACAAAATTAAATg           >  1:341167/1‑143 (MQ=255)
                                                                                                                                  aagTGAGTTAAAACTCACTTCTTATCTATACAACTTGATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATg           >  1:681370/1‑143 (MQ=255)
                                                                                                                                      gAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTGTTACCGTATATGTTGGGGATACTGTTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAg       <  1:289400/143‑1 (MQ=255)
                                                                                                                                        gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAggg     <  1:211844/143‑1 (MQ=255)
                                                                                                                                           aaaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTAAATTACGAGTGTTTATAGATTTATGATGATTCGGGGCCCAAATTAAATGCGAGGGatg  >  1:240242/1‑142 (MQ=255)
                                                                                                                                           |                                                                                                                                                
TGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATA  >  NC_000913/3798095‑3798376

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