Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,477,214 C→T 7.8% intergenic (+149/‑133) PP_1290 → / → PP_1291 polysaccharide deacetylase family protein/PhoH family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,477,2140CT 7.8% 145.0 / 2.2 51intergenic (+149/‑133)PP_1290/PP_1291polysaccharide deacetylase family protein/PhoH family protein
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base T (2/2);  total (28/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.99e-01

GGTATGACTTCTGAGCGTAGTCCCCTCTGACGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGAC  >  NC_002947/1477068‑1477355
                                                                                                                                                    |                                                                                                                                             
ggTATGACTTCTGAGCGTAGTCCCCTCTGACGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTgcagca                                                                                                                                              >  1:26682/1‑148 (MQ=255)
 gTATGACTTCTGAGCGTAGTCCCCTCTGACGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGcaca                                                                                                                                            <  2:183054/149‑1 (MQ=255)
        ttCTGAGCGTAGTCCCCTCTGACGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCg                                                                                                                                     >  6:203733/1‑149 (MQ=255)
        ttCTGAGCGTAGTCCCCTCTGACGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCg                                                                                                                                     <  3:324252/149‑1 (MQ=255)
        ttCTGAGCGTAGTCCCCTCTGACGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAAATTAAAAAAGCCGTAAATTTGGTTTTTACGGACATGGTTTTGGCGGAATCAAGAAACAACACAGCCCATGCCGCAGGCACAGTTGGCg                                                                                                                                     >  8:341424/1‑149 (MQ=255)
        ttCTGAGCGTAGTCCCCTCTGACGCTGGCGCCAAGGCCTCGTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCg                                                                                                                                     >  4:74388/1‑149 (MQ=255)
             aGCGTAGTCCCCTCTGACGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATc                                                                                                                                >  5:151625/1‑149 (MQ=255)
                 tattCCCCTCTGACGCTTGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAACTGCTTTTTTCGGTCATGGTTTTTGCAGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAg                                                                                                                             <  1:94147/146‑1 (MQ=255)
                      cccTCTGACGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCAc                                                                                                                       >  3:100988/1‑149 (MQ=255)
                         tctGACGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCAcc                                                                                                                      >  1:362777/1‑147 (MQ=255)
                           tGACGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGc                                                                                                                  >  8:48843/1‑149 (MQ=255)
                             aCGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGc                                                                                                                  >  5:185640/1‑147 (MQ=255)
                              acccggCGCCAATGCCTCTTCCGTAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTAGGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAg                                                                                                               <  7:341424/145‑1 (MQ=255)
                                cTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGt                                                                                                             >  6:314161/1‑149 (MQ=255)
                                     gcCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCt                                                                                                        <  1:99725/149‑1 (MQ=255)
                                       cAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTc                                                                                                      >  5:76567/1‑149 (MQ=255)
                                             ctctTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTcccgccc                                                                                                <  6:76567/149‑1 (MQ=255)
                                                   gTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATc                                                                                                                                >  3:93696/1‑111 (MQ=255)
                                                   gTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATc                                                                                                                                <  4:93696/111‑1 (MQ=255)
                                                   gTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTcc                                                                                                     <  6:330394/138‑1 (MQ=255)
                                                   gTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTcc                                                                                                     >  5:330394/1‑138 (MQ=255)
                                                       aTCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAAcgc                                                                                      <  3:74388/149‑1 (MQ=255)
                                                       aTCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAAcgc                                                                                      <  2:362777/149‑1 (MQ=255)
                                                               aaTAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACgg                                                                              <  7:48843/149‑1 (MQ=255)
                                                               aaTAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACgg                                                                              <  5:203733/149‑1 (MQ=255)
                                                               aaTAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACgg                                                                              <  4:100988/149‑1 (MQ=255)
                                                                aTAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACggg                                                                             <  6:79306/149‑1 (MQ=255)
                                                                      ttCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGTACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGc                                                                       >  3:252101/1‑149 (MQ=255)
                                                                      ttCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGTACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGc                                                                       >  3:247264/1‑149 (MQ=255)
                                                                      ttCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGAGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCACGCCCGTAACAACGCGGATACGGGGACAGc                                                                       >  2:295291/1‑149 (MQ=255)
                                                                          aaaaTGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGTACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCa                                                                   <  4:252101/149‑1 (MQ=255)
                                                                          aaaaTGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGTACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCa                                                                   <  4:247264/149‑1 (MQ=255)
                                                                           aaaTGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAg                                                                                                                                           >  5:108807/1‑76 (MQ=255)
                                                                           aaaTGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAg                                                                                                                                           <  6:108807/76‑1 (MQ=255)
                                                                           aaaTGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAg                                                                  >  4:317318/1‑149 (MQ=255)
                                                                                  tAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCActct                                                           >  6:70091/1‑149 (MQ=255)
                                                                                                   gTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCAc                                          <  6:185640/149‑1 (MQ=255)
                                                                                                      aTGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGCTTCTCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACAACGcc                                       >  8:309566/1‑149 (MQ=255)
                                                                                                                     aTGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGCCa                         >  3:72181/1‑149 (MQ=255)
                                                                                                                      tGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACa                         >  7:204869/1‑148 (MQ=255)
                                                                                                                      tGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACa                         >  8:24643/1‑148 (MQ=255)
                                                                                                                        aagaaACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACAt                        >  8:344561/1‑147 (MQ=255)
                                                                                                                        aagaaACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATa                       <  8:204869/148‑1 (MQ=255)
                                                                                                                          gaaACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCgtgt                                   >  4:257899/1‑134 (MQ=255)
                                                                                                                          gaaACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCgtgt                                   <  3:257899/134‑1 (MQ=255)
                                                                                                                          gaaACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAgg                    <  3:317318/149‑1 (MQ=255)
                                                                                                                              ccAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGAc                <  7:24643/149‑1 (MQ=255)
                                                                                                                              ccAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGAAACGGGGACAGCGGCAGTCACTCTGCGTCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGAc                <  7:309566/149‑1 (MQ=255)
                                                                                                                                      ccGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATg        >  8:303973/1‑149 (MQ=255)
                                                                                                                                        gATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGt                   <  3:168894/136‑1 (MQ=255)
                                                                                                                                        gATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGt                   >  4:168894/1‑136 (MQ=255)
                                                                                                                                        gATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGa      >  3:334609/1‑149 (MQ=255)
                                                                                                                                             cTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGAc  >  2:319027/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GGTATGACTTCTGAGCGTAGTCCCCTCTGACGCTGGCGCCAAGGCCTCTTCGTCAATCGGAAAAATAAACTTCAAAAATGCGTAAAAGTGCTTTTTTCGGTCATGGTTTTTGCGGTATGAAGAAACCAGACAGCCGATTCCTGCAGCACAGGTGGCGTCATCCAGGCCCACCTAGCCAGGTTCGCTTCCCGCCCGTAACAACGCGGATACGGGGACAGCGGCAGTCACTCTGCGGCGTAACAGACCACGCCGTGTGGCTTCGACATAAGGTGACCGAGTATGGATGAC  >  NC_002947/1477068‑1477355

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: