Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,822,277 G→A 7.7% A186T (GCC→ACC)  PP_2475 → TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,822,2770GA 7.7% 148.1 / 2.2 52A186T (GCC→ACC) PP_2475TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base A (2/2);  total (25/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 = 1.84e-01

TGAATGTATTGAGCAGGGGCAGCGAGACGGTTCTCTGGCAATGGGCGATGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTG  >  NC_002947/2822132‑2822424
                                                                                                                                                 |                                                                                                                                                    
tGAATGTATTGAGCAGGGGCAGCGAGACGGTTCTCTGGCAATGGGCGATGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCca                                                                                                                                                   >  7:237528/1‑149 (MQ=255)
 gAATGTATTGAGCAGGGGCAGCGAGACGGTTCTCTGGCAATGGGCGATGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGATGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCca                                                                                                                                                   >  5:8536/1‑148 (MQ=255)
    tGTATTGAGCAGGGGCAGCGAGACGGTTCTCTGGCAATGGGCGATGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAAc                                                                                                                                               >  8:26209/1‑149 (MQ=255)
        ttGAGCAGGGGCAGCGAGACGGTTCTCTGGCAATGGGCGATGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCa                                                                                                                                            <  5:91569/148‑1 (MQ=255)
          gAGCAGGGGCAGCGAGACGGTTCTCTGGCAATGGGCGATGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGTCACAACTCAGa                                                                                                                                          <  8:48710/148‑1 (MQ=255)
          gAGCAGGGGCAGCGAGACGGTTCTCTGGCAATGGGCGATGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTc                                                                                                                                             <  4:127185/145‑1 (MQ=255)
          gAGCAGGGGCAGCGAGACGGTTCTCTGGCAATGGGCGATGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTc                                                                                                                                             >  3:127185/1‑145 (MQ=255)
           aGCAGGGGCAGCGAGACGGTTCTCTGGCAATGGGCGATGCATTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGAGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGAc                                                                                                                                         >  1:205231/1‑148 (MQ=255)
              aGGTGCAGCGAGACGGTTCTCTGGCAATGGGCGATGCTTTCGACTCTGCGACTGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCAtt                                                                                                                                     <  6:8536/149‑1 (MQ=255)
                                 tctGGCAATGGGCGATGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGTTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCt                                                                                                                  >  2:98720/1‑149 (MQ=255)
                                       aaTGGGCGATGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGgc                                                                                                           >  3:178253/1‑149 (MQ=255)
                                               aTGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGc                                                                                                   >  5:217017/1‑149 (MQ=255)
                                                tGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGc                                                                                                   <  6:217017/148‑1 (MQ=255)
                                                  cTTTCCACACTGCGAGGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTTCTGAGCAAGCTTCACCGTAATAGTCGGGCGCTGGAAATCGGCATGGCCACAACTGAGACGATTCTGTCCTCTGACAACATATGATGGCGCGGTTGCCCa                                                                                                <  2:205231/149‑1 (MQ=255)
                                                  cTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCa                                                                                                >  4:141487/1‑149 (MQ=255)
                                                   tttCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCAAAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCAc                                                                                               >  8:32142/1‑149 (MQ=255)
                                                   tttCCACACTCCGACGGCCCTTTACCAGTTGTGGCACGGTTCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCAc                                                                                               <  5:43126/149‑1 (MQ=255)
                                                        acacTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGtt                                                                                          >  3:12628/1‑149 (MQ=255)
                                                         cacTGCGACGGCCCTTTACCAGTTGTGGATCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTtg                                                                                         <  4:4453/149‑1 (MQ=255)
                                                            tGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGcc                                                                                                                                                    >  6:21468/1‑88 (MQ=255)
                                                            tGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGcc                                                                                                                                                    <  5:21468/88‑1 (MQ=255)
                                                              cGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTtggtgg                                                                                     <  3:141487/148‑1 (MQ=255)
                                                              cGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTtggtgg                                                                                     <  8:237528/148‑1 (MQ=255)
                                                                       tttACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCAtttttt                                                                            >  2:41500/1‑148 (MQ=255)
                                                                         tACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTc                                                                          >  7:31534/1‑148 (MQ=255)
                                                                          aCCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCa                                                                         <  7:32142/148‑1 (MQ=255)
                                                                          aCCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCa                                                                         <  5:300577/147‑1 (MQ=255)
                                                                          aCCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCa                                                                         >  6:300577/1‑147 (MQ=255)
                                                                            cAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTtggt                                                                                       >  4:319396/1‑132 (MQ=255)
                                                                            cAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTtggt                                                                                       <  3:319396/132‑1 (MQ=255)
                                                                            cAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACt                                                                       >  6:251963/1‑148 (MQ=255)
                                                                             aGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGa                                                                     <  5:251963/149‑1 (MQ=255)
                                                                              gTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTtg                                                                                         >  4:207840/1‑128 (MQ=255)
                                                                              gTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTtg                                                                                         <  3:207840/128‑1 (MQ=255)
                                                                                  tgGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCaa                                                                 <  4:178253/148‑1 (MQ=255)
                                                                                      tCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGt                                                            >  4:302826/1‑149 (MQ=255)
                                                                                      tCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATATTTTTCACTGACCAagtagt                                                            >  4:298095/1‑149 (MQ=255)
                                                                                        ggTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCGACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGt                                                            <  3:298095/147‑1 (MQ=255)
                                                                                        ggTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGt                                                            <  3:302826/147‑1 (MQ=255)
                                                                                          tGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGACCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGAc                                                        >  2:100570/1‑149 (MQ=255)
                                                                                                      gAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACt                                                               >  3:309657/1‑130 (MQ=255)
                                                                                                      gAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACt                                                               <  4:309657/130‑1 (MQ=255)
                                                                                                        gCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGTACGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCgg                                           <  1:98720/149‑1 (MQ=255)
                                                                                                           aGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGa                                       >  8:198046/1‑149 (MQ=255)
                                                                                                                cacCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCga                                   <  3:122636/148‑1 (MQ=255)
                                                                                                                   cGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGAc                                >  2:24209/1‑148 (MQ=255)
                                                                                                                    gTAATAGTCGGCCGCTGGAAAACGCCATGACCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACaa                              <  1:100570/149‑1 (MQ=255)
                                                                                                                           tCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAAc                       <  1:24209/149‑1 (MQ=255)
                                                                                                                             ggCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACt                      <  7:198046/148‑1 (MQ=255)
                                                                                                                                     gAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAAc              <  7:26209/148‑1 (MQ=255)
                                                                                                                                            ccATGACCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAAc                       <  4:13641/132‑1 (MQ=255)
                                                                                                                                            ccATGACCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAAc                       >  3:13641/1‑132 (MQ=255)
                                                                                                                                                 gCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTATTTCGGGTGAGCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTg  >  8:282410/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGAATGTATTGAGCAGGGGCAGCGAGACGGTTCTCTGGCAATGGGCGATGCTTTCCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTG  >  NC_002947/2822132‑2822424

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: