Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,088,664 C→G 6.7% intergenic (+160/‑94) pmbA → / → PP_0943 regulatory protease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,088,6640CG 6.7% 188.0 / 3.9 60intergenic (+160/‑94)pmbA/PP_0943regulatory protease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/29);  new base G (2/2);  total (29/31)
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.55e-01

GCCCTTTCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGC  >  NC_002947/1088521‑1088802
                                                                                                                                                |                                                                                                                                          
gCCCTTTCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCgg                                                                                                                                       <  7:274323/149‑1 (MQ=255)
         cgGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAAACCCGGCGCGAGGGCCGGGtcttcgga                                                                                                                              <  1:145359/149‑5 (MQ=37)
              aaaCCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCAGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATACGAAGACCCGGACCTCGCGCCGGGTTTTTATTtcat                                                                                                                          >  8:16987/1‑148 (MQ=255)
               aaCCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTtcatca                                                                                                                        <  2:215686/149‑1 (MQ=255)
               aaCCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTtcatc                                                                                                                         <  2:115295/148‑1 (MQ=255)
                  ccGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTGAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCAcc                                                                                                                      >  1:241616/1‑149 (MQ=255)
                       ccTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGtttt                                                                                                                                  <  1:307416/131‑1 (MQ=255)
                       ccTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGtttt                                                                                                                                  >  2:307416/1‑131 (MQ=255)
                        cTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTACAAATCCGAAGACCCGGCGCGAGGGCCGGGtcttc                                                                                                                                 <  4:146020/131‑2 (MQ=25)
                        cTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTACAAATCCGAAGACCCGGCGCGAGGGCCGGGtcttc                                                                                                                                 >  3:146020/1‑130 (MQ=25)
                          acacGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGttt                                                                                                                                   <  2:54478/127‑1 (MQ=255)
                          acacGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGttt                                                                                                                                   >  1:54478/1‑127 (MQ=255)
                           cacGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTg                                                                                                             <  2:317250/148‑1 (MQ=255)
                              ggtggtGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGt                                                                                                         >  3:174886/1‑149 (MQ=255)
                                     gACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTa                                                                                                  >  7:237510/1‑149 (MQ=255)
                                                gggTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTta                                                                                       <  2:216131/149‑1 (MQ=255)
                                                 ggTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTAATCACCTTTTTCTTGAAGTGATTCTATTTATCACTtaa                                                                                      >  5:206698/1‑149 (MQ=255)
                                                       acacAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATg                                                                                 >  7:23430/1‑148 (MQ=255)
                                                            ggCCGTTGTAGGTGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATc                                                                           <  2:256841/149‑1 (MQ=255)
                                                                gTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTa                                                                       >  1:189889/1‑149 (MQ=255)
                                                                     tcttcgCGGGTGAACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATc                                                                     >  7:152396/6‑147 (MQ=255)
                                                                        gAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACt                                                                                         <  8:93771/124‑1 (MQ=255)
                                                                        gAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACt                                                                                         >  7:93771/1‑124 (MQ=255)
                                                                            gggTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAAccc                                                            <  2:241616/149‑1 (MQ=255)
                                                                                   cccGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCTGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCgatg                                                     <  7:16987/149‑1 (MQ=255)
                                                                                    ccGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTtaa                                                                                      >  6:16625/1‑115 (MQ=255)
                                                                                    ccGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTtaa                                                                                      <  5:16625/115‑1 (MQ=255)
                                                                                    ccGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCgatga                                                    >  7:64783/1‑149 (MQ=255)
                                                                                    ccGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCgatga                                                    >  5:130107/1‑149 (MQ=255)
                                                                                       cgAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGt                                                 <  4:174886/149‑1 (MQ=255)
                                                                                              gCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGa                                                                                                       >  1:243432/1‑88 (MQ=255)
                                                                                              gCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGa                                                                                                       <  2:243432/88‑1 (MQ=255)
                                                                                               ccGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGaaa                                         <  8:152396/149‑1 (MQ=255)
                                                                                                cGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAAt                                        >  3:310572/1‑149 (MQ=255)
                                                                                                 cgAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCCCGCATGGTTTTTATTTCATCACCTTTTTCTTGAAGTCATTCTATTTATCAATTAATAATGAATATAATTATCCAGTACCCCGACGATGATGTTCATGAAATc                                       >  4:62826/2‑149 (MQ=255)
                                                                                                       gCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGtcct                                 <  6:195122/149‑1 (MQ=255)
                                                                                                        cACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTg                                                                                                             >  2:20373/1‑72 (MQ=255)
                                                                                                        cACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTg                                                                                                             <  1:20373/72‑1 (MQ=255)
                                                                                                        cACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGtcctc                                <  8:64783/149‑1 (MQ=255)
                                                                                                           tACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCgatg                                                     >  2:250814/1‑125 (MQ=255)
                                                                                                           tACACCTTTCCCTGCAAATCCGAAGACCCGGCCATCGCGCCCGGTTTTTATTTCAGCACCTTTTTCTTGAAGTTATTCTATTAAACACTTAATAAAAAATAACATTATCCATTAAGCCGTCGATCATGTTCATTAAAACCGTCCTCcgc                             >  8:188565/1‑149 (MQ=255)
                                                                                                                cTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGTTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACt                        <  8:237510/149‑1 (MQ=255)
                                                                                                                cTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCTCTTAATAATGAATATCATTATCCAGTAACCCgg                                                          <  2:270852/115‑1 (MQ=255)
                                                                                                                cTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCTCTTAATAATGAATATCATTATCCAGTAACCCgg                                                          >  1:270852/1‑115 (MQ=255)
                                                                                                                cTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACt                        >  6:201814/1‑149 (MQ=255)
                                                                                                                        gCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCcgcg                            >  5:263797/1‑137 (MQ=255)
                                                                                                                        gCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCcgcg                            <  6:263797/137‑1 (MQ=255)
                                                                                                                         cAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTAc                 <  5:201814/147‑1 (MQ=255)
                                                                                                                         cAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTAc                 <  4:201072/147‑1 (MQ=255)
                                                                                                                          aaaTCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTg              <  2:189889/149‑1 (MQ=255)
                                                                                                                          aaaTCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTg              <  6:130107/149‑1 (MQ=255)
                                                                                                                          aaaTCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGATCATGAAATCCGTCCTCGGCGAACTGCCCTCCCTg              >  8:123986/1‑149 (MQ=255)
                                                                                                                           aaTCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTgg             <  7:67547/149‑1 (MQ=255)
                                                                                                                             tCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGaa           <  4:310572/149‑1 (MQ=255)
                                                                                                                              ataaaaaCCCGGCGCGAGGGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGaa           >  8:267758/7‑148 (MQ=255)
                                                                                                                                   gACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTgg     >  4:257537/1‑149 (MQ=255)
                                                                                                                                   gACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTgg     <  3:62826/149‑1 (MQ=255)
                                                                                                                                    aCCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGgc    <  1:320202/149‑1 (MQ=255)
                                                                                                                                     cccGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGgcg   >  6:263726/1‑149 (MQ=255)
                                                                                                                                      ccGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCATCCGCGAACTGCCCTACCTGGAAAACTGgcgc  >  2:237981/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                          
GCCCTTTCGCGGGTAAACCCGCTCCTACACGGTGGTGGACGTACATCGGGGTTGGACACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGCACTACACCTTTCCCTGCAAATCCGAAGACCCGGCCCTCGCGCCGGGTTTTTATTTCATCACCTTTTTCTTGAAGTGATTCTATTTATCACTTAATAATGAATATCATTATCCAGTAACCCGGCGATGATGTTCATGAAATCCGTCCTCCGCGAACTGCCCTACCTGGAAAACTGGCGC  >  NC_002947/1088521‑1088802

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: