Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,067,198 G→C 7.5% intergenic (+14/‑160) phoR → / → PP_5322 two‑component system sensor histidine kinase/metal ion transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,067,1980GC 7.5% 164.9 / 4.3 53intergenic (+14/‑160)phoR/PP_5322two‑component system sensor histidine kinase/metal ion transporter
Reads supporting (aligned to +/- strand):  ref base G (25/24);  new base C (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GGGCTGGGGCTCGCAATCGTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGAC  >  NC_002947/6067050‑6067336
                                                                                                                                                    |                                                                                                                                          
gggCTGGGGCTCGCAATCGTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACC‑GGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTgcggc                                                                                                                                           >  4:157473/1‑149 (MQ=255)
            gCAATCGTCAAGCATGTGCTGATGCTCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGCCGCAGcctcaa                                                                                                                                 <  4:26361/148‑5 (MQ=255)
              aaTCGTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACAGGGGCATGGCAGCACCTGTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCGTGAGATTGAGGCTGCGGCGCAGCCCTTcgcg                                                                                                                              >  6:258534/1‑149 (MQ=255)
                  gTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCa                                                                                                                          >  2:167861/1‑149 (MQ=255)
                      aGCATGTGCTGATGAGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGcc                                                                                                                      >  6:8668/1‑149 (MQ=255)
                            tgCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGc                                                                                                                <  5:281193/149‑1 (MQ=255)
                              cTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCa                                                                                                              <  3:182247/149‑1 (MQ=255)
                                    cgcCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTggg                                                                                                         <  8:142029/148‑1 (MQ=255)
                                    cgcCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGc                                                                                                        <  3:191258/149‑1 (MQ=255)
                                       cACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCg                                                                                                     >  2:220849/1‑149 (MQ=255)
                                                      gAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCt                                                                                        <  5:8668/147‑1 (MQ=255)
                                                        aaTCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGccc                                                                                    >  4:226356/1‑149 (MQ=255)
                                                                 tgtACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCttt                                                                            <  6:320550/148‑1 (MQ=255)
                                                                 tgtACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCttt                                                                            >  5:320550/1‑148 (MQ=255)
                                                                            aTGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGcc                                                                     <  3:137876/144‑1 (MQ=255)
                                                                            aTGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGcc                                                                     <  3:143160/144‑1 (MQ=255)
                                                                            aTGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGcc                                                                     >  4:137876/1‑144 (MQ=255)
                                                                            aTGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGcc                                                                     >  4:143160/1‑144 (MQ=255)
                                                                                cagcaCCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAg                                                                                              >  3:157696/1‑115 (MQ=255)
                                                                                cagcaCCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAg                                                                                              <  4:157696/115‑1 (MQ=255)
                                                                                 agcaCCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCtt                                                                                  <  2:350596/126‑1 (MQ=255)
                                                                                 agcaCCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCtt                                                                                  >  1:350596/1‑126 (MQ=255)
                                                                                      cTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTg                                                      >  6:241756/1‑149 (MQ=255)
                                                                                             tGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATcccc                                               >  3:283550/1‑149 (MQ=255)
                                                                                              gTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTgg                                                     >  7:163683/1‑142 (MQ=255)
                                                                                              gTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTgg                                                     <  8:163683/142‑1 (MQ=255)
                                                                                              gTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCt                                              >  2:260597/1‑149 (MQ=255)
                                                                                               tCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCg                                                                                                     >  5:327976/1‑93 (MQ=255)
                                                                                               tCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCg                                                                                                     <  6:327976/93‑1 (MQ=255)
                                                                                               tCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCGGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCt                                              <  8:97763/148‑1 (MQ=255)
                                                                                                 aCTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCa                                                                  >  4:6473/1‑126 (MQ=255)
                                                                                                 aCTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCa                                                                  <  3:6473/126‑1 (MQ=255)
                                                                                                           aGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAg                                 <  5:241756/149‑1 (MQ=255)
                                                                                                              acaATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTg                              <  4:283550/149‑1 (MQ=255)
                                                                                                              acaATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTg                              <  3:226356/149‑1 (MQ=255)
                                                                                                                    ggCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCTTTTGCTGGCAca                         <  1:167861/148‑2 (MQ=255)
                                                                                                                     gCCAACGAACAGCCGAGAGATTGAGGCACCCGCGCAGCCCTTCGCCGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACttt                       <  4:199635/149‑1 (MQ=255)
                                                                                                                       cAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTg                              >  6:114272/1‑140 (MQ=255)
                                                                                                                       cAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTg                              <  5:114272/140‑1 (MQ=255)
                                                                                                                               aGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTctctct             <  1:260597/149‑1 (MQ=255)
                                                                                                                                 gCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTa            >  3:136252/1‑148 (MQ=255)
                                                                                                                                 gCTTGAGATTGAGGCTGCGGCACAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTAt           >  1:100427/1‑149 (MQ=255)
                                                                                                                                  ctagagaTTGAGGCTGCGGCGCAGCCCTTCGCGGGCAAGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACtttt                      >  6:264413/4‑137 (MQ=255)
                                                                                                                                  cTTGAGATTGTGGCTGGGGATCAGCCCTTAGCGTGCATGCCGTCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCTCTACATTTGATCCCCTGTGCCAGCAAAAGCTTGCACtttt                      <  5:264413/137‑1 (MQ=255)
                                                                                                                                  cTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTAtc          >  7:110130/1‑149 (MQ=255)
                                                                                                                                     gagaTTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTAtctc        >  4:184260/1‑148 (MQ=255)
                                                                                                                                        aTTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTCCCCTTTGCTTTTACGGCCCCACCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGAACTATTTATATCTCTAtaaca       >  2:198086/1‑142 (MQ=255)
                                                                                                                                          aagggCTGCGCCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGAc  >  5:208876/4‑149 (MQ=255)
                                                                                                                                           gAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCa       >  6:84461/1‑143 (MQ=255)
                                                                                                                                           gAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCa       <  5:84461/143‑1 (MQ=255)
                                                                                                                                           agggCTGCGCCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCtgt                                            <  2:109812/104‑1 (MQ=255)
                                                                                                                                           agggCTGCGCCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCtgt                                            >  1:109812/3‑106 (MQ=255)
                                                                                                                                             ggCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCaaaa                                  >  7:129438/1‑114 (MQ=255)
                                                                                                                                             ggCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCaaaa                                  <  8:129438/114‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GGGCTGGGGCTCGCAATCGTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGAC  >  NC_002947/6067050‑6067336

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