Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,054,127 T→G 11.4% S167R (AGC→CGC)  PP_1829 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,054,1270TG11.4% 128.7 / 4.7 50S167R (AGC→CGC) PP_1829alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (24/20);  new base G (3/3);  total (27/23)
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

GCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGC  >  NC_002947/2053990‑2054271
                                                                                                                                         |                                                                                                                                                
gcgcAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAggg                                                                                                                                       <  6:37164/149‑1 (MQ=255)
                acaTCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGgc                                                                                                                       >  2:17657/1‑149 (MQ=255)
                acaTCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCAAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGgc                                                                                                                       >  4:29965/1‑149 (MQ=255)
                        gCTGGTCCGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGtgat                                                                                                                                   <  3:74578/129‑1 (MQ=255)
                        gCTGGTCCGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGtgat                                                                                                                                   >  4:74578/1‑129 (MQ=255)
                                  ccGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCggcggc                                                                                                     >  8:235039/1‑149 (MQ=255)
                                      ccTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAgg                                                                                                 >  7:134399/1‑149 (MQ=255)
                                       cTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGc                                                                                                <  8:134399/149‑1 (MQ=255)
                                           aGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTgg                                                                                             >  3:261347/1‑148 (MQ=255)
                                             gCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGCCTGCGCTGGCGCCAGCCTGCTCGGCGGCAAGGCTGGGaa                                                                                          >  4:193444/1‑149 (MQ=255)
                                             gCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCaa                                                                                          <  7:235039/149‑1 (MQ=255)
                                             gCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCaa                                                                                          >  8:231953/1‑149 (MQ=255)
                                              cGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAAt                                                                                         <  3:29965/149‑1 (MQ=255)
                                                   gAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCGGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTg                                                                                    <  3:161023/149‑1 (MQ=255)
                                                    aGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTgc                                                                                   >  1:41210/1‑149 (MQ=255)
                                                         ggCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCtt                                                                              >  7:95624/1‑149 (MQ=255)
                                                           cACTGCCACCGCCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCg                                                                            >  7:184080/1‑149 (MQ=255)
                                                               gCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTCCAGCACCTCGGCCAGCGGCCGGCCCCTGCCGGCCAGGCTGATGACTGCGCAGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTCCGCCTTCGCtgt                                                                        >  1:133558/1‑149 (MQ=255)
                                                                ccACCGTCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCtgt                                                                        <  8:184080/148‑1 (MQ=255)
                                                                ccACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCGGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCtgt                                                                        <  5:182217/148‑1 (MQ=255)
                                                                  aCCGGCCTGTTGCGCGGGGGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGACAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGGGGCCAGGCTGGCAATCAGTGCGCCTTCGCtgtgtc                                                                     <  1:311545/149‑3 (MQ=255)
                                                                         aGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATGAg                                                              <  8:95624/149‑1 (MQ=255)
                                                                                 gggcggCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCg                                                      >  2:309901/1‑149 (MQ=255)
                                                                                     ggcAGGCGCTGGGCCAGTTGATCGAGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCAGCTCCGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGcc                                                  >  1:53003/1‑149 (MQ=255)
                                                                                          gcgcTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAgc                                             <  1:17657/149‑1 (MQ=255)
                                                                                           cgcTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAgcg                                            >  5:263756/1‑149 (MQ=255)
                                                                                           cgcTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAgcg                                            <  2:118796/149‑1 (MQ=255)
                                                                                              tGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCtgt                                                                        >  7:104514/1‑118 (MQ=255)
                                                                                              tGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCCGCGGGCCGAAGCGCgg                                         >  6:48026/1‑149 (MQ=255)
                                                                                                   cAGTTGTTCGTGCAGCACCTCGGCCAGCGGGCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCgg                                    >  1:29559/1‑149 (MQ=255)
                                                                                                   cAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTgg                                                                                             >  8:229453/1‑92 (MQ=255)
                                                                                                   cAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTgg                                                                                             <  7:229453/92‑1 (MQ=255)
                                                                                                    aGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGc                                   <  7:287640/149‑1 (MQ=255)
                                                                                                         ttCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGa                              <  7:231953/149‑1 (MQ=255)
                                                                                                           cGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAg                             <  1:309901/148‑1 (MQ=255)
                                                                                                            gTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGc                            <  1:133033/148‑1 (MQ=255)
                                                                                                            gTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCt                           >  7:54305/1‑149 (MQ=255)
                                                                                                              gcagcaCCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGc                                                   >  7:5319/1‑123 (MQ=255)
                                                                                                              gcagcaCCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGc                                                   <  8:5319/123‑1 (MQ=255)
                                                                                                               cagcaCCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGc                        <  2:41210/149‑1 (MQ=255)
                                                                                                                  gagcTCGGCCTGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGGCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGATGCGCGGGTCGGCTTTGAGCTTGCGGc                     <  2:133558/146‑1 (MQ=255)
                                                                                                                     cTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCa                                                                                       >  1:328245/1‑80 (MQ=255)
                                                                                                                     cTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTc                   >  3:85327/1‑148 (MQ=255)
                                                                                                                        ggCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAgg               <  8:54305/149‑1 (MQ=255)
                                                                                                                          accctggccggcaGCGGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGc              >  8:65462/6‑148 (MQ=255)
                                                                                                                           cAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGa            >  2:207148/1‑149 (MQ=255)
                                                                                                                            cctggccggcaGCGGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGAc           >  6:267264/4‑149 (MQ=255)
                                                                                                                            cctggccggcaGCGGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGAc           >  5:142852/4‑149 (MQ=255)
                                                                                                                                  cggcagCGGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATc     <  5:267264/149‑1 (MQ=255)
                                                                                                                                     ccgCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGc  <  3:193444/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                
GCGCAGCGGCGCTGGCACATCCAGGCTGGTCTGCCCGGCCTGCAGGCGGTCGAGCAGGGCACTGCCACCGGCCAGTTGCGCGGGCGGCAGGCGCTGGGCCAGTTGTTCGTGCAGCACCTCGGCCAGCGGCCGGCCGCTGCCGGCCAGGGTGATGACTGCGCTGGCGCCAGCCTGCTCGGCGGCCAGGCTGGCAATCAGTGCGCCTTCGCTGTGGCCAATCAGGATCAGCGGGCCGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGC  >  NC_002947/2053990‑2054271

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