Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,863,469 G→C 11.8% A140A (GCC→GCG)  PP_5138 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,863,4690GC11.8% 136.7 / 10.1 51A140A (GCC→GCG) PP_5138LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (21/24);  new base C (3/3);  total (24/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 = 8.28e-01

CACCACCTCGATCGACACCGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATGTAGTAAGGCGCAGTGGCAGCAATGCGC  >  NC_002947/5863343‑5863613
                                                                                                                              |                                                                                                                                                
caccacCTCGATCGACACCGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAATATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCt                                                                                                                            <  5:128055/149‑1 (MQ=255)
   tatagagatctaCAGCGTCTGGCGATGGGCCAGCGGATGGTTGCTGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCGCGa                                                                                                                         <  8:293859/143‑1 (MQ=255)
      cTCGATCGACACCGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCCGCGTCGTCGCAATggg                                                                                                                                          <  6:297636/129‑3 (MQ=255)
      cTCGATCGACACCGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCCGCGTCGTCGCAATggg                                                                                                                                          >  5:297636/1‑127 (MQ=255)
      cTCGATCGACACCGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCCGCGTCGTCGCAATggg                                                                                                                                          <  4:227654/129‑3 (MQ=255)
      cTCGATCGACACCGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCCGCGTCGTCGCAATggg                                                                                                                                          >  3:227654/1‑127 (MQ=255)
             gaCACCGACTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGa                                                                                                                         >  4:254342/1‑139 (MQ=255)
             gaCACCGACTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGa                                                                                                                         <  3:254342/139‑1 (MQ=255)
               cacCGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTccacca                                                                                                             <  3:303171/149‑1 (MQ=255)
                  cGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACct                                                                                                          >  5:128020/1‑149 (MQ=255)
                   gTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGc                                                                                                                       <  1:302546/135‑1 (MQ=255)
                   gTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGc                                                                                                                       >  2:302546/1‑135 (MQ=255)
                   gTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACctg                                                                                                         >  1:15050/1‑149 (MQ=255)
                                 cAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTcc                                                                                                                 <  2:121636/127‑1 (MQ=255)
                                 cAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTcc                                                                                                                 >  1:121636/1‑127 (MQ=255)
                                 cAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGa                                                                                           >  4:75373/1‑149 (MQ=255)
                                  ctttggtggtggCGGTGCACCGCAACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGa                                                                                           >  8:149450/7‑148 (MQ=255)
                                      gATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGc                                                                                              >  3:310109/1‑141 (MQ=255)
                                      gATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGc                                                                                              <  4:310109/141‑1 (MQ=255)
                                           ttGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGc                                                                                 <  1:219624/149‑1 (MQ=255)
                                              cGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCCCCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCAc                                                                              >  3:4762/1‑149 (MQ=255)
                                               ggTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCAc                                                                                    <  6:295466/142‑1 (MQ=255)
                                               ggTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCAc                                                                                    >  5:295466/1‑142 (MQ=255)
                                                          accaccaAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTTCGTCCTTCACCCATTTCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTTCTGCTAGTTTCCGATCTCCACCGCCACTTCAACCTGc                                                                    <  7:52078/147‑1 (MQ=255)
                                                          accaccaAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATACTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGc                                                                    >  8:52078/1‑147 (MQ=255)
                                                            caccaAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCgg                                                                  <  4:329872/147‑1 (MQ=255)
                                                            caccaAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCgg                                                                  >  3:329872/1‑147 (MQ=255)
                                                             accaAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAg                                                               >  4:322639/1‑149 (MQ=255)
                                                                aaaGGGTCAGTGCCCAGTATCCGACGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATTTACCCCCGGTAATCCACCAGCATATCCAGCACCTGATGCGAGTTTCCGATCTCCACCGCCACTTCAACCTGCGGAAGGCg                                                            >  8:132081/1‑149 (MQ=255)
                                                                          tGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGa                                                     >  1:210731/1‑146 (MQ=255)
                                                                          tGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGa                                                     >  6:148395/1‑146 (MQ=255)
                                                                          tGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATa                                                   <  5:148395/148‑1 (MQ=255)
                                                                          tGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAg                                                  <  2:15050/149‑1 (MQ=255)
                                                                              cAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCCGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCt                                              >  1:64152/1‑149 (MQ=255)
                                                                                gTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCtt                                             >  7:81250/1‑148 (MQ=255)
                                                                                   tCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCCGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGAc                                          <  2:64152/148‑1 (MQ=255)
                                                                                       ccgcACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGt                                     <  7:149450/149‑1 (MQ=255)
                                                                                          cACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCga                                  <  8:81250/149‑1 (MQ=255)
                                                                                              aaCCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAAt                              <  3:322639/149‑1 (MQ=255)
                                                                                                 cGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATgta                           <  6:128020/149‑1 (MQ=255)
                                                                                                 cGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATgta                           <  4:4762/149‑1 (MQ=255)
                                                                                                 cGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATgta                           <  3:75373/149‑1 (MQ=255)
                                                                                                     cGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATgtagta                        <  2:210731/148‑1 (MQ=255)
                                                                                                      gTCCTCCACCCATTGCGACGACGCGGCAATGTGCACCCGGTAATGCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATGTAGTAAg                      <  7:132081/149‑1 (MQ=255)
                                                                                                        cctccACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGt                                                           >  6:72126/1‑110 (MQ=255)
                                                                                                        cctccACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGt                                                           <  5:72126/110‑1 (MQ=255)
                                                                                                          tccACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATgtag                          <  7:268653/141‑1 (MQ=255)
                                                                                                          tccACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATgtag                          >  8:268653/1‑141 (MQ=255)
                                                                                                            cACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATGTAGTAAGGCGCat                >  5:320321/1‑148 (MQ=255)
                                                                                                               ccATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCgaga                                >  7:260880/1‑130 (MQ=255)
                                                                                                               ccATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCgaga                                <  8:260880/130‑1 (MQ=255)
                                                                                                                              ggCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATGTAGTAAGGCGCAGTGGCAGCAATgcgc  <  7:107704/145‑1 (MQ=255)
                                                                                                                              ggCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATGTAGTAAGGCGCAGTGGCAGCAATgcgc  >  8:107704/1‑145 (MQ=255)
                                                                                                                              |                                                                                                                                                
CACCACCTCGATCGACACCGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATGTAGTAAGGCGCAGTGGCAGCAATGCGC  >  NC_002947/5863343‑5863613

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