Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,633,545 T→G 15.1% V165G (GTC→GGC)  lepB → signal peptidase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,633,5450TG15.1% 135.0 / 15.6 53V165G (GTC→GGC) lepBsignal peptidase I
Reads supporting (aligned to +/- strand):  ref base T (26/19);  new base G (4/4);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

AAGAAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCA  >  NC_002947/1633400‑1633689
                                                                                                                                                 |                                                                                                                                                
aagaagGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCaa                                                                                                                                               >  3:58649/1‑149 (MQ=255)
         aTCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACgggctggc                                                                                                                                      >  7:115197/1‑143 (MQ=255)
             aGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGACTGTTCGTCAACGGCCAGCCCGtt                                                                                                                                  >  6:45291/1‑149 (MQ=255)
               gTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGc                                                                                                                                <  4:58649/149‑1 (MQ=255)
                  ggTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACg                                                                                                                              <  2:212552/148‑1 (MQ=255)
                   gTGACCCGCCCCGCGGGGATGTAAGGGGGTTCCGCTACCCGCGTGACCCGAATGTCAACTACATCACGCGGGTGGGTGGCCGGCCGGGCGACGTGGTCCGCTACACCCACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaa                                                                                                                            <  8:115197/149‑1 (MQ=255)
                        ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACt                                                                                                                       >  8:236245/1‑149 (MQ=255)
                        ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGTCGTTGACGa                                                                                                                             >  1:163447/1‑143 (MQ=255)
                        ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGTCGTTGACGa                                                                                                                             <  2:163447/143‑1 (MQ=255)
                        ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGa                                                                                                                             >  4:237114/1‑143 (MQ=255)
                        ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGa                                                                                                                             <  3:237114/143‑1 (MQ=255)
                            aaCGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGTCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggtg                                                                                                                   <  5:45291/149‑1 (MQ=255)
                                         aaTGGTGTTACGCTACACGAGTGACACGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCACCGGCCAGCCCTTTGCCTACAACCTGGTGGTCAGCGAGCCgg                                                                                                      >  7:240983/1‑149 (MQ=255)
                                          aTGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCggg                                                                                                     >  4:232635/1‑149 (MQ=255)
                                                ttCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggtg                                                                                                                   <  2:7684/129‑1 (MQ=255)
                                                ttCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggtg                                                                                                                   >  1:7684/1‑129 (MQ=255)
                                                   cGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTggg                                                                                            <  7:236245/149‑1 (MQ=255)
                                                   cGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTgg                                                                                             <  3:1357/148‑1 (MQ=255)
                                                     cTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCa                                                                                          <  3:241348/149‑1 (MQ=255)
                                                              tGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGc                                                                                 <  3:232635/149‑1 (MQ=255)
                                                                 cccGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACTAGCGGCTGTTCGTCAACGGCCAGCCCGTAGCCGAACAACTGGTGGGCAGCGAGCCGGGTCCCCTGGGCAGTGCTCAGCTg                                                                               >  8:232895/1‑148 (MQ=255)
                                                                  ccGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGc                                                                                                                                >  6:36935/1‑98 (MQ=255)
                                                                  ccGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGc                                                                                                                                <  5:36935/98‑1 (MQ=255)
                                                                  ccGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCCACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCAGAACAACTGGTGGGCAGCGAGACGGGTACACTGGGCAGTGCTAAGCTGt                                                                              >  1:178230/1‑148 (MQ=255)
                                                                  ccGAATGTCAACGACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCATCGAGCCGGGTACCGAGGGCAGGGCTCAGCGGTa                                                                             >  8:144458/1‑149 (MQ=255)
                                                                     aaTGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTAc                                                                            >  8:31971/1‑147 (MQ=255)
                                                                      aTGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAg                                                                         >  2:36893/1‑149 (MQ=255)
                                                                                  tCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGgc                                                             >  1:175879/1‑149 (MQ=255)
                                                                                  tCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGgc                                                             >  1:259898/1‑149 (MQ=255)
                                                                                             gTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAaca                                                  >  2:118460/1‑149 (MQ=255)
                                                                                             gTTGGCCTGCCGGGCGACGTGGTCAGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAaca                                                  >  8:17773/1‑149 (MQ=255)
                                                                                                 gCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCt                                               <  1:118460/148‑1 (MQ=255)
                                                                                                 gCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCt                                               <  7:31971/148‑1 (MQ=255)
                                                                                                  ccTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGc                                                                                      <  3:113373/108‑1 (MQ=255)
                                                                                                  ccTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGc                                                                                      <  3:109296/108‑1 (MQ=255)
                                                                                                  ccTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGc                                                                                      >  4:113373/1‑108 (MQ=255)
                                                                                                  ccTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGc                                                                                      >  4:109296/1‑108 (MQ=255)
                                                                                                    tGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATc                                           >  4:235613/1‑149 (MQ=255)
                                                                                                          gCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAg                                     <  3:235613/149‑1 (MQ=255)
                                                                                                          gCGACGGGGTCAGCTACACCACCGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGTTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCTAACACCTGATCCGCAAg                                     <  8:240983/149‑1 (MQ=255)
                                                                                                                   tCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGc                            >  7:144646/1‑149 (MQ=255)
                                                                                                                     cGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCg                          <  8:144646/149‑1 (MQ=255)
                                                                                                                          cacCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCTAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATc                     <  1:223209/149‑1 (MQ=255)
                                                                                                                             cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCa                  >  4:58132/1‑149 (MQ=255)
                                                                                                                             cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCa                  >  2:108538/1‑149 (MQ=255)
                                                                                                                              aaCGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCaa                                      <  3:226413/128‑1 (MQ=255)
                                                                                                                              aaCGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCaa                                      >  4:226413/1‑128 (MQ=255)
                                                                                                                              aaCGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCAt                 >  6:107911/1‑149 (MQ=255)
                                                                                                                                  aCAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATgccg             >  3:155598/1‑149 (MQ=255)
                                                                                                                                        ggCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGAc       <  2:259898/149‑1 (MQ=255)
                                                                                                                                            gTTCGGCAACGGGCTGGCCGTTGACGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCAGCCGGACcag    >  5:145789/1‑148 (MQ=255)
                                                                                                                                             ttCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagca  <  4:155598/149‑1 (MQ=255)
                                                                                                                                             ttCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGAGCCGCAAGGAAATGAGCCGGTAGCGCCTGCCGCCGGCCcagcc  >  8:200958/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AAGAAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCA  >  NC_002947/1633400‑1633689

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