Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,976,212 C→G 7.4% D187H (GAC→CAC)  flgJ ← peptidoglycan hydrolase FlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,976,2120CG 7.4% 170.3 / 2.9 54D187H (GAC→CAC) flgJpeptidoglycan hydrolase FlgJ
Reads supporting (aligned to +/- strand):  ref base C (28/22);  new base G (2/2);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GCTCGGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATGGCA  >  NC_002947/4976065‑4976357
                                                                                                                                                   |                                                                                                                                                  
gCTCGGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCa                                                                                                                                                   >  7:159758/1‑149 (MQ=255)
     gCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACggg                                                                                                                                               >  2:65937/1‑148 (MQ=255)
     gCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggc                                                                                                                                              <  5:182437/149‑1 (MQ=255)
      ccATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggca                                                                                                                                             >  8:4389/1‑149 (MQ=255)
      ccATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggca                                                                                                                                             <  3:14285/149‑1 (MQ=255)
         tttgcagcaTGGTGGCGACGAAGGCGGCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGCCCGTGACGgtcaggc                                                                                                                                          <  2:251320/146‑1 (MQ=255)
              agcaTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGCGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggcaggc                                                                                                                                          >  8:288472/1‑144 (MQ=255)
              agcaTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGCGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggcaggc                                                                                                                                          <  7:288472/144‑1 (MQ=255)
                 aTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACggcggc                                                                                                                                  >  6:341461/1‑149 (MQ=255)
                 aTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACggcggc                                                                                                                                  <  5:341461/149‑1 (MQ=255)
                 aTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACggcggc                                                                                                                                  <  2:291546/149‑1 (MQ=255)
                         gacgaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTgg                                                                                                                          >  6:100529/1‑149 (MQ=255)
                          acgaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCAGGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGc                                                                                                                         <  5:150047/149‑1 (MQ=255)
                               ggCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGttt                                                                                                                    >  5:126198/1‑149 (MQ=255)
                                gCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTc                                                                                                                   >  4:188720/1‑149 (MQ=255)
                                    cGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggc                                                                                                                                              >  3:41359/1‑118 (MQ=255)
                                    cGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggc                                                                                                                                              <  4:41359/118‑1 (MQ=255)
                                     gCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTccgcc                                                                                                               >  1:233092/1‑148 (MQ=255)
                                      cTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGa                                                                                                             <  5:100529/149‑1 (MQ=255)
                                             tGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCAc                                                                                                       >  7:309455/1‑148 (MQ=255)
                                                            gcgcCCGGGGTGTCGGGGCCTCGCCACGGCCCACAATGCGCAGCCCGTGGCCCGGCCCGGCTAGGGCCTGGGCTAGGGCCTCCCCGCCCCGGGCTGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAgg                                                                                       <  8:125028/149‑1 (MQ=255)
                                                            gcgcCAGGGGGGGCGGGGCCCGGGCAGGGCCGACAATGGGCCGGCCGTTGTCCGGCCGGGCTAAGGCCTTGGCTACTGCGTGCCCGCCCCGGGCTGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAgg                                                                                       <  3:188720/149‑1 (MQ=255)
                                                            gcgcCAGCGGTGGCTGGGCCACGGCAGGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCTAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAgg                                                                                       <  7:4389/149‑1 (MQ=255)
                                                                  gCGGGGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcggtcg                                                                                 <  8:163329/149‑1 (MQ=255)
                                                                    ggTGGCTGGGCCACGGCACGGCCTACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGTCGGCACGATACCGGCCAACAGGcggtcggt                                                                               <  8:78840/149‑1 (MQ=255)
                                                                    ggTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcggtcggt                                                                               <  8:159758/149‑1 (MQ=255)
                                                                     gTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGAc                                                                                                            <  3:126076/119‑1 (MQ=255)
                                                                     gTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGAc                                                                                                            >  4:126076/1‑119 (MQ=255)
                                                                     gTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcggtcggt                                                                               >  5:157289/1‑148 (MQ=255)
                                                                        ctggcgCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGt                                                                           <  2:245826/144‑1 (MQ=255)
                                                                         cggggCCCCGGCCCGCCCGACATGCCCCCGGCCGGTGTCGGCCACGGCGACGCCCTTGGCTACTCCCTCCCCGCCCCGGCCCGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGtt                                                                          <  8:309455/147‑1 (MQ=255)
                                                                          tGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGc                                                                                      <  8:174069/136‑1 (MQ=255)
                                                                          tGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGc                                                                                      >  7:174069/1‑136 (MQ=255)
                                                                             gCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCt                                                                      >  2:227548/1‑149 (MQ=255)
                                                                                  ggcacggcCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCATCAGGCGGTCGGTGAGTTTGCTTgg                                                                   >  8:306392/1‑147 (MQ=255)
                                                                                         ccGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGCCCGTCACGGGCAGGCACGGCGGCGTTAGTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGGTTGCTTGGCAGGGc‑ccg                                                         >  2:224905/1‑147 (MQ=255)
                                                                                                   gcAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCcgcgcg                                               >  1:187579/1‑149 (MQ=255)
                                                                                                               ccGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTgcagc                                    >  6:381153/1‑148 (MQ=255)
                                                                                                                  gCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCACTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAAcg                                >  5:300275/1‑149 (MQ=255)
                                                                                                                       gCGAAGGCCTTGGCTACTGCCTGCCCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCCACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGtt                           >  2:87226/1‑149 (MQ=255)
                                                                                                                       gCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGtt                           <  2:290632/149‑1 (MQ=255)
                                                                                                                        cGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTg                          >  7:150575/1‑149 (MQ=255)
                                                                                                                            ggCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCACAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTgcgc                       <  1:65937/149‑1 (MQ=255)
                                                                                                                              ccTTGGCTACTGCCTGACCGTCAAGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGcgcg                      >  7:235638/1‑147 (MQ=255)
                                                                                                                                   gCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGa               >  6:221148/1‑149 (MQ=255)
                                                                                                                                    cTACTGCCTTACCGGCCCGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGAc              <  6:21616/149‑1 (MQ=255)
                                                                                                                                     ccgtgcctgccCGTGACGGTCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCACCAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACg             >  5:235707/4‑149 (MQ=255)
                                                                                                                                     ccgtgcctgccCGTGACGGTCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGc                                                             <  2:215351/99‑1 (MQ=39)
                                                                                                                                     ccgtgcctgccCGTGACGGTCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGc                                                             >  1:215351/4‑102 (MQ=39)
                                                                                                                                          ccTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGc         >  2:369398/1‑148 (MQ=255)
                                                                                                                                          ccTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGc         >  2:266590/1‑148 (MQ=255)
                                                                                                                                          ccTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCAGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGc         >  6:1434/1‑148 (MQ=255)
                                                                                                                                              aCCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATg     <  6:126198/148‑1 (MQ=255)
                                                                                                                                                cGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATGGCa  >  2:361379/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTCGGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGC‑CAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATGGCA  >  NC_002947/4976065‑4976357

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