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.96.1% 218.7 / 4.5 51coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (24/25);  total (25/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.98e-01

GCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801023‑3801117
                                                 |                                             
gCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               >  1:203171/1‑50 (MQ=255)
gCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               >  1:1488982/1‑50 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              >  1:2337803/1‑50 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              >  1:441480/1‑50 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              >  1:4188335/1‑50 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               <  1:3278416/49‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                            <  1:1864577/49‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:899037/48‑1 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:4066163/50‑2 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          >  1:3308861/1‑49 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          >  1:1616181/1‑49 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          >  1:1450228/1‑49 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                        <  1:192832/50‑3 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                        <  1:3592805/50‑3 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                        <  1:3781442/50‑3 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         >  1:1412553/1‑48 (MQ=255)
       ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         >  1:2845985/1‑48 (MQ=255)
        tCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                       >  1:262286/1‑50 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                       <  1:2461514/49‑1 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:3626491/50‑1 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:2533261/50‑1 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:2769439/1‑50 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:1159817/1‑50 (MQ=255)
           tGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:1324815/1‑49 (MQ=39)
           tGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:2011041/1‑49 (MQ=39)
           tGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:2871027/1‑49 (MQ=39)
            gagaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   >  1:3016302/1‑49 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  1:2860357/49‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  1:260617/49‑1 (MQ=39)
              gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:482529/50‑1 (MQ=255)
              gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:2497656/50‑1 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              <  1:2316845/50‑1 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              <  1:2281265/50‑1 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              <  1:2758546/50‑1 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             <  1:3094950/49‑1 (MQ=39)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:2447320/1‑50 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            <  1:211654/50‑1 (MQ=255)
                   aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:1175972/1‑49 (MQ=39)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:1087363/49‑1 (MQ=39)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:168/49‑1 (MQ=39)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:1172968/49‑1 (MQ=39)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         >  1:2530505/1‑50 (MQ=255)
                      gatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        >  1:2995236/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       >  1:3221108/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       <  1:2102090/50‑1 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        <  1:1362600/49‑1 (MQ=39)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:2458310/1‑50 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      <  1:3899466/50‑1 (MQ=255)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     <  1:4245244/49‑1 (MQ=39)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     <  1:3061260/49‑1 (MQ=39)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGc                    <  1:2578733/50‑1 (MQ=255)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGc                    <  1:3738160/50‑1 (MQ=255)
                            gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  >  1:3723085/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 >  1:781379/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 <  1:1283825/50‑1 (MQ=255)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                >  1:1172459/1‑50 (MQ=255)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               <  1:4140776/50‑1 (MQ=255)
                                tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              >  1:1651604/1‑50 (MQ=255)
                                tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              >  1:127512/1‑50 (MQ=255)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:1912593/1‑50 (MQ=255)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:235156/1‑48 (MQ=39)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          <  1:3115718/50‑1 (MQ=255)
                                       tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       >  1:3954738/1‑50 (MQ=255)
                                          ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:3479225/1‑50 (MQ=255)
                                           tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    <  1:1047311/49‑1 (MQ=39)
                                           tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   >  1:1089099/1‑50 (MQ=255)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:2104791/1‑50 (MQ=255)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:2274627/1‑50 (MQ=255)
                                                 |                                             
GCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801023‑3801117

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