Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,129,394 G→C 7.9% intergenic (+20/+20) murJ → / ← flgN putative lipid II flippase/export chaperone for FlgK and FlgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,129,3940GC 7.9% 128.0 / 3.2 50intergenic (+20/+20)murJ/flgNputative lipid II flippase/export chaperone for FlgK and FlgL
Reads supporting (aligned to +/- strand):  ref base G (27/19);  new base C (2/2);  total (29/21)
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.78e-01

ATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCC  >  NC_000913/1129248‑1129540
                                                                                                                                                  |                                                                                                                                                  
aTTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtg                                                                                                                                                  >  1:890839/1‑149 (MQ=255)
    ggTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAgg                                                                                                                                              <  1:426815/149‑1 (MQ=255)
    ggTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTGCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAgg                                                                                                                                              <  1:422827/149‑1 (MQ=255)
        ccATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCgg                                                                                                                                          <  2:323770/149‑1 (MQ=255)
         cATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGga                                                                                                                                         >  1:790038/1‑149 (MQ=255)
         cATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCGga                                                                                                                                         <  2:861948/149‑1 (MQ=255)
         cATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAAGTTGCCCGCCGGCCGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGga                                                                                                                                         >  1:559180/1‑149 (MQ=255)
         cATGCCCTGGCGTTTACTGCGTTTAATGGCGGGCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTAGCCCGCCGGACGGGGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGga                                                                                                                                         <  2:138715/149‑1 (MQ=255)
          aTGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGTTTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCGgaa                                                                                                                                        <  1:204547/149‑2 (MQ=255)
           tGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGgaga                                                                                                                                       >  2:518823/1‑149 (MQ=255)
               cTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAAt                                                                                                                                   >  2:357364/1‑149 (MQ=255)
                      ttACTGCGTTTAATGGAGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTAATGGGCTTCAAAGTTATAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGACGGTGATAATCTTCAGa                                                                                                                            >  2:674095/1‑149 (MQ=255)
                             gTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAAt                                                                                                                     >  2:965965/1‑149 (MQ=255)
                                 aaTGGCGGTCGTGCTGGGTTGGATTGCAGCGTACTTCGCTACACTGGCGGTACTTGGCTTCATATTTATAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCttt                                                                                                                 <  1:197300/149‑1 (MQ=255)
                                   tGGCGGTCGTGCTGGCTGGGATTGCCGCGTACTTCGCTGCTCTTGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGACGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCttttt                                                                                                               <  1:158402/149‑1 (MQ=255)
                                    ggCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTa                                                                                                              >  2:435470/1‑149 (MQ=255)
                                     gCGGTGGTGCAGGCGGGGTTTGCCGAGGACTGCGCTGCACTGGAGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCGGGCCTGCAGTGCAGGCCGGAGATAATCTGCAGATCGAAATCATTTTAc                                                                                                             <  2:190656/149‑1 (MQ=255)
                                             gCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTTCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCCTAATCTTTTTACCGCCGCGa                                                                                                     >  2:673768/1‑149 (MQ=255)
                                             gCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACTATGCATTCCGGCCTGCACTGCAGGCCGGAGAGTATCTTCAGAACGAAATCTTTTTACCGACGCGa                                                                                                     >  1:267395/1‑149 (MQ=255)
                                               tGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTCCTGGGCTTCTAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCGGCAGTGCAGGCCGGAGATAATCTTCAGATCGAACTCTTTTTACCGCCGCGATg                                                                                                   >  1:465834/1‑149 (MQ=255)
                                                 gCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCTGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                                                                                                 >  1:800768/1‑149 (MQ=255)
                                                      gATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgttgtt                                                                                            <  1:375980/149‑1 (MQ=255)
                                                         tGCCGCGTATTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAGAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGGGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGGTGTTGAg                                                                                         <  2:299357/149‑1 (MQ=255)
                                                           ccGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTTCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGCCCGGAGATAATCATCAGATCTAAATCTTTTTACCGACGCGATGGGTTCTTGAgtt                                                                                       >  2:696998/1‑147 (MQ=255)
                                                                 aCTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTCCCGCCGCGATGGGTTGTTGAGGTCTGAcc                                                                                 >  1:582394/1‑149 (MQ=255)
                                                                       cTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGc                                                                           >  1:2203/1‑149 (MQ=255)
                                                                           aCTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATATTCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGaccaa                                                                       >  2:616024/1‑147 (MQ=255)
                                                                           aCTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCAGGAGATAATCTTCAGATCGAAATCTTTTTACCGACGCGATGGGTTGTTGAGGTCTGAACGTTCGACCCa                                                                       >  1:744637/1‑149 (MQ=255)
                                                                            cTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACTGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCat                                                                      <  1:389174/149‑1 (MQ=255)
                                                                             tGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTTACCGTTCGCCCCata                                                                     >  1:922159/1‑149 (MQ=255)
                                                                                    aCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTc                                                              <  1:965965/149‑1 (MQ=255)
                                                                                    aCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCACATCGACATCGTTTTACCGCAGCGATCGGTTGTTGATGTATGACCATTCGCCCCATAAAGCCTc                                                              >  1:236573/1‑149 (MQ=255)
                                                                                    aCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATACCGGTCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTGGTTGAGGTCTGACCGTTCGCCCCATATAGCGTc                                                              <  1:918197/149‑1 (MQ=255)
                                                                                       gggCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGAGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGGCGGt                                                           >  1:553666/1‑149 (MQ=255)
                                                                                         gCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGAATTCCGGCCTGCAGTGCAGGACGGAGATAAACTTCAGATCGAAAACTTTTTACAGCCGCGATGGGTTGTTGAGGTATGACCTTTAGACACATATATCGTCGGTTc                                                         >  2:4967/1‑149 (MQ=255)
                                                                                         gCTTCAAAGGAAAAGAAGTTGCCCGCCGGACGGGGTATCACTGCATTCCGGCCTGCAGTGCAGGCCGGCGATAAGCTGCAGATCGAAAGCTTTTTACCGCCGAGACGCGGAGTTGAGGTCTGACAGTTCGCCCGTAATAGCGTCGGtgc                                                         >  2:668250/1‑147 (MQ=255)
                                                                                              aaaGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGGTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGAt                                                    <  1:123326/149‑1 (MQ=255)
                                                                                                 gTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCg                                                 >  2:396027/1‑149 (MQ=255)
                                                                                                   tAAAGAATTTGCCCGCCGGACGGTTTAACAATGCATTCCGGCCTGCAGTGCAGGCATGAGATAATCTTCACATCGAAATCTTTTTACCGCCTATATGGGTTGTTGAGTTCTGACCGTTAGCCCCATATAGAGTCGGAGACTGATgagtt                                               >  2:689131/1‑145 (MQ=255)
                                                                                                    aaaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGCGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGtt                                              >  2:285005/1‑149 (MQ=255)
                                                                                                       gAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCa                                           >  1:774373/1‑149 (MQ=255)
                                                                                                             gCCCGCCGGACGGGGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACAttt                                     <  1:47061/149‑1 (MQ=255)
                                                                                                                cgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCa                                  >  1:293285/1‑149 (MQ=255)
                                                                                                                     gACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGcc                             >  1:39455/1‑149 (MQ=255)
                                                                                                                        ggTGTCACAATGCATTCAGGCCTGCAGTGGAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGTTCTGACCGTTGGCCCCATTTAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGt                          <  2:766486/149‑1 (MQ=255)
                                                                                                                           gtAACAATGCATTCCGGCCTGCAGGGCTGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGGTGa                       <  2:465834/149‑1 (MQ=255)
                                                                                                                               cAATGCATTCCGGCCTGTAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTCTTGATTgc                   >  1:92197/1‑149 (MQ=255)
                                                                                                                                 aTGCATTCCGGCCTGCAGTGCAGGGCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcgc                 <  2:790038/149‑1 (MQ=255)
                                                                                                                                   gCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTc               <  1:518823/149‑1 (MQ=255)
                                                                                                                                      ttCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAAt            <  2:890839/149‑1 (MQ=255)
                                                                                                                                         cGGCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTg         >  2:628805/1‑149 (MQ=255)
                                                                                                                                               gcagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTc                                            <  2:366332/108‑1 (MQ=255)
                                                                                                                                               gcagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTc                                            >  1:366332/1‑108 (MQ=255)
                                                                                                                                               gcagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTc   >  2:888982/1‑149 (MQ=255)
                                                                                                                                                cagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTCCCGCCGCGATGGGTTGTTGAGGTCTGACCGTGCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACAATTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTcc  <  2:591616/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCC  >  NC_000913/1129248‑1129540

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑