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.93.1% 408.1 / 31.2 102coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base . (50/45);  total (52/50)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 129 total aligned reads displayed.

CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801024‑3801117
                                                |                                             
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              <  1:4355281/50‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              <  1:4670032/50‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              <  1:5461444/50‑1 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              >  1:7104037/1‑50 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              >  1:172959/1‑50 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              >  1:7985049/1‑50 (MQ=255)
cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               <  1:4076237/49‑1 (MQ=255)
 tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  1:7898600/1‑50 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttt                                           <  1:2941774/50‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                            <  1:1406266/49‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                            <  1:1743942/49‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:6676256/48‑1 (MQ=255)
   taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:101188/48‑1 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:3392689/50‑2 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:5598011/50‑2 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:3607976/50‑2 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          <  1:8066862/49‑1 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          >  1:6643517/1‑49 (MQ=255)
     aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          <  1:5070345/49‑1 (MQ=255)
      ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                        <  1:2515969/50‑3 (MQ=255)
      ttCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttttttaa                                        <  1:6155682/50‑3 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      >  1:752743/1‑49 (MQ=255)
        ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      >  1:1048120/1‑49 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:1231531/49‑2 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:1314126/49‑2 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:1145124/50‑1 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:5721307/50‑1 (MQ=255)
         cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    <  1:6094059/50‑1 (MQ=255)
           gagaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   >  1:5983004/1‑49 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:3848269/48‑1 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:1466446/48‑1 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:1929291/48‑1 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:7963229/48‑1 (MQ=39)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 <  1:3352117/50‑1 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:815871/1‑50 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:6165755/1‑50 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 <  1:4033242/50‑1 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 <  1:538303/50‑1 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:1041197/1‑50 (MQ=255)
            agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  1:2375945/49‑1 (MQ=39)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:3159831/50‑1 (MQ=255)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:2276844/50‑1 (MQ=255)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:8272786/1‑49 (MQ=39)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:7768228/1‑49 (MQ=39)
             gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaa                                 >  1:1943315/1‑49 (MQ=39)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                >  1:311136/1‑49 (MQ=39)
              atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                >  1:8057733/1‑49 (MQ=39)
               tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              >  1:7178841/1‑50 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:3539679/1‑50 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:7162620/1‑50 (MQ=255)
                ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:4108449/1‑50 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:2177490/1‑50 (MQ=255)
                 taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            <  1:4178073/50‑1 (MQ=255)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          <  1:4585586/50‑1 (MQ=255)
                   atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          <  1:733964/50‑1 (MQ=255)
                    tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         >  1:4222151/1‑50 (MQ=255)
                    tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         >  1:4041/1‑50 (MQ=255)
                    tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         >  1:4363419/1‑50 (MQ=255)
                     gatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        >  1:6422623/1‑50 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       >  1:4366376/1‑50 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       <  1:5803755/50‑1 (MQ=255)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        <  1:4988246/49‑1 (MQ=39)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        <  1:2392337/49‑1 (MQ=39)
                      atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        <  1:6536651/49‑1 (MQ=39)
                       tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      <  1:210125/50‑1 (MQ=255)
                       tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:4577332/1‑50 (MQ=255)
                        gTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:6466872/1‑50 (MQ=255)
                         ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:7231960/1‑49 (MQ=39)
                         ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:277696/1‑49 (MQ=39)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   <  1:3395842/49‑1 (MQ=39)
                           gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  >  1:4361497/1‑50 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 <  1:7414472/50‑1 (MQ=255)
                            cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:1843027/49‑1 (MQ=39)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                <  1:2422196/50‑1 (MQ=255)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                <  1:4834991/50‑1 (MQ=255)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                <  1:1955791/50‑1 (MQ=255)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                >  1:6771230/1‑50 (MQ=255)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                >  1:1623790/1‑50 (MQ=255)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                <  1:607184/50‑1 (MQ=255)
                             aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                >  1:4293574/1‑50 (MQ=255)
                              cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               <  1:7731961/50‑1 (MQ=255)
                              cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               <  1:1482585/50‑1 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             <  1:5234301/49‑1 (MQ=39)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            <  1:6059190/50‑1 (MQ=255)
                                 aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:2986997/1‑50 (MQ=255)
                                  aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:612972/1‑49 (MQ=39)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:1049622/1‑48 (MQ=39)
                                   aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          <  1:4523221/50‑1 (MQ=255)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          >  1:7203304/1‑48 (MQ=39)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          >  1:5939068/1‑48 (MQ=39)
                                     aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          >  1:2946904/1‑48 (MQ=39)
                                      tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       <  1:5355838/50‑1 (MQ=255)
                                      tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       <  1:3751207/50‑1 (MQ=255)
                                      tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       <  1:5707007/50‑1 (MQ=255)
                                        tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:6439952/1‑49 (MQ=39)
                                        tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:8231383/1‑49 (MQ=39)
                                          tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   >  1:4718749/1‑50 (MQ=255)
                                          tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   >  1:8235750/1‑50 (MQ=255)
                                           atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:7824438/1‑50 (MQ=255)
                                           atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:1186803/1‑50 (MQ=255)
                                                |                                             
CTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801024‑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: