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.97.8% 201.2 / ‑1.2 46coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (17/28);  total (17/29)
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.99e-01

CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCA  >  NC_000913/3801024‑3801115
                                                |                                           
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                            <  1:2095453/50‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                            >  1:472452/1‑50 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                             <  1:1096812/49‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                             <  1:1944814/49‑1 (MQ=255)
 tataATTCCTGTGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                           >  1:1675227/1‑50 (MQ=255)
  ataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                          >  1:320229/1‑50 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                         >  1:1304936/1‑50 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                         >  1:2124358/1‑50 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                          <  1:1873677/49‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                           <  1:1337859/48‑1 (MQ=255)
   taATTCCGGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                         <  1:730272/50‑1 (MQ=255)
    aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                         >  1:323169/1‑49 (MQ=255)
     actCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                        <  1:1857008/47‑1 (MQ=255)
      ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                       >  1:1262085/1‑48 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaa                                     <  1:1701908/49‑1 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                 <  1:2284673/48‑1 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                               <  1:1076450/50‑1 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                <  1:1810615/49‑1 (MQ=39)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                              <  1:897991/50‑1 (MQ=255)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                              <  1:606542/50‑1 (MQ=255)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                              <  1:77051/50‑1 (MQ=255)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                             <  1:2193248/50‑1 (MQ=255)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                             <  1:1441678/50‑1 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                           >  1:743265/1‑50 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                           <  1:217645/49‑1 (MQ=39)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                           <  1:1551866/49‑1 (MQ=39)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                          <  1:1251349/50‑1 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                          >  1:695072/1‑50 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                          <  1:2718450/50‑1 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                          <  1:2425996/50‑1 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                          <  1:1131356/50‑1 (MQ=255)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                         <  1:2693113/49‑1 (MQ=39)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                        >  1:1820045/1‑50 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                     <  1:2776041/50‑1 (MQ=255)
                       tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                    <  1:1661495/50‑1 (MQ=255)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                >  1:863500/1‑50 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc               >  1:1632753/1‑50 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                <  1:989151/49‑1 (MQ=39)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                <  1:1275281/49‑1 (MQ=39)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct              <  1:1447335/50‑1 (MQ=255)
                              cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc             <  1:1519303/50‑1 (MQ=255)
                               tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg            >  1:2670424/1‑50 (MQ=255)
                               tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg            >  1:1903453/1‑50 (MQ=255)
                                aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt           >  1:1738565/1‑50 (MQ=255)
                                aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt           >  1:385774/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg            <  1:159625/48‑1 (MQ=39)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg          >  1:2580637/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg          <  1:1863126/50‑1 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg          <  1:1777974/50‑1 (MQ=255)
                                  aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg          >  1:1233337/1‑49 (MQ=39)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg          >  1:2078291/1‑48 (MQ=39)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga        <  1:367490/50‑1 (MQ=255)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga        <  1:1272563/50‑1 (MQ=255)
                                    gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg          >  1:1707492/1‑47 (MQ=39)
                                    gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg          >  1:1737495/1‑47 (MQ=39)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga        >  1:2692126/1‑48 (MQ=39)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga        >  1:239607/1‑48 (MQ=39)
                                       gTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt    >  1:1652929/1‑50 (MQ=255)
                                        tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt    <  1:677799/49‑1 (MQ=39)
                                          tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa  <  1:1521179/49‑1 (MQ=39)
                                                |                                           
CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCA  >  NC_000913/3801024‑3801115

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: