Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,004,644 G→A 6.9% G258D (GGC→GAC)  fliD → flagellar filament capping protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,004,6440GA 6.9% 178.0 / 2.7 58G258D (GGC→GAC) fliDflagellar filament capping protein
Reads supporting (aligned to +/- strand):  ref base G (31/23);  new base A (2/2);  total (33/25)
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.98e-01

GCCAGTGCCAGCAGCAACGGTATGGAGGTCTCGGTTGCCGCCCAGAATGCGCAGCTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAGCTGATAGCCAGAGTT  >  NC_000913/2004496‑2004790
                                                                                                                                                    |                                                                                                                                                  
gCCAGTGCCAGCAGCAACGGTATGGAGGTCTCGGTTGCCGCCCAGAATGCGCAGCTGACAGTCAACAACGTCGCCATCGAGAACCGCAGCAACACCATCCGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACggg                                                                                                                                                    >  1:258570/1‑149 (MQ=255)
     tGCCAGCAGCAACGGTATGGAGGTCTCGGTTGCCGCCCAGAATGCGCAGCTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAAcc                                                                                                                                               <  4:120838/149‑1 (MQ=255)
            agcaACGGTATGGAGGTCTCGGTTGCCGCCCAGAATGCGCAGCTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCCGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGc                                                                                                                                         >  5:194903/1‑148 (MQ=255)
             gcaACGGTATGGAGGTCTCGGTTGCCGCCCAGAATGCGCAGCTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCt                                                                                                                                        >  7:244143/1‑148 (MQ=255)
               aaCGGTATGGAGGTCTCGGTTGCCGCCCAGAATGCGCAGCTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAAc                                                                                                                                     >  1:141109/1‑149 (MQ=255)
                  ggTATGGAGGTCTCGGTTGCCGCCCAGAATGCGCAGCTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCAt                                                                                                                                  <  6:177770/149‑1 (MQ=255)
                                             aaTGCGCAGCTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAAc                                                                                                       <  7:128112/149‑1 (MQ=255)
                                                  gcAGCTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGACAACCAGACGCTAACCATCACTCAGGac                                                                                                                        <  8:228134/127‑1 (MQ=255)
                                                  gcAGCTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGACAACCAGACGCTAACCATCACTCAGGac                                                                                                                        >  7:228134/1‑127 (MQ=255)
                                                   cAGCTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCCCCCTGAACCTGAACGATGTCACCCCGGGCCACCAGACGCTAACCATCACTCAGGACCCCTCCAAAGCGCAACCGGCGag                                                                                                 >  4:149635/1‑148 (MQ=255)
                                                    aGCGGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGAtt                                                                                                <  3:149635/149‑1 (MQ=255)
                                                      cTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTa                                                                                               >  1:161638/1‑148 (MQ=255)
                                                         aCAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTaaa                                                                                             <  2:161638/147‑1 (MQ=255)
                                                         aCAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTaaa                                                                                             <  2:141109/147‑1 (MQ=255)
                                                           aGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGATAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACt                                                                                         >  8:332901/1‑149 (MQ=255)
                                                            gTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTACCCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTg                                                                                        >  3:125670/1‑149 (MQ=255)
                                                               aacaacGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGt                                                                                     <  5:310633/149‑1 (MQ=255)
                                                               aacaacGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGt                                                                                     >  4:311542/1‑149 (MQ=255)
                                                               aacaacGTCGCAATCGATCAAAGCAGCAACAACATCAGCGACGCGCTGGAAAACATCAACCTGAACCTTAACTATGTCAACACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGt                                                                                     <  6:194903/149‑1 (MQ=255)
                                                                 caacGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTg                                                                                    >  4:46561/1‑148 (MQ=255)
                                                                   acGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAAt                                                                                 >  2:50752/1‑149 (MQ=255)
                                                                   acGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAAt                                                                                 <  8:244143/149‑1 (MQ=255)
                                                                     gTCGCCATCGAGAATATCAGCAACACCATCAGCTACGCGCTGGAAAAAATCACCCTGAACCTGAACTATGTCACCACGGGCAACCAGCCGCTAACCATCACACAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGc                                                                               >  4:8874/1‑149 (MQ=255)
                                                                      tCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAA‑CAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCt                                                                             <  8:213144/149‑1 (MQ=255)
                                                                              gagaACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTc                                                                      >  6:63244/1‑149 (MQ=255)
                                                                                gaACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCg                                                                     >  3:269830/1‑148 (MQ=255)
                                                                                gaACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGc                                                                    <  5:63244/149‑1 (MQ=255)
                                                                                        gcaACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGAt                                                            <  3:46561/149‑1 (MQ=255)
                                                                                          aaCACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATAc                                                          >  6:220578/1‑149 (MQ=255)
                                                                                            cacCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCt                                                        >  4:225154/1‑149 (MQ=255)
                                                                                               cATCAGCGACGCGTTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTca                                                     <  3:322069/149‑1 (MQ=255)
                                                                                               cATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCaa                                                                                                         >  4:254458/1‑97 (MQ=255)
                                                                                               cATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCaa                                                                                                         <  3:254458/97‑1 (MQ=255)
                                                                                               cATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCtt                                                       >  3:168661/1‑147 (MQ=255)
                                                                                                 tCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTcagc                                                   <  7:131050/149‑1 (MQ=255)
                                                                                                        cgcgCTGGAAAACATCACCCTGACCCTGAACGATGTCACCACGGGCACCCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGa                                            >  8:274460/1‑149 (MQ=255)
                                                                                                          cgcTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGAcc                                          <  2:258570/149‑1 (MQ=255)
                                                                                                                       cACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACAccgccg                             >  8:55070/1‑149 (MQ=255)
                                                                                                                        aCCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACAccgccg                             <  4:125670/148‑1 (MQ=255)
                                                                                                                        aCCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACAccgccg                             >  7:201973/1‑148 (MQ=255)
                                                                                                                          cctgaacctgaaCGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACAccgccg                             >  6:179196/1‑146 (MQ=255)
                                                                                                                          cctgaacctgaaCGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACAccg                                <  6:253930/143‑1 (MQ=255)
                                                                                                                          cctgaacctgaaCGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACAccg                                >  5:253930/1‑143 (MQ=255)
                                                                                                                                cctgaacGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCgg                    <  7:332901/149‑1 (MQ=255)
                                                                                                                                cctgaacGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCgg                    <  3:311542/149‑1 (MQ=255)
                                                                                                                                 ctgaacGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCgg                    >  8:282431/1‑148 (MQ=255)
                                                                                                                                   gaacGATGTCACCACGGGCAACCAGACGCTACCCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAg                 >  3:190915/1‑149 (MQ=255)
                                                                                                                                   gaacGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACCCCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAAGGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAAGACACCGCCGTAGATGCGGGAg                 >  4:283986/1‑149 (MQ=255)
                                                                                                                                   gaacGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAg                 >  1:28856/1‑149 (MQ=255)
                                                                                                                                    aacGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAGc                >  8:326547/1‑149 (MQ=255)
                                                                                                                                       gATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTca                                                     >  1:50675/1‑109 (MQ=255)
                                                                                                                                       gATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTca                                                     <  2:50675/109‑1 (MQ=255)
                                                                                                                                         tGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAGCTGATa           <  5:220578/149‑1 (MQ=255)
                                                                                                                                            caccacGGACAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAGCTGATAGc         >  2:131907/1‑148 (MQ=255)
                                                                                                                                             accacGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAGCTGATAGc         >  4:97859/1‑147 (MQ=255)
                                                                                                                                              ccacGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAGCTGATAGc         >  8:230477/1‑146 (MQ=255)
                                                                                                                                               cacGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAGCTGATAGCCaga     <  1:50752/149‑1 (MQ=255)
                                                                                                                                               cacGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAGCTGATAGCCa       >  6:130486/1‑147 (MQ=255)
                                                                                                                                               cacGGACAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAGCTGATAGCCaga     <  1:131907/149‑1 (MQ=255)
                                                                                                                                                  gggCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAGCTGATAGCCAGAGtt  <  7:230477/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCCAGTGCCAGCAGCAACGGTATGGAGGTCTCGGTTGCCGCCCAGAATGCGCAGCTGACAGTCAACAACGTCGCCATCGAGAACAGCAGCAACACCATCAGCGACGCGCTGGAAAACATCACCCTGAACCTGAACGATGTCACCACGGGCAACCAGACGCTAACCATCACTCAGGACACCTCCAAAGCGCAAACGGCGATTAAAGACTGGGTGAATGCCTACAACTCGCTAATAGATACCTTCAGCAGCCTGACCAAATACACCGCCGTAGATGCGGGAGCTGATAGCCAGAGTT  >  NC_000913/2004496‑2004790

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