Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,858,150 G→C 10.7% V67V (GTG→GTC)  PP_4273 → MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,858,1500GC10.7% 164.7 / 9.8 56V67V (GTG→GTC) PP_4273MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/24);  new base C (2/4);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

CTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGCCCGAAATTCCA  >  NC_002947/4858006‑4858293
                                                                                                                                                |                                                                                                                                               
ctacGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGAc                                                                                                                                              <  7:4284/148‑1 (MQ=255)
 tacCAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGCCCAGCGTCGACAGTATGCAGACGATGTGGCAGACCGGTTCCCGCTGATCCCCCTGGGGCTGGCTGCGCGGGAAACGCTGGATGAACTGGGGGCGATGCTGGTGGACCt                                                                                                                                            >  4:170620/1‑149 (MQ=255)
    gagaAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTg                                                                                                                                           >  8:181816/1‑147 (MQ=255)
        agaaAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTc                                                                                                                                     <  1:180543/149‑1 (MQ=255)
         gaaAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCg                                                                                                                                    <  2:208990/149‑1 (MQ=255)
          aaaGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTACGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGa                                                                                                                                   >  6:248848/1‑149 (MQ=255)
             ggCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCCGGTCGACCTGCAGGTCCAcca                                                                                                                                <  3:109954/149‑2 (MQ=255)
                   cTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGCTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCCGGTCGACCTGCAGGTCCACCAGCatc                                                                                                                           <  7:109047/148‑3 (MQ=255)
                    tCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGTTCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATg                                                                                                                         >  6:105396/1‑149 (MQ=255)
                    tCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGAcc                                                                                                                                             <  2:37138/129‑1 (MQ=255)
                    tCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGAcc                                                                                                                                             >  1:37138/1‑129 (MQ=255)
                     cAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGc                                                                                                                        >  6:249010/1‑149 (MQ=255)
                        gacgctGGCCGCGCCCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTATCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGa                                                                                                                     >  5:337193/6‑149 (MQ=255)
                         tCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGAt                                                                                                                    >  3:2451/1‑149 (MQ=255)
                           cTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATgg                                                                                                                  >  5:247980/1‑149 (MQ=255)
                              ggCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGTTCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGcc                                                                                                                <  5:105396/148‑1 (MQ=255)
                              ggCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGc                                                                                                                 <  7:181816/147‑1 (MQ=255)
                                  gggCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTACGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGc                                                                                                           <  5:248848/149‑1 (MQ=255)
                                           cAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCt                                                                                                  <  5:63630/149‑1 (MQ=255)
                                           cAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCCGGTCGACCTGCAGGTCCACCAGc                                                                                                                              <  8:139277/121‑1 (MQ=255)
                                           cAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCCGGTCGACCTGCAGGTCCACCAGc                                                                                                                              >  7:139277/1‑121 (MQ=255)
                                                tCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATg                                                                                             <  4:2451/149‑1 (MQ=255)
                                                   aCAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATgcgc                                                                                          <  6:225378/149‑1 (MQ=255)
                                                    cAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGcgcg                                                                                         <  5:249010/149‑1 (MQ=255)
                                                          gCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATc                                                                                     >  4:272421/1‑147 (MQ=255)
                                                            aGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCCGGTCGACCTGCAGGTCCACCAGCa                                                                                                                             <  3:207214/105‑1 (MQ=39)
                                                            aGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCCGGTCGACCTGCAGGTCCACCAGCa                                                                                                                             >  4:207214/1‑105 (MQ=39)
                                                             gCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCg                                                                                >  7:229850/1‑149 (MQ=255)
                                                              ccGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGAt                                                                                      >  1:205882/1‑142 (MQ=255)
                                                              ccGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGAt                                                                                      <  2:205882/142‑1 (MQ=255)
                                                               cGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGt                                                                              <  1:136338/149‑1 (MQ=255)
                                                                    tgGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTg                                                                            >  2:53773/1‑146 (MQ=255)
                                                                    tgGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCa                                                                          <  1:53773/148‑1 (MQ=255)
                                                                      gCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGCTCAAGCGGTTGCAGGc                                                                       >  6:321305/1‑149 (MQ=255)
                                                                      gCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGc                                                                       >  8:252619/1‑149 (MQ=255)
                                                                      gCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGc                                                                       >  6:55882/1‑149 (MQ=255)
                                                                              gTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGc                                                                >  5:320073/1‑148 (MQ=255)
                                                                                        aTCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCg                                                     <  2:41168/149‑1 (MQ=255)
                                                                                              cTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGc                                               <  6:320073/149‑1 (MQ=255)
                                                                                                 gggCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGcc                                              <  6:247980/147‑1 (MQ=255)
                                                                                                   gCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGc                                                      <  4:280146/137‑1 (MQ=255)
                                                                                                        cTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGc                                     >  2:84928/1‑149 (MQ=255)
                                                                                                           cGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGcc                                    <  8:167741/147‑1 (MQ=255)
                                                                                                            ggggTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGTTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCt                                  >  1:310067/1‑148 (MQ=255)
                                                                                                                ttttCGCTGGATGAAGGGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGTCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGaa                             <  3:170620/149‑1 (MQ=255)
                                                                                                                   tCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTgaga                              >  6:114252/1‑145 (MQ=255)
                                                                                                                    cGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGa                          <  3:272421/148‑1 (MQ=255)
                                                                                                                        ggATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGatcatc                     >  5:342606/1‑149 (MQ=255)
                                                                                                                                tGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGc             >  6:220822/1‑149 (MQ=255)
                                                                                                                                   gggCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGccc           >  6:345775/1‑148 (MQ=255)
                                                                                                                                    ggCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGc             <  3:254104/145‑1 (MQ=255)
                                                                                                                                    ggCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGc             >  4:254104/1‑145 (MQ=255)
                                                                                                                                     gCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGCCCGaa        <  5:220822/149‑1 (MQ=255)
                                                                                                                                       gATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGCCCGAAAt      >  5:72567/1‑149 (MQ=255)
                                                                                                                                         tGCTGGTGGACCTGCAGGTCGAACGGCAGATGCTGATCGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGCTTGCGGCATGCCGCTCAATGCCCTGAGAAAGATCATCTGGCATGCCCGAAATTc    >  7:285544/1‑149 (MQ=255)
                                                                                                                                           cTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGCCCGAAATTCca  <  8:229850/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGCCCGAAATTCCA  >  NC_002947/4858006‑4858293

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑