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.90.9% 325.6 / 39.0 88coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base . (52/28);  total (56/32)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 114 total aligned reads displayed.

GCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAT  >  NC_000913/3801023‑3801116
                                                 |                                            
gCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                              >  1:3153615/1‑50 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                             >  1:3872833/1‑50 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                             <  1:5090770/50‑1 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                              <  1:6934376/49‑1 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                              <  1:2602305/49‑1 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                              <  1:2255249/49‑1 (MQ=255)
   ataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                           >  1:5151679/1‑50 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                          <  1:7001824/50‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                           <  1:3675270/49‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                           <  1:4325374/49‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                           <  1:448315/49‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                            <  1:8213120/48‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                            <  1:4368956/48‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                            <  1:1656499/48‑1 (MQ=255)
     aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                          >  1:6306584/1‑49 (MQ=255)
     aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                          >  1:1615298/1‑49 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                        <  1:8615321/50‑2 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                        <  1:3250597/50‑2 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                        <  1:5449534/50‑2 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                        <  1:5691139/50‑2 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                         <  1:1908240/49‑1 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                        >  1:1527664/1‑48 (MQ=255)
        tCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                      >  1:5858823/1‑50 (MQ=255)
        tCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                      >  1:3540057/1‑50 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                     >  1:4477478/1‑49 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                      <  1:5579054/49‑1 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                      <  1:1038449/49‑1 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                      <  1:7905938/49‑1 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                     <  1:2680629/49‑2 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                   >  1:8071157/1‑50 (MQ=255)
           tGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                   >  1:6972985/1‑49 (MQ=39)
            gagaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                  >  1:507858/1‑49 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                  <  1:1270043/48‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                >  1:7589741/1‑50 (MQ=255)
              gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                >  1:427904/1‑49 (MQ=39)
               atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                              <  1:1098623/50‑1 (MQ=255)
               atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                              <  1:8041960/50‑1 (MQ=255)
               atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                               >  1:7256630/1‑49 (MQ=39)
                 ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                            >  1:2817282/1‑50 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                           <  1:2433613/50‑1 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                           >  1:1270648/1‑50 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                           >  1:3498383/1‑50 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                           >  1:6905722/1‑50 (MQ=255)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                         >  1:8363482/1‑50 (MQ=255)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                         >  1:465984/1‑50 (MQ=255)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                         <  1:1738793/50‑1 (MQ=255)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                        >  1:7918237/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                      >  1:3967909/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                       <  1:7144235/49‑1 (MQ=39)
                       aagTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                       <  1:6529072/47‑1 (MQ=38)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                     >  1:5506885/1‑50 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                     <  1:5500103/50‑1 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                     <  1:5994322/50‑1 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                     >  1:8182489/1‑50 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                     <  1:4380687/50‑1 (MQ=255)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                    <  1:5189428/49‑1 (MQ=39)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                    <  1:6340515/49‑1 (MQ=39)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                    >  1:3156771/1‑49 (MQ=39)
                           tGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                  >  1:6810580/1‑50 (MQ=255)
                            gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                  <  1:3849572/49‑1 (MQ=39)
                            gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                 >  1:8431857/1‑50 (MQ=255)
                            gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                 >  1:4666033/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                >  1:8429510/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                <  1:8429031/50‑1 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                >  1:5044692/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                <  1:3345618/50‑1 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                >  1:5451420/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                >  1:2622185/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                >  1:1032218/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                 <  1:8334551/49‑1 (MQ=39)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct               >  1:243820/1‑50 (MQ=255)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct               >  1:199194/1‑50 (MQ=255)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc              >  1:4702090/1‑50 (MQ=255)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc              >  1:1961251/1‑50 (MQ=255)
                                tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg             >  1:6406327/1‑50 (MQ=255)
                                tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg             >  1:8124911/1‑50 (MQ=255)
                                 aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt            >  1:5073415/1‑50 (MQ=255)
                                 aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt            >  1:606502/1‑50 (MQ=255)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt            <  1:4758368/49‑1 (MQ=39)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt            <  1:36904/49‑1 (MQ=39)
                                   aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg           >  1:5650965/1‑49 (MQ=39)
                                   aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg           >  1:7854218/1‑49 (MQ=39)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg           >  1:7186392/1‑48 (MQ=39)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg           >  1:1668942/1‑48 (MQ=39)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg           >  1:436079/1‑48 (MQ=39)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg           >  1:8447131/1‑48 (MQ=39)
                                     gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg           >  1:4950806/1‑47 (MQ=39)
                                      aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga         >  1:8623323/1‑48 (MQ=39)
                                      aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga         >  1:2493035/1‑48 (MQ=39)
                                        gTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt     >  1:7645946/1‑50 (MQ=255)
                                        gTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt     >  1:6434539/1‑50 (MQ=255)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt     >  1:2345831/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt     <  1:1871211/49‑1 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt     >  1:8537864/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt     >  1:1566757/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTc    <  1:2207286/50‑1 (MQ=255)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTc    <  1:7061779/50‑1 (MQ=255)
                                           tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt  >  1:2313167/1‑50 (MQ=255)
                                           tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt  >  1:651244/1‑50 (MQ=255)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt  <  1:4211367/49‑1 (MQ=39)
                                                 |                                            
GCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAT  >  NC_000913/3801023‑3801116

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