Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,858,167 G→A 9.8% R73Q (CGG→CAG)  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,1670GA 9.8% 135.7 / 5.4 51R73Q (CGG→CAG) PP_4273MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base A (2/3);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

CTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGCCCGAAATTCCAGCAGCTGATGAAGCTGTCGG  >  NC_002947/4858025‑4858313
                                                                                                                                              |                                                                                                                                                  
cTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGCTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCCGGTCGACCTGCAGGTCCACCAGCatc                                                                                                                                               <  7:109047/148‑3 (MQ=255)
 tCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGTTCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATg                                                                                                                                             >  6:105396/1‑149 (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)
                                                                                                                              gACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGATGCAGGCGATGAGCAAGGGGATGCGGAATGCCGCGAAATGCCATGAGAAAGATCATCTGGCATGCCCGAAATTCCAGAAGCt                >  7:162364/1‑149 (MQ=255)
                                                                                                                                   gCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGAGCATCTGGCATGCCCGAAATTCCAGCAGCtgatg            >  8:168859/1‑148 (MQ=255)
                                                                                                                                     aGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGCCCGAAATTCCAGCAGCTGATGAAg         <  1:84928/149‑1 (MQ=255)
                                                                                                                                            ccGGCAGATGCTGATGCCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGCCCGAAATTCCAGCAGCTGATGAAGCTGTcgg  <  7:252619/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGATGTGGCAGACCGGTTGGCGCTGATCGCCCTGGGGCAGGCTGCGGGGTTTTCGCTGGATGAAGTGGGGGCGATGCTGGTGGACCTGCAGGTCGACCGGCAGATGCTGATGGCCAAGGCGGACGAGCTGGATGCGCGGATCAAGCGGTTGCAGGCGATGAGCAAGGGGTTGCGGCATGCCGCGCAATGCCCTGAGAAAGATCATCTGGCATGCCCGAAATTCCAGCAGCTGATGAAGCTGTCGG  >  NC_002947/4858025‑4858313

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