Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,526,729 G→T 10.5% 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,7290GT10.5% 170.4 / 8.9 57V306L (GTG→TTG) ddlBD‑alanine‑‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base G (25/25);  new base T (3/3);  total (29/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 = 8.27e-01

ATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCTGGTGGCCGAC  >  NC_002947/1526585‑1526862
                                                                                                                                                |                                                                                                                                     
aTCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCt                                                                                                                                   <  5:88560/149‑1 (MQ=255)
   gAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTgg                                                                                                                                 >  3:165277/1‑148 (MQ=255)
    aaGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCACCAGTTGCTggaa                                                                                                                               <  6:240108/149‑3 (MQ=255)
      ggCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAAt                                                                                                                             >  7:42959/1‑149 (MQ=255)
            gggCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGc                                                                                                                       <  7:148532/149‑1 (MQ=255)
             ggCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGcc                                                                                                                      <  2:172406/149‑1 (MQ=255)
                      aCGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTg                                                                                                             >  2:9377/1‑149 (MQ=255)
                          gATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGGTCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGa                                                                                                         >  2:102190/1‑149 (MQ=255)
                           aTGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGa                                                                                                         <  3:317557/148‑1 (MQ=255)
                             gCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGAcgc                                                                                                      >  3:204/1‑149 (MQ=255)
                             gCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGAcgc                                                                                                      >  4:20350/1‑149 (MQ=255)
                              cAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCCCCCGGCCTGACAGATCCTAGCCTGGGGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGAcgcg                                                                                                     <  5:103821/149‑1 (MQ=255)
                                ggACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAg                                                                                                   <  7:121082/149‑1 (MQ=255)
                                ggACGAGCAGGGCCCGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATCGCGCCCCGTGCGGATGGCCTGGACTTCCAGCAACTGCTGCTGGCAATCCTGGCCGAAAGCGTGGCGAAGCGAg                                                                                                   >  6:332352/1‑149 (MQ=255)
                                 gACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAgg                                                                                                  >  6:289538/1‑149 (MQ=255)
                                      gCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTtaat                                                                                             <  6:10065/149‑1 (MQ=255)
                                       cAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTtaat                                                                                             >  6:182799/1‑148 (MQ=255)
                                           gCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCa                                                                                        <  4:165277/149‑1 (MQ=255)
                                                        gctCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCgat                                                                           >  1:227904/1‑149 (MQ=255)
                                                          tCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCgatga                                                                         >  3:54845/1‑149 (MQ=255)
                                                           cGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCgatga                                                                         >  1:136906/1‑148 (MQ=255)
                                                            gAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCgatgat                                                                        <  1:9377/148‑1 (MQ=255)
                                                             aaGTCCACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGCCTTCCAGCAACTGGGGCTGGCAATCCGGGCCGCAAGCGTg                                                                                                             >  5:188704/1‑110 (MQ=255)
                                                                tCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGt                                                                    >  5:132687/1‑148 (MQ=255)
                                                                  aaCACCGCACCAGGCATGACAGGTCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTca                                                                  <  1:102190/148‑1 (MQ=255)
                                                                  aaCACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTcag                                                                 >  5:49170/1‑149 (MQ=255)
                                                                     acCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTcagcag                                                              >  7:236973/1‑149 (MQ=255)
                                                                           ccAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAAccc                                                        <  4:204/149‑1 (MQ=255)
                                                                                   gACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCgg                                                >  2:48285/1‑149 (MQ=255)
                                                                                   gACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCgg                                                >  2:43030/1‑149 (MQ=255)
                                                                                     cAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGcc                                              >  1:109269/1‑149 (MQ=255)
                                                                                      aGATCATAGCCTTGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCTGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTTATAACCATGCAAGGCGCGATGATACGTCAGCAGGAACCCGTTACCGGCCg                                             <  5:332352/149‑1 (MQ=255)
                                                                                      aGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAAccc                                                        <  4:53615/138‑1 (MQ=255)
                                                                                      aGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAAccc                                                        >  3:53615/1‑138 (MQ=255)
                                                                                      aGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCg                                             <  3:20350/149‑1 (MQ=255)
                                                                                        aTCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCg                                             >  8:131477/1‑147 (MQ=255)
                                                                                         tCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTTATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAg                                          >  1:293899/1‑149 (MQ=255)
                                                                                                cTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCg                                   <  8:42959/149‑1 (MQ=255)
                                                                                                   gTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGt                                 <  2:136906/148‑1 (MQ=255)
                                                                                                    tGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGc                               <  8:236973/149‑1 (MQ=255)
                                                                                                       ccATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGAGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCcgc                            >  1:40189/1‑149 (MQ=255)
                                                                                                         aTGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgt                          <  2:227904/149‑1 (MQ=255)
                                                                                                         aTGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCAAGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgt                          >  8:169328/1‑149 (MQ=255)
                                                                                                          tGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtg                         <  5:289538/149‑1 (MQ=255)
                                                                                                            gcggcCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtgg                        <  6:49170/148‑1 (MQ=255)
                                                                                                              ggcCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGc                     <  4:59957/149‑1 (MQ=255)
                                                                                                                      cGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGc             <  8:233501/149‑1 (MQ=255)
                                                                                                                         cTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCtgg          <  3:318793/149‑1 (MQ=255)
                                                                                                                           ggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAg                                          >  7:37293/1‑115 (MQ=25)
                                                                                                                           ggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAg                                          <  8:37293/115‑1 (MQ=25)
                                                                                                                           ggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCtggt         >  1:205562/1‑148 (MQ=38)
                                                                                                                           ggCCAGGATTGCCAGCACCAGTTGCTGGAAGTCCAGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCtggt         <  2:205562/148‑1 (MQ=38)
                                                                                                                                  attgccAGCACCAGTTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGTAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCTGGTGGCcgac  >  6:242015/5‑148 (MQ=255)
                                                                                                                                      ccAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtggt                       >  6:299847/1‑123 (MQ=255)
                                                                                                                                      ccAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCgtggt                       <  5:299847/123‑1 (MQ=255)
                                                                                                                                       cAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGc                 >  8:209254/1‑128 (MQ=255)
                                                                                                                                       cAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGc                 <  7:209254/128‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                     
ATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTAGCAAGCCGGTGCCGCGTGGTGCCAGCCGGCTGGTGGCCGAC  >  NC_002947/1526585‑1526862

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: ‑