Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,950,973 C→A 8.9% intergenic (‑123/+34) fliL ← / ← fliK flagellar protein FliL/flagellar hook‑length control protein FliK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,950,9730CA 8.9% 126.4 / 2.8 52intergenic (‑123/+34)fliL/fliKflagellar protein FliL/flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base C (27/19);  new base A (2/3);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

CTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGC  >  NC_002947/4950832‑4951111
                                                                                                                                             |                                                                                                                                          
cTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCcgcg                                                                                                                                     <  5:4357/149‑1 (MQ=255)
   cTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaag                                                                                                                                  <  5:335203/149‑1 (MQ=255)
   cTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCACgaa                                                                                                                                   <  1:114042/148‑1 (MQ=255)
       gCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagc                                                                                                                              >  2:82772/1‑149 (MQ=255)
           ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCTGTGCTGGCTTCTTAGAGGGCTCGCCCGAGAAGAagccagc                                                                                                                          >  8:153740/1‑149 (MQ=255)
              gCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagccagcc                                                                                                                         >  5:186092/1‑147 (MQ=255)
              gCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAagccagcc                                                                                                                         >  7:269336/1‑147 (MQ=255)
                cATGCCAATAATCCGTCGTCCAGCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGGGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTCCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAGGCAGCCAGACCAgg                                                                                                                     >  6:304724/1‑149 (MQ=255)
                  tGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGcc                                                                                                                   <  1:82772/149‑1 (MQ=255)
                      aataatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCaa                                                                                                                 >  3:82155/1‑147 (MQ=255)
                             cgtcgtCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCa                                                                                                        <  3:356269/149‑1 (MQ=255)
                                   ccATGGGGGGTTCAAAGGAGAATGGGTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCTGGGGATGGATTCCAGCAGGTTTTGCGGGGGGGGCTTCTTCGCGGGCTGGCCCGCGAAGAAGCCAGCCCGGGCCAACAGTTATCAGGCgta                                                                                                  <  1:168024/149‑1 (MQ=255)
                                    cctccggggGTCAACGGCGAATGGCTTGGCTTGCGCAAGTGTAATGGAAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                  <  5:304724/143‑1 (MQ=255)
                                        ggggtTTCAAAGGAGTATGGCTTGGCTTTAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAAt                                                                                             <  7:153740/149‑1 (MQ=255)
                                          gggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATc                                                                                            >  6:105426/1‑148 (MQ=255)
                                           ggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGa                                                                                          <  8:269336/149‑1 (MQ=255)
                                               cAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCg                                                                                       <  6:186092/148‑1 (MQ=255)
                                                aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAc                                                                                                                <  2:312770/122‑1 (MQ=255)
                                                aaaGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAAc                                                                                                                >  1:312770/1‑122 (MQ=255)
                                                  aGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACt                                                                                   >  3:247054/1‑149 (MQ=255)
                                                                   tgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGAc                                                                                         <  1:302091/126‑1 (MQ=255)
                                                                   tgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGAc                                                                                         >  2:302091/1‑126 (MQ=255)
                                                                    gAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCCGCCCCGGCCAACAGTTATGAGGCGGAGTAATAGACCGCACTGTCGCNGACAACCNCctg                                                                 >  8:33397/1‑149 (MQ=255)
                                                                      gCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgct                                                               >  7:290217/1‑149 (MQ=255)
                                                                       cAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCGGGCGAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACctgctg                                                              >  4:237509/1‑149 (MQ=255)
                                                                          gtgtTATGCCAGATTTTGCAGGGGAGGGATTCCAGCAGGTTTTGCGGGGCGGGCGTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTc                                                           <  4:330544/149‑1 (MQ=255)
                                                                           tgtTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCCCGCCAACAGTTATCAGGCGTAGTAATCGACCGCCCTGTCGcc                                                                             >  2:10155/1‑130 (MQ=255)
                                                                           tgtTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTc                                                           >  5:162365/1‑148 (MQ=255)
                                                                           tgtTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGcc                                                                             <  1:10155/130‑1 (MQ=255)
                                                                               aTGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCCGGCTCGCCCGCCAAAAAGCCAGCCCAGGCCAACAGGTATAAGGCGTAGTAATAGACCGCACTGTCGCCGCCAACAACATGCTGCTCCAAAg                                                      >  3:21569/1‑149 (MQ=255)
                                                                                tGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAg                                                      <  4:82155/148‑1 (MQ=255)
                                                                                  ccAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCAAAGTCGCTGACAACCCACTGCTGCTCCGCAGGcc                                                   >  4:233923/1‑149 (MQ=255)
                                                                                   cAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCa                                                                                     <  6:70898/114‑1 (MQ=255)
                                                                                   cAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCa                                                                                     >  5:70898/1‑114 (MQ=255)
                                                                                             aGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCg                                        >  3:342183/1‑149 (MQ=255)
                                                                                                    ggatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGgc                                               <  6:215508/135‑1 (MQ=255)
                                                                                                    ggatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGgc                                               >  5:215508/1‑135 (MQ=255)
                                                                                                      atTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGCAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCg                               >  7:105671/1‑149 (MQ=255)
                                                                                                        tCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGc                             >  2:14094/1‑149 (MQ=255)
                                                                                                          cagcagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGTccc                           <  8:105671/149‑1 (MQ=255)
                                                                                                          cagcagGTTTTGCGGTGCTGGCTTCTTCGCGGGCGCGCCCGCGAAGAAGCCCGCCCAGGCCAACAGGTTTAAGGCGTAGTAATAGACCCCACAGTCGCCGACAACCACCGGCGGGTGCAAAGGCCCGGCGCTATCGACCGCTTCGGccc                           >  4:223590/1‑149 (MQ=255)
                                                                                                           agcagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGcccc                          >  1:5232/1‑149 (MQ=255)
                                                                                                            gcagGTTTTGCGGTGCGGGCTTCTTCGCGGGCTCGCCCGCTAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGcc                                                   <  4:198081/123‑1 (MQ=255)
                                                                                                             cagGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACAGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCGATACACCGCTTCGGCCCCAc                        >  4:30914/1‑149 (MQ=255)
                                                                                                               ggTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCAccc                      <  4:342183/149‑1 (MQ=255)
                                                                                                                    tGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTc                 >  7:252190/1‑149 (MQ=255)
                                                                                                                        gTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGcc             <  7:33397/149‑1 (MQ=255)
                                                                                                                            tGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCt                                  >  7:288449/1‑124 (MQ=255)
                                                                                                                            tGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCt                                  <  8:288449/124‑1 (MQ=255)
                                                                                                                              gCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcg        <  5:105426/148‑1 (MQ=255)
                                                                                                                                  cttcGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgccg   >  5:17503/1‑149 (MQ=255)
                                                                                                                                   ttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgccgc  >  4:267988/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                          
CTGCTTCGCTCTTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCCGC  >  NC_002947/4950832‑4951111

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