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.94.2% 696.0 / 51.7 173coding (928/1017 nt)waaJUDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (5/5);  new base . (80/82);  total (86/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 211 total aligned reads displayed.

GCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATC  >  NC_000913/3801023‑3801117
                                                 |                                             
gCTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCt                                               >  1:108920/1‑50 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              <  1:354417/50‑1 (MQ=255)
 cTATAATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtt                                              <  1:7403950/50‑1 (MQ=255)
  tataATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             >  1:4497747/1‑50 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                            <  1:6283338/49‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttt                                            <  1:1873097/49‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:5839128/48‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:5575203/48‑1 (MQ=255)
    taATTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCttt                                             <  1:2040287/48‑1 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttta                                         <  1:6899409/50‑2 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          <  1:2503666/49‑1 (MQ=255)
      aTTCCTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCtttttt                                          >  1:1454415/1‑49 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      >  1:3778367/1‑49 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      >  1:2183956/1‑49 (MQ=255)
         ccTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      >  1:6540254/1‑49 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:5930749/49‑2 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTaaat                                      <  1:2583084/49‑2 (MQ=255)
          cTGAGATATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAAtt                                    >  1:7648773/1‑50 (MQ=255)
            gagaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   >  1:6237343/1‑49 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:7854613/48‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:509650/48‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTc                                   <  1:6482925/48‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  1:4999738/49‑1 (MQ=39)
             agaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCa                                  <  1:6696771/49‑1 (MQ=39)
              gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:2299573/50‑1 (MQ=255)
              gaTATAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                <  1:1980362/50‑1 (MQ=255)
               atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                               <  1:5888067/50‑1 (MQ=255)
               atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaata                               <  1:4353525/50‑1 (MQ=255)
               atatAATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaat                                >  1:1057656/1‑49 (MQ=39)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              >  1:7632602/1‑50 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              <  1:6326452/50‑1 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              <  1:4323445/50‑1 (MQ=255)
                tataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataa                              <  1:6939669/50‑1 (MQ=255)
                 ataATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             >  1:1251396/1‑50 (MQ=255)
                  taATGATGTTGCACTAACAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:4963359/1‑50 (MQ=255)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             <  1:945028/49‑1 (MQ=39)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCaataat                             <  1:1625090/49‑1 (MQ=39)
                  taATGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            <  1:6438436/50‑1 (MQ=255)
                   aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:1101424/1‑49 (MQ=39)
                   aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:7598762/1‑49 (MQ=39)
                   aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:7528058/1‑49 (MQ=39)
                   aaTGATGTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAAtt                            >  1:5739677/1‑49 (MQ=39)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:3986991/49‑1 (MQ=39)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:8258419/49‑1 (MQ=39)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:4461607/49‑1 (MQ=39)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTg                           <  1:3283432/49‑1 (MQ=39)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          >  1:1369566/1‑50 (MQ=255)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          <  1:2115729/50‑1 (MQ=255)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          >  1:5995184/1‑50 (MQ=255)
                    atgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGa                          <  1:2373046/50‑1 (MQ=255)
                     tgatgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAt                         >  1:3945440/1‑50 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAttt                       <  1:611908/50‑1 (MQ=255)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        <  1:7964937/49‑1 (MQ=39)
                       atgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGAtt                        <  1:7100180/49‑1 (MQ=39)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:6254514/1‑50 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:4260783/1‑50 (MQ=255)
                        tgTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTc                      >  1:6694647/1‑50 (MQ=255)
                         gTTGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:6644550/1‑50 (MQ=255)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:3265750/1‑49 (MQ=39)
                          ttGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCg                     >  1:120503/1‑49 (MQ=39)
                           tGCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   >  1:1531317/1‑50 (MQ=255)
                            gCACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCa                   <  1:6724959/49‑1 (MQ=39)
                             cACTAGAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 >  1:6601219/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtc                 >  1:3153492/1‑50 (MQ=255)
                             cACTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAt                  <  1:1795205/49‑1 (MQ=39)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                <  1:5027125/50‑1 (MQ=255)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                >  1:3381394/1‑50 (MQ=255)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                <  1:4190704/50‑1 (MQ=255)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                <  1:2480678/50‑1 (MQ=255)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                >  1:505578/1‑50 (MQ=255)
                              aCTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtct                <  1:5335749/50‑1 (MQ=255)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               >  1:920487/1‑50 (MQ=255)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               >  1:7351240/1‑50 (MQ=255)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               >  1:7258069/1‑50 (MQ=255)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               <  1:562417/50‑1 (MQ=255)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               <  1:5427117/50‑1 (MQ=255)
                               cTAAAAGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCAtctc               <  1:2953735/50‑1 (MQ=255)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCg              <  1:2857947/48‑1 (MQ=39)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGt             <  1:243959/49‑1 (MQ=39)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:7994952/1‑50 (MQ=255)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:6822206/1‑50 (MQ=255)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            <  1:1391785/50‑1 (MQ=255)
                                  aaaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:832029/1‑50 (MQ=255)
                                   aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:4171065/1‑49 (MQ=39)
                                   aaGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:3894446/1‑49 (MQ=39)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:4850003/1‑48 (MQ=39)
                                    aGATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          <  1:807860/50‑1 (MQ=255)
                                    aGATGTTTATATATTGTTTAAATTCAATAATTGATTTCGCATCTCGTg            >  1:1755573/1‑48 (MQ=37)
                                     gATGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGgag         <  1:2246886/50‑1 (MQ=255)
                                      aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGgaga        <  1:7680107/50‑1 (MQ=255)
                                      aTGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGga          >  1:4655156/1‑48 (MQ=39)
                                       tGTTTATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGagag       >  1:5788581/1‑50 (MQ=255)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:757047/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:1667869/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGt      >  1:4739047/1‑49 (MQ=39)
                                         tttATATCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTc     <  1:1336813/50‑1 (MQ=255)
                                           tataTCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCAt   >  1:4889107/1‑50 (MQ=255)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:363038/1‑50 (MQ=255)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:4421396/1‑50 (MQ=255)
                                            atatCTTTTTTAAATTCAATAATTGATTTCGCATCTCGTGGAGAGTCATc  >  1:6360775/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: