Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,869,903 C→A 17.0% E354* (GAG→TAG) ‡ PP_3419 ← sigma‑54 dependent transcriptional regulator/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,869,9030CA17.0% 132.3 / 16.2 53E354* (GAG→TAG) ‡PP_3419sigma‑54 dependent transcriptional regulator/response regulator
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base A (4/5);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

CGGCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAGGTCACGCCTG  >  NC_002947/3869757‑3870051
                                                                                                                                                   |                                                                                                                                                    
cGGCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTtcggcctcggc                                                                                                                                                   <  1:1109126/150‑1 (MQ=255)
    tCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGaa                                                                                                                                               >  2:161734/1‑150 (MQ=255)
      gTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTTCGGCCGAGGCCGaaa                                                                                                                                              >  1:1144424/1‑150 (MQ=255)
      gTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGaaaa                                                                                                                                             >  1:1131676/1‑150 (MQ=255)
       tCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCGAGGCCGaaa                                                                                                                                              <  1:1067626/148‑1 (MQ=255)
       tCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCGAGGCCGaaa                                                                                                                                              >  2:1067626/1‑148 (MQ=255)
        cAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAgc                                                                                                                                           <  2:1090925/150‑1 (MQ=255)
         aaCCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAgcg                                                                                                                                          >  2:982598/1‑150 (MQ=255)
            cAAGGTCATGCACAGATTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTTCGGCCGAGGCCGaaagcctgc                                                                                                                                        <  2:838073/150‑7 (MQ=255)
             aaGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTTCGGCCGAGGCCGaaagcctgct                                                                                                                                       <  2:1144424/150‑8 (MQ=255)
                           gTTCAACTGATGCGCTTGCCCCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATa                                                                                                                        <  2:793389/150‑1 (MQ=255)
                           gTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATa                                                                                                                        >  1:249623/1‑150 (MQ=255)
                           gTTCAACTGATGCGCTTGCACCACCACAAGCTGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATa                                                                                                                        >  1:31509/1‑150 (MQ=255)
                            ttCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAg                                                                                                                       >  1:113138/1‑150 (MQ=255)
                             tCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAgc                                                                                                                      <  1:982598/150‑1 (MQ=255)
                                      gcgcTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAt                                                                                                             <  2:241787/150‑1 (MQ=255)
                                      gcgcTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAt                                                                                                             >  2:852742/1‑150 (MQ=255)
                                                  ccTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGa                                                                                                 >  1:829032/1‑150 (MQ=255)
                                                    tCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGaaa                                                                                               <  1:852742/150‑1 (MQ=255)
                                                       ttCGGCGCCAGCAACACCGTCTGCGCCAGCATGTTGCACAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCaa                                                                                                    >  2:322450/1‑142 (MQ=255)
                                                       tgcGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCACAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCaa                                                                                                    <  1:322450/140‑1 (MQ=255)
                                                                aGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGt                                                                                   >  2:611714/1‑150 (MQ=255)
                                                                    acacCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCa                                                                               <  1:386546/150‑1 (MQ=255)
                                                                    acacCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCa                                                                               >  2:1031453/1‑150 (MQ=255)
                                                                        cGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGAt                                                                           <  1:1031453/150‑1 (MQ=255)
                                                                        cGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGAt                                                                           <  2:23159/150‑1 (MQ=255)
                                                                        cGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGAt                                                                           <  2:829032/150‑1 (MQ=255)
                                                                                            gcAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACa                                                       >  1:948783/1‑150 (MQ=255)
                                                                                                  cACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGcgcg                                                 >  2:179293/1‑150 (MQ=255)
                                                                                                    cGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTAGTACAAAAAGCTCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGc                                               <  2:31509/150‑1 (MQ=255)
                                                                                                         aTTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCtt                                          >  2:786585/1‑150 (MQ=255)
                                                                                                                  ccagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAg                                                                              <  2:665151/105‑1 (MQ=255)
                                                                                                                  ccagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAAGGTGGGCAAGGAAATGCCGCGCCAGTAGCAg                                                                              >  1:665151/1‑105 (MQ=255)
                                                                                                                   cagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCcgc                                                                                         >  1:1147688/1‑93 (MQ=255)
                                                                                                                   cagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCcgc                                                                                         <  2:1147688/93‑1 (MQ=255)
                                                                                                                    agCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATaa                               >  2:835585/1‑150 (MQ=255)
                                                                                                                    agCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATaa                               <  1:302151/150‑1 (MQ=255)
                                                                                                                    agCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCAATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATaa                               >  2:767704/1‑150 (MQ=255)
                                                                                                                     gCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCTGCACCTTCAGGGCGATAAt                              <  2:249623/150‑1 (MQ=255)
                                                                                                                            ccGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTg                                                                     <  2:148454/104‑1 (MQ=255)
                                                                                                                            ccGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTg                                                                     >  1:148454/1‑104 (MQ=255)
                                                                                                                             cGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGc                      <  2:1131676/150‑1 (MQ=255)
                                                                                                                              gAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCg                     >  2:534189/1‑150 (MQ=255)
                                                                                                                              gAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCg                     <  1:786585/150‑1 (MQ=255)
                                                                                                                                 cagcagGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaa                  <  1:1060244/150‑1 (MQ=255)
                                                                                                                                 cagcagGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaa                  <  1:534189/150‑1 (MQ=255)
                                                                                                                                 cagcagGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaa                  <  1:161734/150‑1 (MQ=255)
                                                                                                                                    cagGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaagaa               <  2:180400/150‑1 (MQ=255)
                                                                                                                                     gcgcttttCGGCCGAGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAggg                                     >  2:334719/6‑128 (MQ=255)
                                                                                                                                     gcgcttttCGGCCGAGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAggg                                     <  1:334719/123‑1 (MQ=255)
                                                                                                                                      ggCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACa             >  1:916938/1‑150 (MQ=255)
                                                                                                                                       gcttttCGGCCGAGGCCGAAAGGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACa             >  2:171217/4‑150 (MQ=255)
                                                                                                                                       gCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGc                                               >  2:429243/1‑115 (MQ=255)
                                                                                                                                       gCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGc                                               <  1:429243/115‑1 (MQ=255)
                                                                                                                                                  tcggccGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAGGTCACGCCTg  <  2:948783/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAGGTCACGCCTG  >  NC_002947/3869757‑3870051

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