Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 193,289 A→G H86R (CAC→CGC)  rseP → sigma E protease regulator RseP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273193,2890AG87.0% 139.2 / 14.0 54H86R (CAC→CGC) rsePsigma E protease regulator RseP
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base G (30/17);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AATAGGGTTTGGTAAGGCGCTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGT  >  NZ_CP009273/193155‑193424
                                                                                                                                      |                                                                                                                                       
aaTAGGGTTTGGTAAGGCGCTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCgccgcc                                                                                                                                       <  2:931569/137‑1 (MQ=255)
     ggTTTGGTAAGGCGCTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCtt                                                                                                                                >  1:963698/1‑139 (MQ=255)
     ggTTTGGTAAGGCGCTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCtt                                                                                                                                >  2:1074981/1‑139 (MQ=255)
     ggTTTGGTAAGGCGCTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCACCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCtt                                                                                                                                >  2:814892/1‑139 (MQ=255)
      gTTTGGTAAGGCGCTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTc                                                                                                                               <  2:963698/139‑1 (MQ=255)
       tttGGTAAGGCGCTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTc                                                                                                                               >  1:959647/1‑138 (MQ=255)
          ggTAAGGCGCTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCtt                                                                                                                                >  2:382368/1‑134 (MQ=255)
          ggTAAGGCGCTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCtt                                                                                                                                <  1:382368/134‑1 (MQ=255)
          ggTAAGGCGCTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCaat                                                                                                                            >  2:480552/1‑138 (MQ=255)
                gcgcTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAAt                                                                                                                     >  1:1131209/1‑139 (MQ=255)
                gcgcTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAAt                                                                                                                     >  2:1189163/1‑139 (MQ=255)
                  gcTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCt                                                                                                                   <  2:530367/139‑1 (MQ=255)
                      tGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCg                                                                                                               >  2:1222045/1‑139 (MQ=255)
                        gcggcgAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGc                                                                                                             >  2:1067768/1‑139 (MQ=255)
                           gcgAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTc                                                                                                                               <  2:544951/118‑1 (MQ=255)
                           gcgAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTc                                                                                                                               >  1:544951/1‑118 (MQ=255)
                                          aggagccGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCCCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGattat                                                                                            <  2:647016/133‑1 (MQ=255)
                                             cACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGcc                                                                                        >  2:223939/1‑139 (MQ=255)
                                                  aaTATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAgg                                                                                   >  1:232993/1‑139 (MQ=255)
                                                  aaTATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAgg                                                                                   >  2:77211/1‑139 (MQ=255)
                                                           tCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGt                                                                                                                 >  2:635578/1‑100 (MQ=255)
                                                           tCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGt                                                                                                                 <  1:635578/100‑1 (MQ=255)
                                                           ccGCCCTGATCCAGTTGGGCGGTTATGTCAAAATGCTGGATGCGCGCGCAGACCCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTg                                                                           <  2:291381/137‑1 (MQ=255)
                                                              cccTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCaaa                                                                       >  2:1001801/1‑139 (MQ=255)
                                                              cccTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCaaa                                                                       >  1:1162588/1‑139 (MQ=255)
                                                                   aTCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCa                                                                  >  2:153184/1‑139 (MQ=255)
                                                                        gTTGGGCGGTTATGTCAAAATGCTGGATGAGTGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCAttttt                                                             >  2:732522/1‑139 (MQ=255)
                                                                         ttGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTg                                                            >  1:416078/1‑139 (MQ=255)
                                                                         ttGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTg                                                            <  1:1222045/139‑1 (MQ=255)
                                                                             gCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTc                                                                                                                               <  1:747695/68‑1 (MQ=255)
                                                                             gCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTc                                                                                                                               >  2:747695/1‑68 (MQ=255)
                                                                              cGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATc                                                       >  2:583699/1‑139 (MQ=255)
                                                                                   aTGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGTCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGc                                                  >  1:1279855/1‑139 (MQ=255)
                                                                                   aTGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGc                                                  >  1:118923/1‑139 (MQ=255)
                                                                                        aaaaTGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTAct                                             >  1:878389/1‑139 (MQ=255)
                                                                                         aaaTGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTc                                                                                                                               <  2:692443/56‑1 (MQ=255)
                                                                                         aaaTGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTc                                                                                                                               >  1:692443/1‑56 (MQ=255)
                                                                                           aTGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTActggc                                          <  1:637033/139‑1 (MQ=255)
                                                                                            tGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTActggct                                         <  1:1001801/139‑1 (MQ=255)
                                                                                               tGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGt                                      >  1:771431/1‑139 (MQ=255)
                                                                                                 gATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCg                                                                                       >  2:589391/1‑88 (MQ=255)
                                                                                                 gATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCg                                                                                       <  1:589391/88‑1 (MQ=255)
                                                                                                  aTGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGtttt                                   <  2:1015984/139‑1 (MQ=255)
                                                                                                     aGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttat                                <  2:771431/139‑1 (MQ=255)
                                                                                                       cgcgcAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGc                                                  >  1:1286268/1‑119 (MQ=255)
                                                                                                       cgcgcAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGc                                                  <  2:1286268/119‑1 (MQ=255)
                                                                                                        gcgcAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTtattat                             >  1:2615/1‑139 (MQ=255)
                                                                                                         cgcAGAACCGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTattatt                            >  1:82743/1‑139 (MQ=255)
                                                                                                                 cGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGcc                    <  2:878389/139‑1 (MQ=255)
                                                                                                                 cGGTCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGcc                    <  1:995118/139‑1 (MQ=255)
                                                                                                                  ggTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCt                   >  2:1257230/1‑139 (MQ=255)
                                                                                                                    tCGTTCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTgg                 <  1:153184/139‑1 (MQ=255)
                                                                                                                        tCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGcc                    >  2:883923/1‑132 (MQ=255)
                                                                                                                        tCCGGAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGcc                    <  1:883923/132‑1 (MQ=255)
                                                                                                                            gAACTCCGCCGCCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTc         <  2:833526/139‑1 (MQ=255)
                                                                                                                                  cgccgcCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTActggc                                          <  2:855959/100‑1 (MQ=255)
                                                                                                                                  cgccgcCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTActggc                                          >  1:855959/1‑100 (MQ=255)
                                                                                                                                  cgccgcCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttatta                              >  1:1010641/1‑112 (MQ=255)
                                                                                                                                  cgccgcCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTttatta                              <  2:1010641/112‑1 (MQ=255)
                                                                                                                                  cgccgcCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCggtgg   >  1:305410/1‑139 (MQ=255)
                                                                                                                                  cgccgcCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCggtgg   >  1:622007/1‑139 (MQ=255)
                                                                                                                                  cgccgcCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCggtgg   >  1:1073768/1‑139 (MQ=255)
                                                                                                                                   gccgccATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCggtggt  <  1:814892/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
AATAGGGTTTGGTAAGGCGCTCTGGCGGCGAACTGATAAGCTCGGCACCGAATATGTTATCGCCCTGATCCCGTTGGGCGGTTATGTCAAAATGCTGGATGAGCGCGCAGAACCGGTCGTTCCGGAACTCCGCCACCATGCCTTCAATAATAAATCTGTCGGCCAACGAGCGGCGATTATTGCCGCAGGTCCGGTTGCAAACTTCATTTTTGCTATCTTTGCCTACTGGCTGGTTTTTATTATTGGTGTGCCTGGCGTACGTCCGGTGGT  >  NZ_CP009273/193155‑193424

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