Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,798,235 A→T 23.3% intergenic (+28/+4) waaL → / ← waaU O‑antigen ligase/lipopolysaccharide core biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,798,2350AT23.3% 127.7 / 41.3 61intergenic (+28/+4)waaL/waaUO‑antigen ligase/lipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base A (23/23);  new base T (7/7);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01

TGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATT  >  NC_000913/3798116‑3798378
                                                                                                                       |                                                                                                                                               
tgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCAct                                                                                                                                         >  1:303346/1‑128 (MQ=255)
tgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCAct                                                                                                                                         <  2:303346/128‑1 (MQ=255)
                  gcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                                                                >  4:166481/1‑103 (MQ=255)
                  gcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                                                                <  3:166481/103‑1 (MQ=255)
                          aTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                                                                <  1:816891/95‑1 (MQ=255)
                          aTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                                                                >  2:816891/1‑95 (MQ=255)
                                 ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActt                                                                                                                                        <  4:870564/96‑1 (MQ=255)
                                 ttATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActt                                                                                                                                        >  3:870564/1‑96 (MQ=255)
                                  tataTCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCa                                                                                                                                           <  4:344223/92‑1 (MQ=255)
                                  tataTCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCa                                                                                                                                           >  3:344223/1‑92 (MQ=255)
                                                 ccTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                                                                >  1:551601/1‑72 (MQ=255)
                                                 ccTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                                                                <  2:551601/72‑1 (MQ=255)
                                                    ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                                                                           <  4:620309/74‑1 (MQ=255)
                                                    ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                                                                           >  3:620309/1‑74 (MQ=255)
                                                    ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCt                                                                                                                                 <  1:645950/84‑1 (MQ=255)
                                                    ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCt                                                                                                                                 >  2:645950/1‑84 (MQ=255)
                                                         tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTGAATCTCTTTACTGAAtt                                                                                                       <  1:328577/105‑1 (MQ=255)
                                                         tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAAAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAtt                                                                                                       >  2:328577/1‑105 (MQ=255)
                                                          cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg                                                           >  1:428833/1‑148 (MQ=255)
                                                               ttAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct                                                                                                                                         >  3:949909/1‑65 (MQ=255)
                                                               ttAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct                                                                                                                                         <  4:949909/65‑1 (MQ=255)
                                                                               tacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTat                                      <  2:428833/148‑1 (MQ=255)
                                                                                acaatTAATTAAGAATAAACAAGTTTAATAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaa                                                                                                 <  3:662087/88‑1 (MQ=255)
                                                                                acaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTACAACTCACTTCTTATCTATACAACTCAATCTCTTTACTGAATTAGTTaa                                                                                                 >  4:662087/1‑88 (MQ=255)
                                                                                   attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaa                                                                                                                      <  1:1055290/64‑1 (MQ=255)
                                                                                   attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaa                                                                                                                      >  2:1055290/1‑64 (MQ=255)
                                                                                    ttaattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                                                               >  1:340489/1‑70 (MQ=255)
                                                                                    ttaattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                                                               <  2:340489/70‑1 (MQ=255)
                                                                                         taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGt                                                >  1:505370/1‑128 (MQ=255)
                                                                                         taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGt                                                <  2:505370/128‑1 (MQ=255)
                                                                                           aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt                                                               <  4:794232/111‑1 (MQ=255)
                                                                                           aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt                                                               >  3:794232/1‑111 (MQ=255)
                                                                                              aTAAACAATTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATCGATTTATGATGATTCGGGGAcc                       >  1:966738/1‑148 (MQ=255)
                                                                                              aTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTATTTAAAGTTTc                                                                                           <  3:652696/80‑1 (MQ=255)
                                                                                              aTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                           >  4:652696/1‑80 (MQ=255)
                                                                                                              agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                                                                        >  1:776972/1‑51 (MQ=25)
                                                                                                              agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                                                                        <  2:776972/51‑1 (MQ=25)
                                                                                                              agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat                                                                    <  4:608648/87‑1 (MQ=255)
                                                                                                              agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat                                                                    >  3:608648/1‑87 (MQ=255)
                                                                                                              agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                              <  2:448317/93‑1 (MQ=255)
                                                                                                              agTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                              >  1:448317/1‑93 (MQ=255)
                                                                                                              agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTAt                                                                                      <  3:724985/69‑1 (MQ=255)
                                                                                                              agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTAt                                                                                      >  4:724985/1‑69 (MQ=255)
                                                                                                              agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAggg       <  2:966738/148‑1 (MQ=255)
                                                                                                               gTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa                                       <  4:491734/115‑1 (MQ=255)
                                                                                                               gTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa                                       >  3:491734/1‑115 (MQ=255)
                                                                                                                tGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg                                                                                                     <  1:1091581/52‑1 (MQ=37)
                                                                                                                tGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg                                                                                                     >  2:1091581/1‑52 (MQ=37)
                                                                                                                tGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGat     >  2:790135/1‑148 (MQ=255)
                                                                                                                 gAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGata    >  1:745675/1‑148 (MQ=255)
                                                                                                                   gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATAtt  <  1:790135/148‑1 (MQ=255)
                                                                                                                    ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGa                                                                                   <  1:251775/66‑1 (MQ=255)
                                                                                                                    ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGa                                                                                   >  2:251775/1‑66 (MQ=255)
                                                                                                                    ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt                                                   >  4:105109/1‑98 (MQ=255)
                                                                                                                    ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt                                                   <  3:105109/98‑1 (MQ=255)
                                                                                                                    ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa                                       >  2:758811/1‑110 (MQ=255)
                                                                                                                    ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa                                       <  1:758811/110‑1 (MQ=255)
                                                                                                                     tAAAACTCACTTCTTCTCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCCGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGta                                               <  1:556344/101‑1 (MQ=255)
                                                                                                                     tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtt                                                                             >  3:108751/1‑71 (MQ=255)
                                                                                                                     tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtt                                                                             <  4:108751/71‑1 (MQ=255)
                                                                                                                     tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtcg                                                              <  4:84245/86‑3 (MQ=255)
                                                                                                                     tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                              >  3:84245/1‑86 (MQ=255)
                                                                                                                     tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGta                                               >  2:556344/1‑101 (MQ=255)
                                                                                                                     tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAg                                            <  4:733782/104‑1 (MQ=255)
                                                                                                                     tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAg                                            >  3:733782/1‑104 (MQ=255)
                                                                                                                     tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTa                 <  1:88020/131‑1 (MQ=255)
                                                                                                                     tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTa                 >  2:88020/1‑131 (MQ=255)
                                                                                                                      aaaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg                                                             <  2:817857/86‑1 (MQ=255)
                                                                                                                      aaaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg                                                             >  1:817857/1‑86 (MQ=255)
                                                                                                                       aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                           <  2:267922/55‑1 (MQ=255)
                                                                                                                       aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                           >  1:267922/1‑55 (MQ=255)
                                                                                                                       aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtttt                                                                           <  1:869403/71‑1 (MQ=255)
                                                                                                                       aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtttt                                                                           >  2:869403/1‑71 (MQ=255)
                                                                                                                       aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGt                                                                      <  1:612813/76‑1 (MQ=255)
                                                                                                                       aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGt                                                                      >  2:612813/1‑76 (MQ=255)
                                                                                                                       aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGAttt                                        <  1:760564/106‑1 (MQ=255)
                                                                                                                       aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGAttt                                        >  2:760564/1‑106 (MQ=255)
                                                                                                                       |                                                                                                                                               
TGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATT  >  NC_000913/3798116‑3798378

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