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.91.9% 620.2 / 64.2 160coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (4/9);  new base . (97/50);  total (101/59)
Fisher's exact test for biased strand distribution p-value = 1.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 188 total aligned reads displayed.

CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801024‑3801117
                                                |                                             
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               <  1:3059930/49‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               <  1:2831711/49‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           <  1:4479917/50‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           <  1:4849247/50‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           <  1:1777228/50‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:4319466/48‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:3851168/48‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:2167288/48‑1 (MQ=255)
    aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:941152/1‑49 (MQ=255)
    aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:1552271/1‑49 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:1055093/50‑2 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:312827/50‑2 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          >  1:3465057/1‑49 (MQ=255)
      ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         >  1:610549/1‑48 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      >  1:4329030/1‑49 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      >  1:5064231/1‑49 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                       <  1:5191702/49‑1 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                       <  1:2907553/49‑1 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:4984992/49‑2 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:497093/49‑2 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:1500410/49‑2 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:3388789/50‑1 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:4648351/50‑1 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:5837976/1‑50 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:5408257/1‑50 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:257280/1‑50 (MQ=255)
          tGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:2802716/1‑49 (MQ=39)
           gagaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   >  1:5637427/1‑49 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:3890712/48‑1 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:5844710/48‑1 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:2378905/1‑50 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:2129035/1‑50 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  1:395942/49‑1 (MQ=39)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:3696797/1‑49 (MQ=39)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:3013132/1‑49 (MQ=39)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:1613360/1‑49 (MQ=39)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                               <  1:3606871/50‑1 (MQ=255)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                               <  1:438669/50‑1 (MQ=255)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                               <  1:1702665/50‑1 (MQ=255)
               tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              >  1:6018304/1‑50 (MQ=255)
               tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              >  1:4832570/1‑50 (MQ=255)
               tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              <  1:5979997/50‑1 (MQ=255)
               tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              >  1:2069124/1‑50 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:805922/1‑50 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:3022983/1‑50 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:5393760/1‑50 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:2240026/1‑50 (MQ=255)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:5497010/49‑1 (MQ=39)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          >  1:1028514/1‑50 (MQ=255)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          >  1:38265/1‑50 (MQ=255)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          >  1:5967369/1‑50 (MQ=255)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          >  1:233851/1‑50 (MQ=255)
                    tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         >  1:3142306/1‑50 (MQ=255)
                    tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         <  1:924824/50‑1 (MQ=255)
                    tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         >  1:3275310/1‑50 (MQ=255)
                     gatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        >  1:899407/1‑50 (MQ=255)
                     gatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        >  1:1929432/1‑50 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       >  1:4286200/1‑50 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        <  1:950906/49‑1 (MQ=39)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        <  1:2369563/49‑1 (MQ=39)
                         ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:1401215/1‑49 (MQ=39)
                         ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGc                    <  1:2715572/50‑1 (MQ=255)
                          tGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   >  1:4977508/1‑50 (MQ=255)
                          tGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   >  1:5236538/1‑50 (MQ=255)
                          tGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   >  1:3202592/1‑50 (MQ=255)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   <  1:5700510/49‑1 (MQ=39)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  >  1:3741975/1‑50 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 >  1:3928157/1‑50 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:3224010/49‑1 (MQ=39)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:2057902/49‑1 (MQ=39)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                <  1:84744/50‑1 (MQ=255)
                              cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               >  1:2034201/1‑50 (MQ=255)
                               tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              >  1:704486/1‑50 (MQ=255)
                               tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              >  1:5466719/1‑50 (MQ=255)
                                aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             >  1:24788/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:3807090/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:4708890/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:4260646/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:1146917/1‑50 (MQ=255)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          <  1:3308781/50‑1 (MQ=255)
                                    gATTTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:2195777/1‑47 (MQ=38)
                                    gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:1464005/1‑47 (MQ=39)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          >  1:1565620/1‑48 (MQ=39)
                                      tGTTTATATCTTTTTTTAATTCAATAATTGATTTCGCATCTCGTGGagag       >  1:5746209/1‑50 (MQ=255)
                                      tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       >  1:116538/1‑50 (MQ=255)
                                      tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       >  1:5952008/1‑50 (MQ=255)
                                       gTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:4061252/1‑50 (MQ=255)
                                        tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:5335366/1‑49 (MQ=39)
                                        tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:1821202/1‑49 (MQ=39)
                                         ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:5153894/1‑50 (MQ=255)
                                         ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:2028845/1‑50 (MQ=255)
                                         ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:1841802/1‑50 (MQ=255)
                                         ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:5005989/1‑50 (MQ=255)
                                          tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    <  1:1496278/49‑1 (MQ=39)
                                          tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    <  1:3873228/49‑1 (MQ=39)
                                          tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   >  1:4247476/1‑50 (MQ=255)
                                          tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   >  1:724275/1‑50 (MQ=255)
                                           atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   <  1:5935957/49‑1 (MQ=39)
                                           atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:839365/1‑50 (MQ=255)
                                           atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:5819217/1‑50 (MQ=255)
                                                |                                             
CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801024‑3801117

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