Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,876,566:1 +A 8.0% coding (1014/1212 nt) xcpS → type II secretion pathway protein XcpS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,876,5661.A 8.0% 185.7 / 10.5 50coding (1014/1212 nt)xcpStype II secretion pathway protein XcpS
Reads supporting (aligned to +/- strand):  ref base . (18/28);  new base A (2/2);  total (20/30)
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.43e-01

ACCGCCATGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCCTGCTGATCGGCCTGATCGT  >  NC_002947/3876417‑3876715
                                                                                                                                                       |                                                                                                                                                     
aCCGCCATGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTgcgc                                                                                                                                                        <  1:424450/150‑1 (MQ=255)
   gCCATGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCa                                                                                                                                                    <  2:580088/150‑1 (MQ=255)
      aTGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGATCGGCAGACGCTATCCACACTGCGCGCGTGGGTGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGc                                                                                                                                                 >  1:878694/1‑150 (MQ=255)
      aTGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGc                                                                                                                                                 <  1:1317782/150‑1 (MQ=255)
       tGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCg                                                                                                                                                >  2:895157/1‑150 (MQ=255)
        gggCATGGCCCGCGGCCTGCTCGGCAAGGGAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGa                                                                                                                                               <  1:1072459/150‑1 (MQ=255)
         ggCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAg                                                                                                                                              <  1:1012887/150‑1 (MQ=255)
          gCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGc                                                                                                                                             >  2:1047858/1‑150 (MQ=255)
          gCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGc                                                                                                                                             >  1:121100/1‑150 (MQ=255)
          gCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCGCGCAGCAGGCgcag                                                                                                                                             <  1:1062417/150‑4 (MQ=255)
           cATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCa                                                                                                                                            <  2:124419/150‑1 (MQ=255)
            aTGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAg                                                                                                                                           >  1:17801/1‑150 (MQ=255)
                 ccGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTgg                                                                                                                                      <  2:100433/150‑1 (MQ=255)
                            tCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCg                                                                                                                           >  2:1180559/1‑150 (MQ=255)
                              ggCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCgag                                                                                                                         <  1:740958/150‑1 (MQ=255)
                                cAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGAt                                                                                                                       <  2:1173836/150‑1 (MQ=255)
                                       gcagcGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGTGGCCGGGCATCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCt                                                                                                                                        >  2:1100565/1‑126 (MQ=255)
                                       gcagcGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGTGGCCGGGCATCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCt                                                                                                                                        <  1:1100565/126‑1 (MQ=255)
                                              aGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTgcg                                                                                                         >  1:1181042/1‑150 (MQ=255)
                                                    tGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACt                                                                                                   >  1:1392488/1‑150 (MQ=255)
                                                      gAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTc                                                                                                 <  1:821627/150‑1 (MQ=255)
                                                           ggCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGa                                                                                            >  2:319001/1‑150 (MQ=255)
                                                            gccagccaGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGAc                                                                                           <  2:1189178/150‑1 (MQ=255)
                                                            gccagccaGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCACTGCGCCGACTTCCATGAc                                                                                           >  1:9718/1‑150 (MQ=255)
                                                              cagccagcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCgaga                                                                                                                        >  1:1241740/1‑119 (MQ=255)
                                                              cagccagcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCgaga                                                                                                                        <  2:1241740/119‑1 (MQ=255)
                                                              cagccagcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCa                                                                                         <  1:796226/150‑1 (MQ=255)
                                                               agccagcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAg                                                                                        >  1:1187953/1‑150 (MQ=255)
                                                                gccagcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAgg                                                                                       >  1:658798/1‑150 (MQ=255)
                                                                            gCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCACTGCGCCGACTTCCATGACCAGGAAATCGGCCGtt                                                                           <  2:9718/150‑1 (MQ=255)
                                                                                aGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGt                                                                       <  2:768323/150‑1 (MQ=255)
                                                                                        cGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGt                                                               <  1:888426/150‑1 (MQ=255)
                                                                                          cTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTc                                                             <  2:493978/150‑1 (MQ=255)
                                                                                           tttCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCg                                                            <  2:878694/150‑1 (MQ=255)
                                                                                                     tctGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGTTTCGTCAAACTGtt                                                  <  2:1208390/150‑1 (MQ=255)
                                                                                                     tctGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGtt                                                  <  1:1304736/150‑1 (MQ=255)
                                                                                                       ctGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTc                                                 <  2:17801/150‑1 (MQ=255)
                                                                                                              ccGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCgc                                          <  1:895157/150‑1 (MQ=255)
                                                                                                               cGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATg                                                                 <  2:995565/126‑1 (MQ=255)
                                                                                                               cGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATg                                                                 >  1:995565/1‑126 (MQ=255)
                                                                                                                    aCCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCtga                                    <  1:705774/150‑1 (MQ=255)
                                                                                                                         gTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGAcc                               <  2:1181042/150‑1 (MQ=255)
                                                                                                                         gTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGGTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGAGGAcc                               >  1:513670/1‑150 (MQ=255)
                                                                                                                                    aTCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCct                    >  1:251752/1‑150 (MQ=255)
                                                                                                                                      ccccGCGCCTGCTGCGCGCAGGCGCGCAGTATGGCAACCGTTGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCctgc                  <  1:671996/150‑1 (MQ=255)
                                                                                                                                       ccTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCctgct                 <  2:251752/150‑1 (MQ=255)
                                                                                                                                        cTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAGCTGTTCGAACCGCTGCTGATGACCTTCATCGGCctgctg                >  1:1143833/1‑150 (MQ=255)
                                                                                                                                         ctcgcCTGCGCGCAGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGAcc                               <  2:847109/132‑1 (MQ=255)
                                                                                                                                         ctcgcCTGCGCGCAGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGAcc                               >  1:847109/3‑134 (MQ=255)
                                                                                                                                          tcgcCTGCGCGCAGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGTTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCCTGCTGAt              >  1:772953/2‑150 (MQ=255)
                                                                                                                                            gcCTGCTGCGCGCAGGCGAGGAGGCTGGCAACCCGGGGGAGAGGGCGATACAGTGCGGCAACTTTCGTGACCAGTATAGCGGCGGTTGCGTCGACCGGTTCGTCCAACTCTTCGACCACATCCTGTTGGGCTTCAGCGGCCTGGTGAtgg            <  2:666652/150‑3 (MQ=255)
                                                                                                                                                    cccgcAGGCGAGTATTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCCTGCTGATCGGCCTGATc    <  2:658798/148‑1 (MQ=255)
                                                                                                                                                      cgcAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCCTGCTGATCGGCCTGATCGt  >  1:1229143/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
ACCGCCATGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCCTGCTGATCGGCCTGATCGT  >  NC_002947/3876417‑3876715

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: