Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,288,092 G→A 91.1% R57C (CGT→TGT)  arcB ← hybrid sensory histidine kinase in two‑component regulatory system with ArcA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,288,0920G→A91.1% 5.3 ‑0.4 ‑0.256169arcBhybrid sensory histidine kinase in two‑component regulatory system with ArcA
Reads supporting (aligned to +/- strand):  new base (33/18):  ref base (4/1):  total (37/19)
Fisher's exact test strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.85e-01

ACGGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCA  >  REL606/3288059‑3288124
                                 |                                
aCGGCGTAATCAGCAAACCAAAGAAGATAGAACGaa                                >  1:6483560/1‑36 (MQ=255)
aCGGCGTAATCAGCAAACCAAAGAAGATAGAACAaa                                >  1:1419120/1‑36 (MQ=255)
aCGGCGTAATAAGCAAACCAAAGAAGATAGAACAaa                                <  1:1451968/36‑1 (MQ=37)
 cGGCGTAATCAGCAAACCAAAGAAGATAGAACAAAt                               >  1:8296480/1‑36 (MQ=255)
  ggCGTAATCAGCAAACCAAAGAAGATAGAACAAATa                              <  1:7203877/36‑1 (MQ=255)
     gTAATCAGCAAACCAAAGAAGATAGAACAAATAACa                           <  1:5365357/36‑1 (MQ=255)
     gTAATCAGCAAACCAAAGAAGATAGAACAAATAACa                           >  1:7424133/1‑36 (MQ=255)
      tAATCAGCAAACCAAAGAAGATAGAACAAATAACAt                          >  1:5565318/1‑36 (MQ=255)
      tAATCAGCAAACCAAAGAAGATAGAACAAATAACAt                          >  1:8457003/1‑36 (MQ=255)
       caTCAGAAAACAAAAGAAGATAGAACAAATAACATc                         <  1:3078609/35‑1 (MQ=25)
       aaTCAGCAAACCAAAGAAGATAGAACAAATAACATc                         >  1:5031259/1‑36 (MQ=255)
        aTCAGCAAACAAAAGAAGATAGAACAAATAACATCa                        <  1:631758/36‑1 (MQ=37)
        aTCAGCAAAACAAAGAAGATAGAACAAATAACATCa                        <  1:5653472/36‑1 (MQ=37)
         tCAGCAAACAAAAGAAGATAGAACAAATAACATCaa                       <  1:6158576/36‑1 (MQ=37)
          cAGAAAACCAAAGAAGATAGAACAAATAACATCAAt                      <  1:6645534/36‑1 (MQ=37)
          aaGCAAAACAAAGAAGATAGAACAAATAACATCAAt                      <  1:7908326/35‑1 (MQ=37)
           aGCAAACCAAAGAAGATAGAACAAATAACATCAATg                     >  1:6216736/1‑36 (MQ=255)
              aaaCCAAAGAAGATAGAACGAATAACATCAATGCtt                  >  1:7771129/1‑36 (MQ=255)
              aaaCCAAAGAAGATAGAACGAATAACATCAATGCtt                  >  1:3802920/1‑36 (MQ=255)
              aaaCCAAAGAAGATAGAACAAATAACATCAATGCtt                  >  1:2007635/1‑36 (MQ=255)
              aaaCCAAAGAAGATAGAACAAATAACATCAATGCtt                  >  1:8076397/1‑36 (MQ=255)
              aaaCCAAAGAAGATAGAACAAATAACATCAATGCtt                  >  1:3883526/1‑36 (MQ=255)
              aaaACAAAGAAGATAGAACAAATAACATCAATGCtt                  <  1:2963251/36‑1 (MQ=37)
               aaCCAAAGAAGATAGAACAAATAACATCAATGCttt                 >  1:1360775/1‑36 (MQ=255)
               aaCCAAAGAAGATAGAACAAATAACATCAATGCttt                 >  1:2347393/1‑36 (MQ=255)
               aaCCAAAGAAGATAGAACAAATAACATCAATGCttt                 >  1:1066318/1‑36 (MQ=255)
                aCCAAAGAAGATAGAACGAATAACATCAATGCTTTc                <  1:518693/36‑1 (MQ=255)
                aCCAAAGAAGATAGAACAAATAACATCAATGCTTTc                >  1:5689871/1‑36 (MQ=255)
                 ccAAAGAAGATAGAACAAATAACATCAATGCTTTCg               >  1:2306469/1‑36 (MQ=255)
                  cAAAGAAGATAGAACAAATAACATCAATGCTTTCGa              <  1:5133860/36‑1 (MQ=255)
                   aaaGAAGATAGAACAAATAACATCAATGCTTTCGAc             >  1:8413026/1‑36 (MQ=255)
                   aaaGAAGATAGAACAAATAACATCAATGCTTTCGAc             >  1:2692649/1‑36 (MQ=255)
                   aaaGAAGATAGAACAAATAACATCAATGCTTTCGAc             >  1:514067/1‑36 (MQ=255)
                   aaaGAAGATAGAACAAATAACATCAATGCTTTCGAc             >  1:1147349/1‑36 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           >  1:6567312/1‑36 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           >  1:7568554/1‑36 (MQ=255)
                      gaagaTAGAACAAATAACATCAATGCTTTCGACCTg          >  1:7256122/1‑36 (MQ=255)
                       aagaTAGAACAAATAACATCAATGCTTTCGACCTGa         >  1:5703301/1‑36 (MQ=255)
                       aagaTAGAACAAATAACATCAATGCTTTCGACCTGa         >  1:3478911/1‑36 (MQ=255)
                       aagaTAGAACAAATAACATCAATGCTTTCGACCTGa         <  1:222227/36‑1 (MQ=255)
                       aagaTAGAACAAATAACATCAATGCTTTCGACCTGa         <  1:8131288/36‑1 (MQ=255)
                       aagaTAGAACAAATAACATCAATGCTTTCGACCTGa         <  1:8120404/36‑1 (MQ=255)
                        agaTAGAACAAATAACATCAATGCTTTCGACCTGAc        >  1:830932/1‑36 (MQ=255)
                        agaTAGAACAAATAACATCAATGCTTTCGACCTGAc        >  1:1491855/1‑36 (MQ=255)
                         gaTAGAACGAATAACATCAATGCTTTCGACCTGAcc       >  1:2634583/1‑36 (MQ=255)
                          aTAGAACAAATAACATCAATGCTTTCGACCTGACCa      >  1:2091877/1‑36 (MQ=255)
                          aTAGAACAAATAACATCAATGCTTTCGACCTGACCa      >  1:2699835/1‑36 (MQ=255)
                           tAGAACAAATAACATCAATGCTTTCGACCTGACCAt     >  1:1803653/1‑36 (MQ=255)
                           tAGAACAAATAACATCAATGCTTTCGACCTGACCAt     >  1:7118853/1‑36 (MQ=255)
                            aGAACAAATAACATCAATGCTTTCGACCTGACCATg    <  1:7617548/36‑1 (MQ=255)
                             gAACAAATAACATCAATGCTTTCGACCTGACCATgc   <  1:1619357/36‑1 (MQ=255)
                             gAACAAATAACATCAATGCTTTCGACCTGACCATgc   <  1:7342867/36‑1 (MQ=255)
                             gAACAAATAACATCAATGCTTTCGACCTGACCATgc   >  1:697768/1‑36 (MQ=255)
                             gAACAAATAACATCAATGCTTTCGACCTGACCATgc   >  1:2803789/1‑36 (MQ=255)
                             gAACAAATAACATCAATGCTTTCGACCTGACCATgc   <  1:852580/36‑1 (MQ=255)
                              aaCAAATAACATCAATGCTTTCGACCTGACCATgca  >  1:926412/1‑36 (MQ=255)
                                 |                                
ACGGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCA  >  REL606/3288059‑3288124

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