Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,267,833 T→G 5.4% L364W (TTG→TGG)  atoS → sensory histidine kinase in two‑component regulatory system with AtoC

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,267,8330T→G 5.4% 2.4 0.0 0.0551091atoSsensory histidine kinase in two‑component regulatory system with AtoC
Reads supporting (aligned to +/- strand):  new base (1/2):  ref base (18/34):  total (19/36)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGATTTTCTCTGATTTAACTGCCCGCAAAGAA  >  REL606/2267799‑2267867
                                  |                                  
aTTCATAACACGCACGGTGAAATGATAGGTGCTTTg                                   >  1:4735254/1‑36 (MQ=255)
 ttCATAACACGCACGGTGAAATGATAGGTGCTTTgg                                  <  1:5713337/36‑1 (MQ=255)
 ttCATAACACGCACGGTGAAATGATAGGTGCTTTgg                                  <  1:4973567/36‑1 (MQ=255)
  tCATAACACGCACGGTGAAATGATAGGTGCTTTGGt                                 >  1:858821/1‑36 (MQ=255)
  tCATAACACGCACGGTGAAATGATAGGTGCTTTGGt                                 <  1:7412356/36‑1 (MQ=255)
  tCATAACACGCACGGTGAAATGATAGGTGCTTTGGt                                 <  1:4330733/36‑1 (MQ=255)
  tCATAACACGCACGGTGAAATGATAGGTGCTTTGGt                                 <  1:2858704/36‑1 (MQ=255)
   cATAACACGCACGGTGAAATGATAGGTGCTTGGGTg                                <  1:4303012/36‑1 (MQ=255)
    aTAACACGCACGGTGAAATGATAGGTGCTTTGGTGa                               <  1:4540565/36‑1 (MQ=255)
     tAACACGCACGGTGAAATGATAGGTGCTTTGGTGAt                              <  1:736134/36‑1 (MQ=255)
     tAACACGCACGGTGAAATGATAGGTGCTTTGGTGAt                              >  1:3961606/1‑36 (MQ=255)
      aaCACGCACGGTGAAATGATAGGTGCTTTGGTGAtt                             >  1:7869906/1‑36 (MQ=255)
      aaCACGCACGGTGAAATGATAGGTGCTTTGGTGAtt                             <  1:7631060/36‑1 (MQ=255)
      aaCACGCACGGTGAAATGATAGGTGCTTTGGTGAtt                             <  1:3175127/36‑1 (MQ=255)
       acacGCACGGTGAAATGATAGGTGCTTTGGTGAttt                            >  1:6457385/1‑36 (MQ=255)
       acacGCACGGTGAAATGATAGGTGCTTTGGTGAttt                            <  1:1666861/36‑1 (MQ=255)
       acacGCACGGTGAAATGATAGGTGCTTTGGTGAttt                            >  1:4416885/1‑36 (MQ=255)
       acacGCACGGTGAAATGATAGGTGCTTTGGTGAttt                            <  1:1505038/36‑1 (MQ=255)
       acacGCACGGTGAAATGATAGGTGCTTTGGTGAttt                            >  1:1199049/1‑36 (MQ=255)
         acgcacgGTGAAATGATAGGTGCTTTGGTGATTTtc                          <  1:5342102/36‑1 (MQ=255)
         acgcacgGTGAAATGATAGGTGCTTTGGTGATTTtc                          <  1:166558/36‑1 (MQ=255)
         acgcacgGTGAAATGATAGGTGCTTTGGTGATTTtc                          >  1:686023/1‑36 (MQ=255)
         acgcacgGTGAAATGATAGGTGCTTTGGTGATTTtc                          <  1:7109877/36‑1 (MQ=255)
          cgcacgGTGAAATGATAGGTGCTTTGGTGATTTtct                         <  1:2705795/36‑1 (MQ=255)
          cgcacgGTGAAATGATAGGTGCTTTGGTGATTTtct                         <  1:7257945/36‑1 (MQ=255)
           gcacgGTGAAATGATAGGGGCTTTGGTGATTTtctc                        >  1:3859220/1‑36 (MQ=255)
           gcacgGGGAAATGATAGGTGCTTTGGTGATTTtctc                        >  1:2168289/1‑36 (MQ=255)
            cacgGTGAAATGATAGGTGCTTTGGTGATTTTctct                       <  1:3834888/36‑1 (MQ=255)
            cacgGTGAAATGATAGGTGCTTTGGTGATTTTctct                       <  1:6727128/36‑1 (MQ=255)
            cacgGTGAAATGATAGGTGCTTTGGTGATTTTctct                       <  1:1664733/36‑1 (MQ=255)
            cacgGTGAAATGATAGGTGCTTTGGTGATTTTctct                       <  1:4161543/36‑1 (MQ=255)
             acgGTGAAATGATAGGTGCTTTGGTGATTTTCTCTg                      <  1:3321337/36‑1 (MQ=255)
              cgGTGAAATGATAGGTGCTTTGGTGATTTTCTCTGa                     <  1:7609643/36‑1 (MQ=255)
               ggTGAAATGATAGGTGCTTTGGTGATTTTCTCTGAt                    >  1:3192241/1‑36 (MQ=255)
               ggTGAAATGATAGGTGCTTTGGTGATTTTCTCTGAt                    <  1:4987018/36‑1 (MQ=255)
                gTGAAATGATAGGTGCTTTGGTGATTTTCTCTGAtt                   <  1:743910/36‑1 (MQ=255)
                 tGAAATGATAGGTGCTTTGGTGATTTTCTCTGAttt                  <  1:4833869/36‑1 (MQ=255)
                 tGAAATGATAGGTGCTTTGGTGATTTTCTCTGAttt                  <  1:6938265/36‑1 (MQ=255)
                  gAAATGATAGGTGCTTTGGTGATTTTCTCTGATTTa                 <  1:3377359/36‑1 (MQ=255)
                    aaTGATAGGTGCTTTGGTGATTTTCTCTGATTTAAc               <  1:4181684/36‑1 (MQ=255)
                     aTGATAGGTGCTTTGGTGATTTTCTCTGATTTAACt              >  1:4583893/1‑36 (MQ=255)
                     aTGATAGGTGCTTTGGTGATTTTCTCTGATTTAACt              <  1:2783723/36‑1 (MQ=255)
                     aTGATAGGTGCTTTGGTGATTTTCTCTGATTTAACt              >  1:3912977/1‑36 (MQ=255)
                     aTGATAGGTGCTTTGGTGATTTTCTCTGATTTAACt              >  1:366368/1‑36 (MQ=255)
                      tGATAGGTGCTTTGGTGATTTTCTCTGATTTAACTg             <  1:4550805/36‑1 (MQ=255)
                      tGATAGGTGCTTTGGTGATTTTCTCTGATTTAACTg             >  1:3952082/1‑36 (MQ=255)
                         tAGGTGCTTTGGTGATTTTCTCTGATTTAACTTccc          >  1:3958320/1‑36 (MQ=255)
                         tAGGTGCTTTGGTGATTTTCTCTGATTTAACTGccc          <  1:6884331/36‑1 (MQ=255)
                         tAGGTGCTTTGGTGATTTTCTCTGATTTAACTGccc          <  1:4347497/36‑1 (MQ=255)
                          aGGTGCTTTGGTGATTTTCTCTGATTTAACTGCCCg         <  1:5541875/36‑1 (MQ=255)
                          aGGTGCTTTGGTGATTTTCTCTGATTTAACTGCCCg         <  1:1709583/36‑1 (MQ=255)
                           ggTGCTTTGGTGATTTTCTCTGATTTAACTGCCCGc        >  1:3487676/1‑36 (MQ=255)
                           ggTGCTTGGGTGATTTTCTCTGATTTAACTGCCCGc        >  1:6166338/1‑36 (MQ=255)
                           ggTGCTTGGGTGATTTTCTCTGATTTAACTGCCCGc        <  1:2646692/36‑1 (MQ=255)
                             tGCTTTGGTGATTTTCTCTGATTTAACTGCCCGCaa      >  1:5204735/1‑36 (MQ=255)
                                 ttGGTGATTTTCTCTGATTTAACTGCCCGCAAAGaa  >  1:2446606/1‑36 (MQ=255)
                                 ttGGTGATTTTCTCTGATTTAACTGCCCGCAAAGaa  <  1:2427932/36‑1 (MQ=255)
                                  |                                  
ATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGATTTTCTCTGATTTAACTGCCCGCAAAGAA  >  REL606/2267799‑2267867

Base quality scores:  ATCG  < 3 ≤  ATCG  < 21 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG