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.9% 209.2 / 4.7 49coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (22/25);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAT  >  NC_000913/3801024‑3801116
                                                |                                            
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                             >  1:5103663/1‑50 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                             <  1:4517088/50‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                             <  1:4159761/50‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                              <  1:1069990/49‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                          <  1:3962568/50‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                            <  1:1403657/48‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGAAGTTTATATCttttt                                          <  1:1090763/50‑1 (MQ=255)
    aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                          >  1:2090225/1‑49 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                        <  1:3029461/50‑2 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                        <  1:4656212/50‑2 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                     >  1:1986492/1‑49 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                     <  1:1907923/49‑2 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                   >  1:3838415/1‑50 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                   >  1:1053515/1‑50 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                   >  1:1976461/1‑50 (MQ=255)
          tGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                   >  1:4873028/1‑49 (MQ=39)
           gagaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                  >  1:3833725/1‑49 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                  <  1:62675/48‑1 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                <  1:3189186/50‑1 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                 <  1:3454362/49‑1 (MQ=39)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                               <  1:3783948/50‑1 (MQ=255)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                              <  1:2864643/50‑1 (MQ=255)
               tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                             <  1:1235371/50‑1 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                            >  1:1999666/1‑50 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                            >  1:3602358/1‑50 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                            >  1:3859130/1‑50 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                            <  1:1018165/49‑1 (MQ=39)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                            <  1:3907360/49‑1 (MQ=39)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                           <  1:1404578/50‑1 (MQ=255)
                  aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                           >  1:1616406/1‑49 (MQ=39)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                         <  1:1521986/50‑1 (MQ=255)
                    tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                        <  1:2508795/50‑1 (MQ=255)
                     gatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                       >  1:1911007/1‑50 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                      >  1:690897/1‑50 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                      <  1:4791535/50‑1 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                      >  1:2107030/1‑50 (MQ=255)
                       tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                     <  1:3933298/50‑1 (MQ=255)
                       tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                     <  1:792277/50‑1 (MQ=255)
                        gTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                    >  1:4274258/1‑50 (MQ=255)
                         ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                    <  1:1374282/49‑1 (MQ=39)
                           gCACTAAGAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                  <  1:1383681/49‑1 (MQ=38)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                  <  1:432786/49‑1 (MQ=39)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                  <  1:2034214/49‑1 (MQ=39)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                 >  1:2096372/1‑50 (MQ=255)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct               >  1:1162886/1‑50 (MQ=255)
                              cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc              >  1:3336989/1‑50 (MQ=255)
                               tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg             >  1:2972330/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt            <  1:397223/49‑1 (MQ=39)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt            <  1:4670658/49‑1 (MQ=39)
                                  aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTgg          <  1:4485774/50‑1 (MQ=255)
                                  aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTgg          <  1:3089496/50‑1 (MQ=255)
                                  aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTgg          <  1:117320/50‑1 (MQ=255)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGgaga       <  1:1131692/50‑1 (MQ=255)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga         >  1:4912690/1‑48 (MQ=39)
                                      tttttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag      >  1:3238856/3‑50 (MQ=255)
                                      tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag      >  1:4638892/1‑50 (MQ=255)
                                       gTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt     >  1:3931057/1‑50 (MQ=255)
                                        tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt     >  1:4016165/1‑49 (MQ=39)
                                           atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt  <  1:1871852/49‑1 (MQ=39)
                                                |                                            
CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAT  >  NC_000913/3801024‑3801116

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: