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.95.5% 657.1 / 29.9 154coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base . (83/64);  total (84/70)
Fisher's exact test for biased strand distribution p-value = 4.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 192 total aligned reads displayed.

GCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801023‑3801117
                                                 |                                             
gCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               >  1:4138562/1‑50 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               <  1:3153759/49‑1 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               <  1:4690508/49‑1 (MQ=255)
  tataATTCCTGAGATATAATGATGTTGCGCTAAAAGATGTTTATATCttt                                             >  1:1675604/1‑50 (MQ=255)
  tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  1:2411651/1‑50 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           <  1:1016515/50‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:6970460/1‑50 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:6636265/1‑50 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           <  1:1754055/50‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:7002027/48‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:6611687/48‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:5040706/48‑1 (MQ=255)
     aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:2047877/1‑49 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:2261377/50‑2 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          >  1:3367215/1‑49 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          >  1:3340823/1‑49 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                        <  1:2369215/50‑3 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                        <  1:3007563/50‑3 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         >  1:4312970/1‑48 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         >  1:4009711/1‑48 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         >  1:4553582/1‑48 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATGTTTTTTaaat                                      >  1:2067335/1‑49 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                       <  1:6467692/49‑1 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                       <  1:4771760/49‑1 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:1982560/49‑2 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:2391638/49‑2 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:3676672/49‑2 (MQ=255)
            gagaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   >  1:5900535/1‑49 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:4274576/48‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:4112291/48‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:5836804/1‑50 (MQ=255)
              gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:2026773/50‑1 (MQ=255)
               atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                               <  1:405068/50‑1 (MQ=255)
               atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                >  1:3329063/1‑49 (MQ=39)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              >  1:6090794/1‑50 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              >  1:911387/1‑50 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              >  1:6757757/1‑50 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              >  1:1461252/1‑50 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              <  1:2978754/50‑1 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              <  1:256153/50‑1 (MQ=255)
                 ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:6365289/1‑50 (MQ=255)
                 ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:3821570/1‑50 (MQ=255)
                   aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:4729322/1‑49 (MQ=39)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:5084870/49‑1 (MQ=39)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         <  1:5200669/50‑1 (MQ=255)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         <  1:6794700/50‑1 (MQ=255)
                      gatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        >  1:5126258/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTGTATCTTTTTTAAATTCAATAATTGAttt                       >  1:346067/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       >  1:1506841/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       <  1:6262374/50‑1 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       >  1:6540186/1‑50 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:112966/1‑50 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:3438278/1‑50 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:713788/1‑50 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      <  1:5935867/50‑1 (MQ=255)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:6118810/1‑49 (MQ=39)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     <  1:2731788/49‑1 (MQ=39)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     <  1:7435615/49‑1 (MQ=39)
                           tGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   >  1:1419535/1‑50 (MQ=255)
                           tGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   >  1:1262341/1‑50 (MQ=255)
                            gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   <  1:5411799/49‑1 (MQ=39)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:5512217/49‑1 (MQ=39)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:23320/49‑1 (MQ=39)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:2081990/49‑1 (MQ=39)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                >  1:130995/1‑50 (MQ=255)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                <  1:7350350/50‑1 (MQ=255)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               >  1:4888823/1‑50 (MQ=255)
                                tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              >  1:2286008/1‑50 (MQ=255)
                                 aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             >  1:2490355/1‑50 (MQ=255)
                                 aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             >  1:1584148/1‑50 (MQ=255)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             <  1:7220119/49‑1 (MQ=39)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             <  1:1675536/49‑1 (MQ=39)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             <  1:590622/49‑1 (MQ=39)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:1331911/1‑50 (MQ=255)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:1063032/1‑50 (MQ=255)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:4462281/1‑48 (MQ=39)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTCATTTCGCATCTCGTg            >  1:6879303/1‑48 (MQ=38)
                                     gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:6067468/1‑47 (MQ=39)
                                     gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:543134/1‑47 (MQ=39)
                                     gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:1919566/1‑47 (MQ=39)
                                      aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGgaga        <  1:1613434/50‑1 (MQ=255)
                                      aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          >  1:3253534/1‑48 (MQ=39)
                                      aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          >  1:5121405/1‑48 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:5231166/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:6419508/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:6922385/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTc     <  1:7281419/50‑1 (MQ=255)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTc     <  1:5262295/50‑1 (MQ=255)
                                          ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:5583500/1‑50 (MQ=255)
                                          ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:3722696/1‑50 (MQ=255)
                                          ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:2831041/1‑50 (MQ=255)
                                          ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:4631251/1‑50 (MQ=255)
                                           tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   >  1:3309204/1‑50 (MQ=255)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   <  1:1210461/49‑1 (MQ=39)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   <  1:4423581/49‑1 (MQ=39)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   <  1:4808053/49‑1 (MQ=39)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   <  1:2881077/49‑1 (MQ=39)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:5325266/1‑50 (MQ=255)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:7325728/1‑50 (MQ=255)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:5682670/1‑50 (MQ=255)
                                                 |                                             
GCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801023‑3801117

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