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.5% 364.0 / 4.9 81coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (42/37);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801024‑3801117
                                                |                                             
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              <  1:2555834/50‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              >  1:3601689/1‑50 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              >  1:1932307/1‑50 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              <  1:4197791/50‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               <  1:5461185/49‑1 (MQ=255)
 tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  1:2902529/1‑50 (MQ=255)
 tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  1:3349311/1‑50 (MQ=255)
 tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  1:2115837/1‑50 (MQ=255)
  ataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                            >  1:5337109/1‑50 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:2738271/1‑50 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:5626481/48‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:4072092/48‑1 (MQ=255)
    aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:4548061/1‑49 (MQ=255)
    aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:365877/1‑49 (MQ=255)
    aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:3768252/1‑49 (MQ=255)
    aaTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           >  1:4955784/1‑49 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:1877851/50‑2 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:330085/50‑2 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      >  1:2804528/1‑49 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:4983275/49‑2 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:3881520/1‑50 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:4663015/50‑1 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:1787430/50‑1 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:5321066/50‑1 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:6470665/50‑1 (MQ=255)
          tGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:941816/1‑49 (MQ=39)
          tGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:6110907/1‑49 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 <  1:5929916/50‑1 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  1:1756491/49‑1 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  1:2061859/49‑1 (MQ=39)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:4370710/50‑1 (MQ=255)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:4219734/50‑1 (MQ=255)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:4501481/1‑49 (MQ=39)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                               <  1:4037115/50‑1 (MQ=255)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                               <  1:4113183/50‑1 (MQ=255)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                >  1:57450/1‑49 (MQ=39)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                >  1:615331/1‑49 (MQ=39)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:5229587/1‑50 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             <  1:1109507/49‑1 (MQ=39)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             <  1:6547805/49‑1 (MQ=39)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:4535760/1‑50 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            <  1:4978929/50‑1 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            <  1:2577710/50‑1 (MQ=255)
                  aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:3849309/1‑49 (MQ=39)
                  aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:3754649/1‑49 (MQ=39)
                  aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:3598707/1‑49 (MQ=39)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:3844388/49‑1 (MQ=39)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:5943042/49‑1 (MQ=39)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          <  1:2681871/50‑1 (MQ=255)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          <  1:4190989/50‑1 (MQ=255)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          >  1:2800471/1‑50 (MQ=255)
                    tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         <  1:5906752/50‑1 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        <  1:4407884/49‑1 (MQ=39)
                       tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:3728814/1‑50 (MQ=255)
                       tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:5315660/1‑50 (MQ=255)
                       tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      <  1:2839352/50‑1 (MQ=255)
                        gTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:6151503/1‑50 (MQ=255)
                         ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     <  1:1166947/49‑1 (MQ=39)
                         ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:2485998/1‑49 (MQ=39)
                         ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     <  1:6531989/49‑1 (MQ=39)
                          tGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   >  1:4420604/1‑50 (MQ=255)
                          tGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   >  1:6178030/1‑50 (MQ=255)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  >  1:2661341/1‑50 (MQ=255)
                            cACTAAAAGGTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 <  1:3204886/50‑1 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 <  1:404449/50‑1 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 >  1:594996/1‑50 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 <  1:2256421/50‑1 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 >  1:3444330/1‑50 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:4630691/49‑1 (MQ=39)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                >  1:5672685/1‑50 (MQ=255)
                               tAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              >  1:1957524/1‑50 (MQ=255)
                                aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATTTCGt             >  1:5598913/1‑50 (MQ=255)
                                aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             >  1:3495415/1‑50 (MQ=255)
                                aaaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             >  1:1423335/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              <  1:6550922/48‑1 (MQ=39)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              <  1:4714742/48‑1 (MQ=39)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              <  1:2070951/48‑1 (MQ=39)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:5697923/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:1085770/1‑50 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            <  1:6309772/50‑1 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            <  1:1460450/50‑1 (MQ=255)
                                  aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:6091611/1‑49 (MQ=39)
                                  aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:5778889/1‑49 (MQ=39)
                                  aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:2253809/1‑49 (MQ=39)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:2417844/1‑48 (MQ=39)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          <  1:5165438/50‑1 (MQ=255)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          <  1:6062210/50‑1 (MQ=255)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          <  1:224962/50‑1 (MQ=255)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGgaga        <  1:1704475/50‑1 (MQ=255)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          >  1:1803043/1‑48 (MQ=39)
                                      tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       >  1:670193/1‑50 (MQ=255)
                                        tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:4566166/1‑49 (MQ=39)
                                        tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:1456094/1‑49 (MQ=39)
                                        tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:6470858/1‑49 (MQ=39)
                                         ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:1624976/1‑50 (MQ=255)
                                         ttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCa    >  1:2058643/1‑50 (MQ=255)
                                          tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   >  1:2324476/1‑50 (MQ=255)
                                           atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:6248869/1‑50 (MQ=255)
                                           atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:5649890/1‑50 (MQ=255)
                                                |                                             
CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801024‑3801117

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