Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,801,078 (T)7→6 coding (922/1017 nt) waaJ ← UDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,801,0720T.92.9% 282.2 / 19.9 70coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base . (37/28);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

GCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801023‑3801117
                                                 |                                             
gCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               >  1:5172809/1‑50 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              <  1:626548/50‑1 (MQ=255)
  tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  1:3635384/1‑50 (MQ=255)
  tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  1:7201062/1‑50 (MQ=255)
  tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  1:1341949/1‑50 (MQ=255)
  tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  1:6142156/1‑50 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:1339413/1‑50 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           <  1:1950469/50‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                            <  1:3079071/49‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:523270/48‑1 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:7237986/50‑2 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:3711858/50‑2 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          <  1:5542937/49‑1 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                        <  1:4951164/50‑3 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                        <  1:4858623/50‑3 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                        <  1:4534438/50‑3 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         >  1:1573978/1‑48 (MQ=255)
        tCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                       >  1:525962/1‑50 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      >  1:3440981/1‑49 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:5886212/49‑2 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:929523/49‑2 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:7136510/1‑50 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:2091700/1‑50 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:960247/50‑1 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:1261789/1‑50 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:3043862/50‑1 (MQ=255)
            gagaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   >  1:773307/1‑49 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:3072056/48‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:6104906/1‑50 (MQ=255)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:1782320/1‑50 (MQ=255)
              gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:4939638/50‑1 (MQ=255)
              gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:2955402/50‑1 (MQ=255)
              gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:5996549/50‑1 (MQ=255)
              gaTATAATGATGTTGCACAAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:2593138/1‑49 (MQ=38)
               atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                               <  1:1398477/50‑1 (MQ=255)
               atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                >  1:4086977/1‑49 (MQ=39)
                 ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:2238222/1‑50 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             <  1:678452/49‑1 (MQ=39)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             <  1:956157/49‑1 (MQ=39)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:1117416/1‑50 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:1204105/1‑50 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:7324648/1‑50 (MQ=255)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:7440845/49‑1 (MQ=39)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         <  1:2678561/50‑1 (MQ=255)
                      gatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        >  1:1979712/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       >  1:5328298/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       <  1:2992544/50‑1 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      <  1:6138374/50‑1 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:412524/1‑50 (MQ=255)
                         gTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:2057447/1‑50 (MQ=255)
                         gTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:1425426/1‑50 (MQ=255)
                         gTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:2667144/1‑50 (MQ=255)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     <  1:2847562/49‑1 (MQ=39)
                            gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  >  1:5581135/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 >  1:1741098/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 >  1:479963/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:1609459/49‑1 (MQ=39)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:4313148/49‑1 (MQ=39)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:2376817/49‑1 (MQ=39)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:3521906/49‑1 (MQ=39)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                >  1:262821/1‑50 (MQ=255)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               >  1:3575771/1‑50 (MQ=255)
                                tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              >  1:965625/1‑50 (MQ=255)
                                tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              >  1:3593663/1‑50 (MQ=255)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              <  1:1653845/48‑1 (MQ=39)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             <  1:6887185/49‑1 (MQ=39)
                                   aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTgg           <  1:3061826/50‑1 (MQ=255)
                                   aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTgg           <  1:6288784/50‑1 (MQ=255)
                                   aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:6110630/1‑49 (MQ=39)
                                     gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGgag         <  1:939609/50‑1 (MQ=255)
                                      aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          >  1:6352661/1‑48 (MQ=39)
                                      aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          >  1:3254528/1‑48 (MQ=39)
                                       tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       <  1:6371021/50‑1 (MQ=255)
                                       tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       <  1:6774848/50‑1 (MQ=255)
                                       tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       <  1:4031414/50‑1 (MQ=255)
                                       tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       >  1:1136944/1‑50 (MQ=255)
                                        gTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:2235010/1‑50 (MQ=255)
                                        gTTTATATCTTTTTAAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:5113091/1‑50 (MQ=255)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:3072954/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:2563238/1‑49 (MQ=39)
                                          ttATATCTTTTTTAAATTCAATAATTGATTTCGCCTCTCGTGGAGAGTCa    >  1:937624/1‑50 (MQ=255)
                                           tataTCTGTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGTGTCa    <  1:7026153/49‑1 (MQ=255)
                                           tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    <  1:5574083/49‑1 (MQ=39)
                                           tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    <  1:4315263/49‑1 (MQ=39)
                                           tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   >  1:860484/1‑50 (MQ=255)
                                           tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   >  1:3469996/1‑50 (MQ=255)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:1516830/1‑50 (MQ=255)
                                                 |                                             
GCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801023‑3801117

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