Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,693,727 C→G 7.1% P52A (CCT→GCT) ‡ PP_1489 → CheW domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,693,7270CG 7.1% 161.5 / 2.6 52P52A (CCT→GCT) ‡PP_1489CheW domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (22/26);  new base G (2/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACTTGCAACTGCGCACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCA  >  NC_002947/1693579‑1693866
                                                                                                                                                    |                                                                                                                                           
cGACTTGCAACTGCGCACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCAcc                                                                                                                                            >  1:483062/1‑150 (MQ=255)
        aaCTGCGCACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTggg                                                                                                                                    >  2:224506/1‑150 (MQ=255)
          cTGCGCACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTc                                                                                                                                  >  1:652243/1‑150 (MQ=255)
                aCCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGc                                                                                                                            >  1:528352/1‑150 (MQ=255)
                  cGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCgg                                                                                                                              >  2:43846/1‑146 (MQ=255)
                  cGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCgg                                                                                                                              <  1:43846/146‑1 (MQ=255)
                              gAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCg                                                                                                              >  1:775127/1‑150 (MQ=255)
                              gAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGAGTTTTGCCCTGAGCGTGCACGAGGTGATCGAGGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGAGGGCATTCTCGCCCATCg                                                                                                              >  1:100602/1‑150 (MQ=255)
                                acacCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTg                                                                                                            <  1:224506/150‑1 (MQ=255)
                                acacCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTg                                                                                                            <  2:775127/150‑1 (MQ=255)
                                  acCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGAGCGCCCAGTCACCTGCCGGGGTCGCGGGCATTCTCGCCCATCGTGGc                                                                                                          >  2:801374/1‑150 (MQ=255)
                                    cGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCa                                                                                                        >  2:1015929/1‑150 (MQ=255)
                                                   ccTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCAt                                                                                         >  2:611365/1‑150 (MQ=255)
                                                    cTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATc                                                                                        <  1:1015929/150‑1 (MQ=255)
                                                        aGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGAcc                                                                                    <  1:611365/150‑1 (MQ=255)
                                                          tttCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCAGGTGCCTggg                                                                                                                                    <  2:876716/100‑1 (MQ=255)
                                                          tttCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCAGGTGCCTggg                                                                                                                                    >  1:876716/1‑100 (MQ=255)
                                                                 tCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGc                                                                                                                    <  1:2900/109‑1 (MQ=255)
                                                                 tCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGc                                                                                                                    >  2:2900/1‑109 (MQ=255)
                                                                    gagCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGa                                                                        <  2:9128/150‑1 (MQ=255)
                                                                    gagCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGa                                                                        <  1:862783/150‑1 (MQ=255)
                                                                     agCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAg                                                                       <  1:730620/150‑1 (MQ=255)
                                                                          aGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTg                                                                  <  2:608797/150‑1 (MQ=255)
                                                                          aGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTg                                                                  >  2:939372/1‑150 (MQ=255)
                                                                                   cccTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCtt                                                                          >  1:58514/1‑133 (MQ=255)
                                                                                   cccTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCtt                                                                          <  2:58514/133‑1 (MQ=255)
                                                                                       gAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCCTAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGt                                                     <  1:801374/150‑1 (MQ=255)
                                                                                         gCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTgg                                                                                                                                     <  1:111871/68‑1 (MQ=255)
                                                                                         gCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTgg                                                                                                                                     >  2:111871/1‑68 (MQ=255)
                                                                                           gTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCaa                                                 <  1:939372/150‑1 (MQ=255)
                                                                                            tGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAAc                                                >  2:270929/1‑150 (MQ=255)
                                                                                             gCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACg                                               >  2:971369/1‑150 (MQ=255)
                                                                                                  aGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCa                                          <  1:1194107/150‑1 (MQ=255)
                                                                                                        tCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCAccc                                    >  2:1190524/1‑150 (MQ=255)
                                                                                                         cGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGGCTTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCttt                                                                   <  2:709015/118‑1 (MQ=255)
                                                                                                         cGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCttt                                                                   >  1:709015/1‑118 (MQ=255)
                                                                                                          gAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAg                                                                       >  1:692976/1‑113 (MQ=255)
                                                                                                          gAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAg                                                                       <  2:692976/113‑1 (MQ=255)
                                                                                                                     gcgTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctggt                       <  1:1033408/150‑1 (MQ=255)
                                                                                                                      cgTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCAGGcgggtgctggtg                      <  1:971369/150‑1 (MQ=255)
                                                                                                                       gTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCtggtgctggtgc                     <  2:528352/150‑1 (MQ=255)
                                                                                                                        tCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTgg                                                                 >  1:640400/1‑105 (MQ=255)
                                                                                                                        tCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTgg                                                                 <  2:640400/105‑1 (MQ=255)
                                                                                                                        tCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCa                    <  1:905938/150‑1 (MQ=255)
                                                                                                                        tCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCa                    <  2:652243/150‑1 (MQ=255)
                                                                                                                           gCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTa                 >  1:243980/1‑150 (MQ=255)
                                                                                                                              gCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCAGGCgggtgctggtgc                     <  1:181093/143‑1 (MQ=255)
                                                                                                                              gCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGATCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCAGGCgggtgctggtgc                     >  2:181093/1‑143 (MQ=255)
                                                                                                                                          cGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCa  <  2:483062/150‑1 (MQ=255)
                                                                                                                                            cccAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctggt                       <  1:427962/127‑1 (MQ=255)
                                                                                                                                            cccAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctggt                       >  2:427962/1‑127 (MQ=255)
                                                                                                                                             ccAGGCAGGTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCtggtgctggtgcgtt                  <  1:583386/131‑4 (MQ=255)
                                                                                                                                             ccAGGCAGGTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCtggtgctggtgcgtt                  >  2:583386/1‑128 (MQ=255)
                                                                                                                                                    cctgcctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCAGGcgggtgctggt                       >  1:148308/1‑119 (MQ=255)
                                                                                                                                                    cctgcctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCAGGcgggtgctggt                       <  2:148308/119‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CGACTTGCAACTGCGCACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCA  >  NC_002947/1693579‑1693866

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: