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.81.8% 299.0 / 109.8 110coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (8/12);  new base . (67/23);  total (75/35)
Fisher's exact test for biased strand distribution p-value = 6.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 135 total aligned reads displayed.

CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTC  >  NC_000913/3801024‑3801121
                                                |                                                 
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                                  <  1:1351324/50‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                                  >  1:2859320/1‑50 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                                  <  1:1829829/50‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                                  >  2:4836387/1‑50 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                                   <  1:832978/49‑1 (MQ=255)
  ataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                                >  2:1076023/1‑50 (MQ=255)
  ataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                                >  1:718609/1‑50 (MQ=255)
  ataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                                >  1:2271733/1‑50 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                               <  2:2098178/50‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                               <  2:1788073/50‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                                <  2:3637296/49‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                                <  2:5108949/49‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                                 <  2:1627627/48‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                                 <  2:4796061/48‑1 (MQ=255)
    aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                               >  2:1566983/1‑49 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                             <  1:5163932/50‑2 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                             <  2:1656121/50‑2 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                              <  1:651244/49‑1 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                              <  1:5083587/49‑1 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                              <  2:2446744/49‑1 (MQ=255)
      ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                            <  2:2039484/50‑3 (MQ=255)
      ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                            <  1:5324434/50‑3 (MQ=255)
      ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTa                                            <  1:3460086/50‑1 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                          >  2:5267313/1‑49 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                          >  2:2614773/1‑49 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                           <  1:2402162/49‑1 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                          <  2:5231759/49‑2 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                          <  2:217372/49‑2 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCa                                      <  1:908596/50‑1 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                       <  2:2474754/48‑1 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                       <  2:5023285/48‑1 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                      <  2:2621157/49‑1 (MQ=39)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                     >  2:618423/1‑49 (MQ=39)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                     >  1:4354727/1‑49 (MQ=39)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                    >  1:2478465/1‑49 (MQ=39)
               tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                                  >  2:3764246/1‑50 (MQ=255)
               tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                                  <  2:113437/50‑1 (MQ=255)
               tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                                  <  2:4374885/50‑1 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCaataa                                  >  2:390482/1‑50 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                                 >  2:4177918/1‑50 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                                 >  1:3557055/1‑50 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                                 <  1:4443229/49‑1 (MQ=39)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                                >  1:968020/1‑50 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                                <  2:2277144/50‑1 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                                >  1:444211/1‑50 (MQ=255)
                  aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                                >  1:4840672/1‑49 (MQ=39)
                  aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                                >  1:4688667/1‑49 (MQ=39)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                              <  2:2729200/50‑1 (MQ=255)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                              >  2:1479468/1‑50 (MQ=255)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                              <  1:3402201/50‑1 (MQ=255)
                    tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                             >  1:1972488/1‑50 (MQ=255)
                    tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                             <  2:45233/50‑1 (MQ=255)
                     gatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                            >  2:2572498/1‑50 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                           >  1:1970118/1‑50 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                           <  2:2313855/50‑1 (MQ=255)
                        gTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                         >  2:3887731/1‑50 (MQ=255)
                         ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                         >  2:2019621/1‑49 (MQ=39)
                         ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                         >  1:2961505/1‑49 (MQ=39)
                          tGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCg                         >  2:3802818/1‑49 (MQ=255)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                       <  2:4717834/49‑1 (MQ=39)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                      >  2:2178762/1‑50 (MQ=255)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                      >  1:5168148/1‑50 (MQ=255)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                      >  1:587480/1‑50 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                     >  2:145278/1‑50 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                     >  1:1891908/1‑50 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                     >  1:5230702/1‑50 (MQ=255)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                    <  2:2335213/50‑1 (MQ=255)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                    >  2:2066309/1‑50 (MQ=255)
                              cTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCAtct                    <  1:2788987/50‑1 (MQ=255)
                              cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc                   >  1:1977807/1‑50 (MQ=255)
                              cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc                   >  2:4172571/1‑50 (MQ=255)
                               tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg                  >  1:4244868/1‑50 (MQ=255)
                                aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt                 >  2:3912387/1‑50 (MQ=255)
                                aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCACCTCGt                 >  2:1056611/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCg                  <  1:1186592/49‑1 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg                >  1:4863965/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg                >  2:1231997/1‑50 (MQ=255)
                                  aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg                >  1:1283546/1‑49 (MQ=39)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg                >  1:4659892/1‑48 (MQ=39)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg                >  1:4494258/1‑48 (MQ=39)
                                    gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg                >  2:3461587/1‑47 (MQ=39)
                                    gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg                >  1:1759139/1‑47 (MQ=39)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga              >  2:2861115/1‑48 (MQ=39)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga              >  2:2764360/1‑48 (MQ=39)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga              >  2:3684720/1‑48 (MQ=39)
                                      tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag           >  1:4099304/1‑50 (MQ=255)
                                      tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag           <  2:5219581/50‑1 (MQ=255)
                                      tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag           >  2:5064866/1‑50 (MQ=255)
                                       gTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt          >  1:2445641/1‑50 (MQ=255)
                                       gTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt          >  1:4914916/1‑50 (MQ=255)
                                        tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt          >  1:3747297/1‑49 (MQ=39)
                                         ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa        >  1:3580773/1‑50 (MQ=255)
                                         ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa        >  1:4703525/1‑50 (MQ=255)
                                          tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt       >  1:3066325/1‑50 (MQ=255)
                                          tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt       >  2:4842325/1‑50 (MQ=255)
                                           atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt       <  2:2461289/49‑1 (MQ=39)
                                           atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt       <  2:1769517/49‑1 (MQ=39)
                                           atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc      >  2:4231732/1‑50 (MQ=255)
                                              tCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCt     <  1:2296680/49‑1 (MQ=255)
                                                tttttttAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTc  >  2:4499007/1‑50 (MQ=255)
                                                |                                                 
CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATCTTTC  >  NC_000913/3801024‑3801121

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