Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,929,327 C→A 8.6% P239P (CCG→CCT)  PP_4334 ← ParA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,929,3270CA 8.6% 189.6 / 2.6 55P239P (CCG→CCT) PP_4334ParA family protein
Reads supporting (aligned to +/- strand):  ref base C (27/23);  new base A (3/2);  total (30/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 = 1.00e+00

GCACTGCGTTCTTGTAGGTGAGCAGGTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACGATCTGGTACGGCAGCGCCTGCTTGCGCGAACGG  >  NC_002947/4929200‑4929459
                                                                                                                               |                                                                                                                                    
gCACTGCGTTCTTGTAGGTGAGCAGGTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTa                                                                                                                              <  1:168570/136‑1 (MQ=255)
gCACTGCGTTCTTGTAGGTGAGCAGGTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTa                                                                                                                              <  1:17307/136‑1 (MQ=255)
gCACTGCGTTCTTGTAGGTGAGCAGGTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTa                                                                                                                              <  2:42070/136‑1 (MQ=255)
gCACTGCGTTCTTGTAGGTGAGCAGGTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTa                                                                                                                              <  2:630195/136‑1 (MQ=255)
      cGTTCTTGTAGGTGAGCAGGTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTg                                                                                                                        >  1:402180/1‑136 (MQ=255)
      cGTTCTTGTAGGTGAGCAGGTGCTTGAGCAGCGCCATGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTg                                                                                                                        <  1:261829/136‑1 (MQ=255)
        ttCTTGTAGGTGAGCAGGTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGcc                                                                                                                      >  2:71019/1‑136 (MQ=255)
         tCTTGTAGGTGAGCAGGTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCa                                                                                                                     >  1:329766/1‑136 (MQ=255)
                    cgCAGGTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGcc                                                                                                          <  1:260787/135‑1 (MQ=255)
                    aGCAGGTGCTTGAGCAGCGCGCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGcc                                                                                                          <  2:627090/136‑1 (MQ=255)
                       aGGTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGt                                                                                                       >  1:90807/1‑136 (MQ=255)
                        ggTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCCCGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTa                                                                                                      >  2:55572/1‑136 (MQ=255)
                        ggTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTa                                                                                                      >  2:559416/1‑136 (MQ=255)
                           gCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGgt                                                                                                   >  2:604733/1‑136 (MQ=255)
                           gCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACAGGAATGTAGCCCTGCCATACATGCTGCCCGTAGgt                                                                                                   >  2:59052/1‑136 (MQ=255)
                           gCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACAGGAATGTAGCCCTGCCATACATGCTGCCCGTAGgt                                                                                                   >  2:55792/1‑136 (MQ=255)
                             ttGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGgtgt                                                                                                 >  1:427828/1‑136 (MQ=255)
                             ttGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGgtgt                                                                                                 >  1:183710/1‑136 (MQ=255)
                                 gcagcGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACg                                                                                             <  2:311361/136‑1 (MQ=255)
                                 gcagcGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACg                                                                                             <  1:164959/136‑1 (MQ=255)
                                       cccTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGcagcag                                                                                       >  1:619221/1‑136 (MQ=255)
                                        ccTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTa                                                                                                                              >  2:662871/1‑96 (MQ=255)
                                        ccTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTa                                                                                                                              <  1:662871/96‑1 (MQ=255)
                                        ccTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCagcagc                                                                                      >  2:445047/1‑136 (MQ=255)
                                        ccTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCagcagc                                                                                      >  2:392738/1‑136 (MQ=255)
                                                     aaCCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCa                                                                         >  2:102980/1‑136 (MQ=255)
                                                     aaCCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCa                                                                         >  2:472049/1‑136 (MQ=255)
                                                     aaCCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCa                                                                         >  1:273384/1‑136 (MQ=255)
                                                     aaCCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGACCTGCCATACATGCAGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCa                                                                         >  2:172715/1‑136 (MQ=255)
                                                      aCCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAg                                                                        <  1:392738/136‑1 (MQ=255)
                                                            ggCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGc                                                                  >  2:109570/1‑136 (MQ=255)
                                                            ggCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGc                                                                  >  2:174672/1‑136 (MQ=255)
                                                                 ttGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTggg                                                             <  2:596922/136‑1 (MQ=255)
                                                                   gCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTa                                                           <  1:579166/136‑1 (MQ=255)
                                                                   gCTGTCGAACTGCGAAGGCGTGACGCCCTCGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTa                                                           <  2:374195/136‑1 (MQ=255)
                                                                    cTGTCGAACTGCGAAGGCGTGACGCCCTTTCGACTGGCATCGCGCAGGCGCGTATCCACAGGAATGTATCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTaa                                                          >  1:125593/1‑135 (MQ=255)
                                                                         gAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACgg                                                     >  2:264356/1‑136 (MQ=255)
                                                                            cTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCg                                                  <  2:303139/136‑1 (MQ=255)
                                                                                 aaGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAg                                             <  2:329766/136‑1 (MQ=255)
                                                                                  aGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAgg                                            >  1:189056/1‑136 (MQ=255)
                                                                                      gTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCATGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCg                                        <  2:22321/136‑1 (MQ=255)
                                                                                         aCGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTa                                     <  1:264356/136‑1 (MQ=255)
                                                                                           gCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACg                                   <  1:604733/136‑1 (MQ=255)
                                                                                           gCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACAGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACg                                   <  1:55792/136‑1 (MQ=255)
                                                                                           gCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACAGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACg                                   <  1:59052/136‑1 (MQ=255)
                                                                                                 gCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACGATCTgg                             >  2:7600/1‑136 (MQ=255)
                                                                                                   gACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACGATCTGGTa                           >  1:33506/1‑136 (MQ=255)
                                                                                                          aTCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACGATCTGGTACGGCAgc                    >  2:618739/1‑136 (MQ=255)
                                                                                                            cgcgCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACGATCTGGTACGGCAgcgc                  >  1:501554/1‑136 (MQ=255)
                                                                                                            cgcgCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACGATCTGGTACGGCAgcgc                  <  1:559416/136‑1 (MQ=255)
                                                                                                              cgcAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACGATCTGGTACGGCAGCGCCt                <  1:472049/136‑1 (MQ=255)
                                                                                                                  ggCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTAGGTACGATCTGGTACGGCAGCGCCTGCtt            >  2:309601/1‑136 (MQ=255)
                                                                                                                   gcgcGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACGATCTGGTACGGCAGCGCCTGCTTg           <  2:501554/136‑1 (MQ=255)
                                                                                                                    cgcgTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACGATCTGGTACGGCAGCGCCTGCTTgc          <  1:459261/136‑1 (MQ=255)
                                                                                                                            cACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACGATCTGGTACGGCAGCGCCTGCTTGCGCGAACgg  >  1:395021/1‑136 (MQ=255)
                                                                                                                               |                                                                                                                                    
GCACTGCGTTCTTGTAGGTGAGCAGGTGCTTGAGCAGCGCCCTGTAGGCGATCAACCCACGGCTCTTGCTGTCGAACTGCGAAGGCGTGACGCCCTTGCGACTGGCATCGCGCAGGCGCGTATCCACCGGAATGTAGCCCTGCCATACATGCTGCCCGTAGGTGTCACGCAGCAGCTTGAGGGTGCCCAGCGACGCCTGGGTACGACGGTCGAACAGGGTCGGTACGATCTGGTACGGCAGCGCCTGCTTGCGCGAACGG  >  NC_002947/4929200‑4929459

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