Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,240,105 G→C 6.9% A62G (GCC→GGC)  PP_3717 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,240,1050GC 6.9% 183.7 / 4.3 58A62G (GCC→GGC) PP_3717LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base C (2/2);  total (29/29)
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.26e-01

GCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAG  >  NC_002947/4239962‑4240246
                                                                                                                                               |                                                                                                                                             
gcagGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCg                                                                                                                                          >  7:238053/1‑149 (MQ=255)
gcagGTAGCCGGCGCTGAAATAGCGCGCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCTGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCTACTCCACCTGCACCAACTGGCCGATGGCAGCg                                                                                                                                          <  4:264700/149‑1 (MQ=255)
     tAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCcac                                                                                                                                     <  1:16787/149‑1 (MQ=255)
         cGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATAc                                                                                                                                 >  5:207924/1‑149 (MQ=255)
              cTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAg                                                                                                                            <  2:69910/149‑1 (MQ=255)
               tGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGCTGCCATCGGCcagttg                                                                                                                                  <  4:85541/142‑5 (MQ=255)
               tGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGCTGCCATCGGCcagttg                                                                                                                                  >  3:85541/1‑138 (MQ=255)
                 aaaTAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAgc                                                                                                                           <  7:85336/147‑1 (MQ=255)
                        gcTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCgggg                                                                                                                  <  6:191419/149‑1 (MQ=255)
                        gcTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCGCgggg                                                                                                                  >  6:165010/1‑149 (MQ=255)
                         cTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCggggg                                                                                                                 >  6:143/1‑149 (MQ=255)
                                  cgcgGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCAt                                                                                                        <  8:238053/149‑1 (MQ=255)
                                   gcgGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATg                                                                                                       <  1:255132/149‑1 (MQ=255)
                                                      ggATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAc                                                                                    <  8:292813/149‑1 (MQ=255)
                                                            cTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGcc                                                                                                                                       <  7:30146/92‑1 (MQ=255)
                                                            cTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGcc                                                                                                                                       >  8:30146/1‑92 (MQ=255)
                                                            cTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCgg                                                                                                                    >  1:108200/1‑111 (MQ=255)
                                                            cTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCgg                                                                                                                    <  2:108200/111‑1 (MQ=255)
                                                            cTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCa                                                                                >  6:265709/1‑147 (MQ=255)
                                                            cTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCa                                                                                <  5:265709/147‑1 (MQ=255)
                                                            cTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACg                                                                              >  4:287074/1‑149 (MQ=255)
                                                            cTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACg                                                                              >  2:42131/1‑149 (MQ=255)
                                                              tGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTg                                                                            >  5:14130/1‑149 (MQ=255)
                                                               gCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCCCCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGc                                                                           >  8:76560/1‑149 (MQ=255)
                                                                 tGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATgg                                                                                                      >  1:263046/1‑120 (MQ=255)
                                                                 tGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATgg                                                                                                      <  2:263046/120‑1 (MQ=255)
                                                                    tACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCActg                                                                      >  8:76152/1‑149 (MQ=255)
                                                                     aCAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCCCGTGCCActgc                                                                     >  7:152545/1‑149 (MQ=255)
                                                                         cAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCActgctg                                                                   >  2:240407/1‑147 (MQ=255)
                                                                         cAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCa                                                                 <  3:287074/149‑1 (MQ=255)
                                                                             tGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAgcgc                                                             >  7:204800/1‑149 (MQ=255)
                                                                                ggcgCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGcc                                                          <  7:76152/149‑1 (MQ=255)
                                                                                   gcccgcccCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCAt                                                       >  6:25749/1‑149 (MQ=255)
                                                                                        ccccTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGa                                                  <  5:25749/149‑1 (MQ=255)
                                                                                        ccccTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGa                                                  <  1:240407/149‑1 (MQ=255)
                                                                                         cccTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAg                                                 <  7:260723/149‑1 (MQ=255)
                                                                                          ccTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAgg                                                <  1:42131/149‑1 (MQ=255)
                                                                                                aCTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAgg                                                <  1:99379/143‑1 (MQ=255)
                                                                                                aCTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAgg                                                >  2:99379/1‑143 (MQ=255)
                                                                                                              atcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGga                              <  8:29514/147‑1 (MQ=255)
                                                                                                              atcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGga                              >  7:29514/1‑147 (MQ=255)
                                                                                                               tcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGa                                                  <  2:127251/126‑1 (MQ=255)
                                                                                                               tcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGa                                                  >  1:123355/1‑126 (MQ=255)
                                                                                                               tcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGa                                                  >  1:127251/1‑126 (MQ=255)
                                                                                                               tcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGagag                           <  8:204800/149‑1 (MQ=255)
                                                                                                                  tcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGa                                                  <  7:237324/123‑1 (MQ=255)
                                                                                                                  tcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGa                                                  >  8:237324/1‑123 (MQ=255)
                                                                                                                         cACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAgc                                                               <  2:146908/103‑1 (MQ=255)
                                                                                                                         cACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAgc                                                               >  1:146908/1‑103 (MQ=255)
                                                                                                                          aCCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCAt                >  3:89346/1‑149 (MQ=255)
                                                                                                                           ccTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATg               <  6:207924/149‑1 (MQ=255)
                                                                                                                            cTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGc              >  6:131577/1‑149 (MQ=255)
                                                                                                                                aCCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCgg          <  6:14130/149‑1 (MQ=255)
                                                                                                                                 ccAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCgg          >  1:291755/1‑148 (MQ=255)
                                                                                                                                  tatgtGGCCGCTGCCATCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGga                              >  5:168845/5‑127 (MQ=255)
                                                                                                                                  tatgtGGCCGCTGCCATCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGga                              <  6:168845/123‑1 (MQ=255)
                                                                                                                                     caggCCCATGGAAGCGGCCACATACTCCAGCGCCTCGGCCGTATCACCATCGGCCAGGGCGCGCTGCAACGCCACGTGCCACTCCTGCAGCGCGCCCATAGTGAGCGTAAGCGGGGTATCGGCAGAGGCCTGCGTCATGCCCGGCAGt      >  4:105211/3‑148 (MQ=255)
                                                                                                                                        gCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAg  <  7:76560/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
GCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAG  >  NC_002947/4239962‑4240246

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: