Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,863,461 G→C 7.8% S143W (TCG→TGG)  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,4610GC 7.8% 154.7 / 3.1 51S143W (TCG→TGG) PP_5138LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base C (2/2);  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 = 9.70e-01

CACCACCTCGATCGACACCGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATGTAGTAAGGCGCAG  >  NC_002947/5863343‑5863599
                                                                                                                      |                                                                                                                                          
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                                                                                                                            >  3:227654/1‑127 (MQ=255)
      cTCGATCGACACCGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCCGCGTCGTCGCAATggg                                                                                                                            <  4:227654/129‑3 (MQ=255)
             gaCACCGACTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGa                                                                                                           <  3:254342/139‑1 (MQ=255)
             gaCACCGACTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGa                                                                                                           >  4:254342/1‑139 (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                                                                      >  5:295466/1‑142 (MQ=255)
                                               ggTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCAc                                                                      <  6:295466/142‑1 (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)
                                                                                                                      |                                                                                                                                          
CACCACCTCGATCGACACCGTCTGGCGATGGGCCAGCGGATGGTTGCGGTGCACCGCCACCACCAAAGGGTCAGTGCCCAGTATCCGCCGCACCAACCGTGCGTCCTCCACCCATTGCGACGACGCGGCAATGTCCACCCGGTAATCCTCGAGCATTTCCAGCACCTGCTGCGAGTTGCCGATCTCCACCGCCACTTCAACCTGCGGCAGGCGTTCGCGATAGATCTTGACCAGGTCGAGAATGTAGTAAGGCGCAG  >  NC_002947/5863343‑5863599

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