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.2% 348.8 / 15.1 83coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base . (48/31);  total (50/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

Only 100 of 103 total aligned reads displayed.

GCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801023‑3801117
                                                 |                                             
gCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               >  1:283664/1‑50 (MQ=255)
gCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               >  1:2572180/1‑50 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              >  1:1578903/1‑50 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               <  2:5167106/49‑1 (MQ=255)
  tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  2:4181184/1‑50 (MQ=255)
  tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  1:1406640/1‑50 (MQ=255)
   ataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                            >  2:2387417/1‑50 (MQ=255)
   ataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                            >  1:1256457/1‑50 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAGAAGATGTTTATATCttttt                                           >  2:2144212/1‑50 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:1949609/1‑50 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  2:401179/1‑50 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                            <  1:1561160/49‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  2:4691099/48‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  2:3517916/48‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:815368/48‑1 (MQ=255)
     aaTTCCTGAGATATAATGATGTTGCACTAAACGATGTTTATATCttttt                                           >  1:3032908/1‑49 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  2:253185/50‑2 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                        <  1:3797273/50‑3 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         >  1:5139426/1‑48 (MQ=255)
        tCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                       >  1:5655537/1‑50 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      >  1:44741/1‑49 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                       <  2:2160329/49‑1 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:3757751/49‑2 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:4542117/1‑50 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  2:4379018/1‑50 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  2:3045623/50‑1 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  2:2177969/50‑1 (MQ=255)
           tGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:3108068/1‑49 (MQ=39)
           tGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:5440593/1‑49 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:1278598/48‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  2:650913/1‑50 (MQ=255)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:4545029/1‑50 (MQ=255)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 <  1:1802870/50‑1 (MQ=255)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:3248933/1‑50 (MQ=255)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  2:5725623/49‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  2:2504158/49‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  1:5015423/49‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  1:4140990/49‑1 (MQ=39)
              gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  2:4183596/50‑1 (MQ=255)
              gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:4541192/1‑49 (MQ=39)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              <  2:3168904/50‑1 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              <  1:2871594/50‑1 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             <  1:1790999/49‑1 (MQ=39)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            <  1:4583236/50‑1 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            <  2:4551954/50‑1 (MQ=255)
                   aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:5341212/1‑49 (MQ=39)
                   aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:4123114/1‑49 (MQ=39)
                   aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  2:4819760/1‑49 (MQ=39)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  2:898085/49‑1 (MQ=39)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          <  2:3810549/50‑1 (MQ=255)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          >  1:4999165/1‑50 (MQ=255)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         <  1:644225/50‑1 (MQ=255)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         <  1:1564659/50‑1 (MQ=255)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         >  2:1154378/1‑50 (MQ=255)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         <  2:4002193/50‑1 (MQ=255)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         >  2:1620627/1‑50 (MQ=255)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         >  2:3804254/1‑50 (MQ=255)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         <  1:4781701/50‑1 (MQ=255)
                      gatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        >  1:97482/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       <  2:2848325/50‑1 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       <  1:3664117/50‑1 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       <  2:5660806/50‑1 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       >  2:1694937/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        <  1:2915252/49‑1 (MQ=39)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:5571691/1‑50 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      <  1:2809483/50‑1 (MQ=255)
                         gTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:3588619/1‑50 (MQ=255)
                         gTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  2:3177341/1‑50 (MQ=255)
                         gTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  2:854268/1‑50 (MQ=255)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGc                    <  2:2019395/50‑1 (MQ=255)
                           tGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   >  1:3022522/1‑50 (MQ=255)
                            gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  >  1:3869477/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 >  2:5610833/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 <  1:5649062/50‑1 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 >  1:3662303/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  2:4080237/49‑1 (MQ=39)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               >  1:3054444/1‑50 (MQ=255)
                                 aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             >  2:1301652/1‑50 (MQ=255)
                                 aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             >  2:4708303/1‑50 (MQ=255)
                                 aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             >  1:180710/1‑50 (MQ=255)
                                   aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTgg           <  2:5150994/50‑1 (MQ=255)
                                   aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  2:5245858/1‑49 (MQ=39)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  2:2017164/1‑48 (MQ=39)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  2:5554520/1‑48 (MQ=39)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:1851059/1‑48 (MQ=39)
                                     gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:3510273/1‑47 (MQ=39)
                                      aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGgaga        <  1:2297910/50‑1 (MQ=255)
                                      aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          >  1:4186313/1‑48 (MQ=39)
                                        gTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:5723886/1‑50 (MQ=255)
                                        gTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:698594/1‑50 (MQ=255)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGTGAGt      >  2:33329/1‑49 (MQ=38)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  2:3594175/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  2:4863498/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:2454097/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  2:3907858/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTc     <  2:2324844/50‑1 (MQ=255)
                                          ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  2:4740805/1‑50 (MQ=255)
                                          ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  2:3453781/1‑50 (MQ=255)
                                          ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:1018702/1‑50 (MQ=255)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:1750087/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: