Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,526,729 G→T 6.9% V306L (GTG→TTG)  ddlB → D‑alanine‑‑D‑alanine ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,526,7290GT 6.9% 184.1 / 3.8 58V306L (GTG→TTG) ddlBD‑alanine‑‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base T (2/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

GAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCTGGTGGCCGACGAGCCCG  >  NC_002947/1526588‑1526869
                                                                                                                                               |                                                                                                                                            
gAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCTCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTgg                                                                                                                                        <  1:98390/149‑1 (MQ=255)
       gggggCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTg                                                                                                                                <  4:4341/149‑1 (MQ=255)
        ggggCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTgg                                                                                                                               <  1:160825/149‑1 (MQ=255)
           gCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGcca                                                                                                                            <  6:121416/149‑2 (MQ=38)
                  gACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGt                                                                                                                     <  1:136470/149‑1 (MQ=255)
                    cGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTgg                                                                                                                   >  6:35908/1‑149 (MQ=255)
                    cGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTgg                                                                                                                   <  1:197509/149‑1 (MQ=255)
                            aGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGAcgc                                                                                                             >  4:299252/1‑147 (MQ=255)
                               aCGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGt                                                                                                        <  6:105038/149‑1 (MQ=255)
                                           ggTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGc                                                                                            <  6:176346/149‑1 (MQ=255)
                                            gTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGc                                                                                            >  6:347934/1‑148 (MQ=255)
                                                 ggCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGgcg                                                                                      >  5:273621/1‑149 (MQ=255)
                                                  gctgctCGAATTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTtaataa                                                                                                  <  2:195285/136‑1 (MQ=255)
                                                  gctgctCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTtaataa                                                                                                  >  1:195285/1‑136 (MQ=255)
                                                    tgctCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTg                                                                                                                    >  4:286324/1‑116 (MQ=255)
                                                    tgctCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTg                                                                                                                    <  3:286324/116‑1 (MQ=255)
                                                    tgctCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCCTGCCGGATGGCCTGGCCTTCAACCACCTGGTGATGCCAATCATGGCCGAAACCGTGCCGAGGCGAGGTTATTAACAATGAAAGccgc                                                                                     >  1:98461/1‑148 (MQ=255)
                                                      ctCGAAGTCAACACCACACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCga                                                                                   <  5:212551/147‑1 (MQ=255)
                                                          aaGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATAc                                                                             <  5:347934/149‑1 (MQ=255)
                                                            gTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGt                                                                           >  5:196390/1‑149 (MQ=255)
                                                               aaCACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAgg                                                                                        <  8:2559/133‑1 (MQ=255)
                                                               aaCACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGAGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACAATGCAAgg                                                                                        >  7:2559/1‑133 (MQ=255)
                                                                acacCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTcagc                                                                       <  5:35908/149‑1 (MQ=255)
                                                                            gCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCg                                                              >  6:243340/1‑146 (MQ=255)
                                                                            gCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCg                                                              <  5:243340/146‑1 (MQ=255)
                                                                                  cAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGcc                                                     <  5:15248/149‑1 (MQ=255)
                                                                                  cAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGcc                                                     >  2:180974/1‑149 (MQ=255)
                                                                                     aTCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCg                                                    >  1:138897/1‑147 (MQ=255)
                                                                                          aGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCga                                                                                   <  5:102453/111‑1 (MQ=255)
                                                                                          aGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCga                                                                                   >  6:102453/1‑111 (MQ=255)
                                                                                          aGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAg                                             >  1:361080/1‑149 (MQ=255)
                                                                                          aGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAg                                             >  2:89051/1‑149 (MQ=255)
                                                                                              cTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCg                                          <  6:196390/149‑1 (MQ=255)
                                                                                                    cccATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCaa                                              <  6:90899/139‑1 (MQ=255)
                                                                                                    cccATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCaa                                              >  5:90899/1‑139 (MQ=255)
                                                                                                     ccATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCcgc                                   >  8:198923/1‑149 (MQ=255)
                                                                                                       aTGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgt                                 >  6:275606/1‑149 (MQ=255)
                                                                                                        tGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtg                                <  1:180974/149‑1 (MQ=255)
                                                                                                        tGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtg                                <  7:198923/149‑1 (MQ=255)
                                                                                                          ggggcCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtggt                              <  7:327960/147‑1 (MQ=255)
                                                                                                          gcggcCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGtt                                                                                                       <  3:45778/76‑1 (MQ=255)
                                                                                                          gcggcCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGtt                                                                                                       >  4:45778/1‑76 (MQ=255)
                                                                                                          gcggcCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtggt                              >  2:70392/1‑149 (MQ=255)
                                                                                                                       ttcggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCcgc                                   <  8:112801/129‑1 (MQ=25)
                                                                                                                       ttcggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCcgc                                   >  7:112801/4‑132 (MQ=25)
                                                                                                                       gCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCtg                  >  1:353256/1‑149 (MQ=255)
                                                                                                                       gCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCtg                  >  1:347667/1‑149 (MQ=255)
                                                                                                                        cTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCtgg                 <  2:361080/149‑1 (MQ=255)
                                                                                                                            ccTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCTGGTGGc             >  7:125488/1‑149 (MQ=255)
                                                                                                                               ggattgccAGCACCAGTTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCTGGTGGCcga          >  4:77133/7‑149 (MQ=255)
                                                                                                                                gACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCTGGTGGCcgac         >  5:332546/1‑149 (MQ=255)
                                                                                                                                  cTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCTGGTGGCcgac         >  5:225516/1‑147 (MQ=255)
                                                                                                                                      cAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCTTTACCGGCCGTAACAAGCCCGTGCCGCgtg                                >  8:273015/1‑120 (MQ=255)
                                                                                                                                      cAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCTGGTGGCCGACGAGccc   >  1:258492/1‑149 (MQ=255)
                                                                                                                                      cAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGAGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACGGGCCGTAGCAAGCCGGTGCCGCgtg                                <  7:273015/120‑1 (MQ=255)
                                                                                                                                        gCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCgatgat                                                                               <  6:95492/71‑1 (MQ=255)
                                                                                                                                        gCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCgatgat                                                                               >  5:95492/1‑71 (MQ=255)
                                                                                                                                        gCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCTGGTGGCCGACGAGCCCg  >  4:64467/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                            
GAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCTGGTGGCCGACGAGCCCG  >  NC_002947/1526588‑1526869

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