Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,109 T→C 9.6% I176V (ATC→GTC)  flhA ← flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,936,1090TC 9.6% 141.3 / 5.3 52I176V (ATC→GTC) flhAflagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base C (3/2);  total (27/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 = 8.63e-01

AGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTTGATGATCATGAGGATGGCGAACACCACGGCAC  >  NC_002947/4935961‑4936258
                                                                                                                                                     |                                                                                                                                                     
aGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATc                                                                                                                                                      <  1:496715/150‑1 (MQ=255)
aGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATc                                                                                                                                                      <  2:1357907/150‑1 (MQ=255)
    tGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGc                                                                                                                                                  <  2:383386/150‑1 (MQ=255)
    cgCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGc                                                                                                                                                  <  1:1249905/149‑1 (MQ=255)
        atctaTCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGc                                                                                                                                              >  1:410023/5‑150 (MQ=255)
              cAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAg                                                                                                                                        <  1:640529/150‑1 (MQ=255)
                       gaACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTc                                                                                                                               <  2:370223/150‑1 (MQ=255)
                       gaACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTc                                                                                                                               >  1:910343/1‑150 (MQ=255)
                         aCAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGGACGAACTTGCTGGCACCGTCCATCGCGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTcaa                                                                                                                             >  2:419293/1‑148 (MQ=255)
                         aCAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTg                                                                                                                                          >  2:382595/1‑137 (MQ=255)
                         aCAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTg                                                                                                                                          <  1:382595/137‑1 (MQ=255)
                          cAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGAt                                                                                                                            >  2:605768/1‑150 (MQ=255)
                           aGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGGCGATCAGGCCGGCATTGAGGTCGGAGTCAATg                                                                                                                           >  2:607086/1‑150 (MQ=255)
                            agATCAGCCCGCCGGCTATTGCGTCGCCGCGCACGGACTTGATGACCTCCCGGTCTCCCTCGTAGCACTCTTCCTCATGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATgg                                                                                                                          <  1:525287/149‑1 (MQ=255)
                                 aGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATc                                                                                                                     <  2:402091/150‑1 (MQ=255)
                                  gCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCt                                                                                                                    <  2:188823/150‑1 (MQ=255)
                                      gCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCtt                                                                                                                >  1:1018842/1‑150 (MQ=255)
                                                       cgcgCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGa                                                                                               >  2:1190848/1‑150 (MQ=255)
                                                        gcgcACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAg                                                                                              <  1:1227074/150‑1 (MQ=255)
                                                         cgcACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAgg                                                                                             <  2:786733/150‑1 (MQ=255)
                                                         cgcACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAgg                                                                                             <  2:910343/150‑1 (MQ=255)
                                                          gcACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAggg                                                                                            <  1:1310831/150‑1 (MQ=255)
                                                                  ttGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAAcgc                                                                                    >  1:1326842/1‑150 (MQ=255)
                                                                  ttGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCAACTCGTCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAAcgc                                                                                    >  1:1326728/1‑150 (MQ=255)
                                                                     cTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGca                                                                                 <  2:1018842/150‑1 (MQ=255)
                                                                      tGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGCGGGTGAAGCGCGcac                                                                                >  1:1292849/1‑150 (MQ=255)
                                                                      tGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGcac                                                                                <  2:292681/150‑1 (MQ=255)
                                                                                          gTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGcc                                                              >  1:1343311/1‑149 (MQ=255)
                                                                                          gTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGcc                                                              <  2:1343311/149‑1 (MQ=255)
                                                                                                   ggCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAAcgc                                                                                    >  1:833102/1‑118 (MQ=255)
                                                                                                   ggCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAAcgc                                                                                    <  2:833102/118‑1 (MQ=255)
                                                                                                    gCCTCCTGCGCCACCTCGGCACGCCGGTGCCGGGCCGGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCCCACTCACTTCCGAGATACGCCCGGCACCCTTg                                                   <  1:534470/150‑1 (MQ=255)
                                                                                                      ctcctGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGt                                                 >  2:1263720/1‑150 (MQ=255)
                                                                                                      ctcctGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCAAGGCATGTCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGt                                                 <  2:410023/150‑1 (MQ=255)
                                                                                                       tcctGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTc                                                >  2:110071/1‑150 (MQ=255)
                                                                                                             gcCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTcaccacc                                          >  1:91645/1‑150 (MQ=255)
                                                                                                                   tCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAg                                    >  2:1228151/1‑150 (MQ=255)
                                                                                                                   tCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAg                                    >  1:31245/1‑150 (MQ=255)
                                                                                                                       cACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTtga                                <  1:110071/150‑1 (MQ=255)
                                                                                                                          ggcggTGCTTGGCCTGGACCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTtgatga                             >  2:59456/1‑150 (MQ=255)
                                                                                                                                  ttGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTTGATGATCATGAgg                     >  1:196822/1‑150 (MQ=255)
                                                                                                                                      tgccggcctgaTCGACCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTTGATGATCATGAGGATgg                 >  2:169098/5‑150 (MQ=255)
                                                                                                                                       gccggcctgaTCGACCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAg                                    <  2:291738/127‑1 (MQ=255)
                                                                                                                                       gccggcctgaTCGACCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAg                                    >  1:291738/4‑130 (MQ=255)
                                                                                                                                       ctgggcctggTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTTGATGATCATGAGGATGGc                >  2:811293/1‑150 (MQ=255)
                                                                                                                                         cggcctgaTCGACCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCga                                                                      >  2:961559/2‑94 (MQ=255)
                                                                                                                                         cggcctgaTCGACCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCga                                                                      <  1:961559/93‑1 (MQ=255)
                                                                                                                                          ggcctggTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTTGATGATCATGAGGATGGCGaa             <  2:31245/150‑1 (MQ=255)
                                                                                                                                           gcctggTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCCCCCTTGGTCACCACCACGAAGTTGATGATCATGAGGATGGCGAac            >  2:472428/1‑150 (MQ=255)
                                                                                                                                              tggTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTTGATGATCATGAGGATGGCGAAcacc         >  2:143635/1‑150 (MQ=255)
                                                                                                                                                gTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGa                                                                                               <  2:147249/62‑1 (MQ=255)
                                                                                                                                                gTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGa                                                                                               >  1:147249/1‑62 (MQ=255)
                                                                                                                                                     tCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTTGATGATCATGAGGAt                   <  2:483355/133‑1 (MQ=255)
                                                                                                                                                     tCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTTGATGATCATGAGGAt                   >  1:483355/1‑133 (MQ=255)
                                                                                                                                                     tCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTTGATGATCATGAGGATGGCGAACACCACGGCAc  <  2:1326728/150‑1 (MQ=255)
                                                                                                                                                     tCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTTGATGATCATGAGGATGGCGAACACCACGGCAc  <  2:1326842/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCATGCCACCGATCAGGTTGATGAACAGGATCAGCAAGCCGGCGATGGCGTCGCCGCGCACGAACTTGCTGGCACCGTCCATCGAGCCGTAGAACTCGGCCTCCTGCGCCACCTCGGCACGGCGGTGCTTGGCCTGGGCCTGGTCGATCAGGCCGGCATTGAGGTCGGCGTCGATGGCCATCTGCTTGCCAGGCATGGCGTCGAGGGTGAAACGCGCACTCACTTCCGAGATACGCCCGGCACCCTTGGTCACCACCACGAAGTTGATGATCATGAGGATGGCGAACACCACGGCAC  >  NC_002947/4935961‑4936258

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