Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,915,015 G→A 8.1% E66K (GAA→AAA)  fklB‑II → FKBP‑type peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,915,0150GA 8.1% 136.2 / 2.4 49E66K (GAA→AAA) fklB‑IIFKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base A (2/2);  total (24/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 = 2.34e-01

TGGCTGGCTGCGATCAACAAGCCAGCAGCCCCGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAGCCGCCTCCGCCGGCAAGAAATTCCT  >  NC_002947/1914869‑1915151
                                                                                                                                                  |                                                                                                                                        
tggctggcTGCGATCAACAAGCCAGCAGCCCCGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCg                                                                                                                                          >  7:240104/1‑147 (MQ=255)
  gctggctgCGATCAACAAGCCAGCAGCCCCGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGa                                                                                                                                      >  5:158985/1‑149 (MQ=255)
  gctggctgCGATCAACAAGCCAGCAGCCCCCAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGa                                                                                                                                      >  1:21578/1‑149 (MQ=255)
   ctggctgCGATCAACAAGCCAGCAGCCCCGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGAc                                                                                                                                     <  6:158985/149‑1 (MQ=255)
        tgCGATCAACAAGCCAGCAGCCCCGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGt                                                                                                                                >  3:228691/1‑149 (MQ=255)
          cGATCAACAAGCCAGCAGCCCCGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGt                                                                                                                                <  1:132018/147‑1 (MQ=255)
          cGATCAACAAGCCAGCAGCCCCGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGt                                                                                                                                <  6:193402/147‑1 (MQ=255)
            aTCAACAAGCCAGCAGCCCCGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAg                                                                                                                             >  3:148659/1‑148 (MQ=255)
                  aacccAGCAGCCCCGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCaagaaa                                                                                                                      <  4:171007/146‑1 (MQ=255)
                     ccAGCAGCCCCGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAagaaagaa                                                                                                                   <  6:129194/149‑1 (MQ=255)
                            ccccGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGc                                                                                                             >  6:119869/1‑148 (MQ=255)
                            ccccGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGACTGGCTCAGAAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCa                                                                                                            >  7:238274/1‑149 (MQ=255)
                              ccGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCAt                                                                                                           >  3:70865/1‑148 (MQ=255)
                                gAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCaa                                                                                                        <  5:119869/149‑1 (MQ=255)
                                gAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCaa                                                                                                        <  3:23345/149‑1 (MQ=255)
                                 aGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAg                                                                                                       <  6:236531/149‑1 (MQ=255)
                                   cTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGa                                                                                                     <  2:21578/149‑1 (MQ=255)
                                           tCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGc                                                                                             >  3:168963/1‑149 (MQ=255)
                                             cGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGc                                                                                             <  2:122888/147‑1 (MQ=255)
                                                   aGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTagag                                                                                     >  3:54207/1‑149 (MQ=255)
                                                     aaaGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGc                                                                                   >  5:171103/1‑149 (MQ=255)
                                                       aGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCt                                                                                 <  4:228691/149‑1 (MQ=255)
                                                          ctcctACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCa                                                                              <  2:37441/149‑1 (MQ=255)
                                                           tcctACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCAc                                                                             >  8:118103/1‑149 (MQ=255)
                                                                           tGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCAAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAAcgc                                                             >  3:131475/1‑149 (MQ=255)
                                                                           tGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCAAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAAcgc                                                             >  1:47583/1‑149 (MQ=255)
                                                                             aaCATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGATAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGACGCCTTCACCGCATTGCAGACAcgcga                                                           >  4:9803/1‑148 (MQ=255)
                                                                               cATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCAAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCg                                                         <  4:131475/149‑1 (MQ=255)
                                                                               cATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCAAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCg                                                         <  2:47583/149‑1 (MQ=255)
                                                                                 tGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCgag                                                       <  6:171103/149‑1 (MQ=255)
                                                                                   ggCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCgagg                                                      <  4:54207/148‑1 (MQ=255)
                                                                                       agagCCGGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAgcgc                                                 <  8:240104/149‑1 (MQ=255)
                                                                                        gagCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGcc                                                >  8:143525/1‑149 (MQ=255)
                                                                                          gCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTg                                              <  7:118103/149‑1 (MQ=255)
                                                                                            cTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTgg                                             <  8:238274/148‑1 (MQ=255)
                                                                                               gCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCaa                                         <  4:148659/149‑1 (MQ=255)
                                                                                                    ggAAGGCATGGAAGACCTTGATTCGAACGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCa                                    <  3:9803/149‑1 (MQ=255)
                                                                                                     gAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAg                                   >  3:175219/1‑149 (MQ=255)
                                                                                                           aTGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCgaag                              >  5:42134/1‑148 (MQ=255)
                                                                                                              gAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAgc                          >  5:82778/1‑149 (MQ=255)
                                                                                                               aaGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAgcc                         <  4:70865/149‑1 (MQ=255)
                                                                                                                            gAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAGCCGCCTCCGCCGGc             >  6:274982/1‑148 (MQ=255)
                                                                                                                             aaaGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAGCCGCCTCCGCCGGc             >  7:116263/1‑147 (MQ=255)
                                                                                                                               aGCAGTCGCCCTCGGCATCGAAGACGCCGATAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAGCCGCCTCCGCCGGCAAg          >  1:157922/1‑148 (MQ=255)
                                                                                                                                    tCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAGCCGCCTCCGCCGGc             <  3:188364/140‑1 (MQ=255)
                                                                                                                                    tCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAGCCGCCTCCGCCGGc             >  4:188364/1‑140 (MQ=255)
                                                                                                                                    tCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAGCCGCCTCCGCCGGCAAGAAATTc    >  3:122275/1‑149 (MQ=255)
                                                                                                                                      gCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAGCCGCCTCCGCCGGCAAGAAATTCCt  <  6:42134/149‑1 (MQ=255)
                                                                                                                                           cGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAGCCGCCTCCGCCGGCAAGAAATTc    <  3:93793/142‑1 (MQ=255)
                                                                                                                                           cGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAGCCGCCTCCGCCGGCAAGAAATTc    >  4:93793/1‑142 (MQ=255)
                                                                                                                                                  |                                                                                                                                        
TGGCTGGCTGCGATCAACAAGCCAGCAGCCCCGAGCTGAAGACTCCGGCGCAGAAAGCCTCCTACGGTATCGGCCTGAACATGGGCAAGAGCCTGGCTCAGGAAGGCATGGAAGACCTTGATTCGAAAGCAGTCGCCCTCGGCATCGAAGACGCCGTTAGCAAGAAAGAACAACGCATCAAGGATGAAGAGCTGGTAGAGGCCTTCACCGCATTGCAGAAACGCGCCGAGGAGCGCCTGGCCAAGGCCAGCGAAGAAGCCGCCTCCGCCGGCAAGAAATTCCT  >  NC_002947/1914869‑1915151

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: