Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,288,092 G→A 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→A100.0% 156.0 46169arcBhybrid sensory histidine kinase in two‑component regulatory system with ArcA
Reads supporting (aligned to +/- strand):  new base (32/14):  ref base (0/0):  total (32/14)

ACGGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCA  >  REL606/3288059‑3288124
                                 |                                
aCGGCGTAATCAGCAAACCAAAGAAGATAGAACAaa                                <  1:1138743/36‑1 (MQ=255)
aCGGCGTAATCAGCAAACCAAAGAAGATAGAACAaa                                >  1:483503/1‑36 (MQ=255)
  ggCGTAATCAGCAAACCAAAGAAGATAGAACAAATa                              >  1:1505084/1‑36 (MQ=255)
  ggCGTAATCAGCAAACCAAAGAAGATAGAACAAATa                              >  2:2262886/1‑36 (MQ=255)
     gTAATCAGCAAACCAAAGAAGATAGAACAAATAACa                           >  1:2098780/1‑36 (MQ=255)
        aTCAGCAAACCAAAGAAGATAGAACAAATAACATCa                        >  2:2575275/1‑36 (MQ=255)
        aTCAGCAAACCAAAGAAGATAGAACAAATAACATCa                        >  2:2296326/1‑36 (MQ=255)
         tCAGCAAACCAAAGAAGATAGAACAAATAACATCaa                       >  2:671595/1‑36 (MQ=255)
         tCAGCAAACCAAAGAAGATAGAACAAATAACATCaa                       >  1:92353/1‑36 (MQ=255)
         tCAGCAAACCAAAGAAGATAGAACAAATAACATCaa                       >  1:3206236/1‑36 (MQ=255)
          cAGCAAACCAAAGAAGATAGAACAAATAACATCAAt                      >  2:2913554/1‑36 (MQ=255)
          cAGCAAACCAAAGAAGATAGAACAAATAACATCAAt                      <  1:3507936/36‑1 (MQ=255)
           aGCAAACCGAAGAAGCTAGAACAAATAACATCAATg                     >  2:2253517/1‑36 (MQ=25)
           aGCAAACCAAAGAAGATAGAACAAATAACATCAATg                     <  2:1349545/36‑1 (MQ=255)
            gCAAACCAAAGAAGATAGAACAAATAACATCAATGc                    <  2:845717/36‑1 (MQ=255)
              aaaCCAAAGAAGATAGAACAAATAACATCAATGCtt                  <  1:2099848/36‑1 (MQ=255)
               aaCCAAAGAAGATAGAACAAATAACATCAATGCttt                 >  2:3127654/1‑36 (MQ=255)
                aCCAAAGAAGATAGAACAAATAACATCAATGCTTTc                >  1:3408056/1‑36 (MQ=255)
                 ccAAAGAAGATAGAACAAATAACATCAATGCTTTCg               <  1:1364436/36‑1 (MQ=255)
                  cAAAGAAGATAGAACAAATAACATCAATGCTTTCGa              >  1:3640128/1‑36 (MQ=255)
                  cAAAGAAGATAGAACAAATAACATCAATGCTTTCGa              >  1:1900865/1‑36 (MQ=255)
                  cAAAGAAGATAGAACAAATAACATCAATGCTTTCGa              >  2:2914030/1‑36 (MQ=255)
                  cAAAGAAGATAGAACAAATAACATCAATGCTATcaa              >  1:1208116/1‑34 (MQ=25)
                   aaaGAAGATAGAACAAATAACATCAATGCTTTCGAc             >  1:2755693/1‑36 (MQ=255)
                    aagaagATAGAACAAATAACATCAATGCTTTCGAcc            >  1:85792/1‑36 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           <  2:718547/36‑1 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           >  1:928301/1‑36 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           >  2:3273206/1‑36 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           >  1:1864204/1‑36 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           >  1:358566/1‑36 (MQ=255)
                       aagaTAGAACAAATAACATCAATGCTTTCGACCTGa         <  2:489471/36‑1 (MQ=255)
                        agaTTGAACAAATAACATCAATGCTTTCGACCTGAc        <  2:309264/36‑1 (MQ=37)
                        agaTAGAACAAATAACATCAATGCTTTCGACCTGAc        >  2:2795069/1‑36 (MQ=255)
                        agaTAGAACAAATAACATCAATGCTTTCGACCTGAc        >  1:536369/1‑36 (MQ=255)
                         gaTAGAACAAATAACATCAATGCTTTCGACCTGAcc       >  2:2906271/1‑36 (MQ=255)
                         gaTAGAACAAATAACATCAATGCTTTCGACCTGAcc       >  1:1583155/1‑36 (MQ=255)
                          aTAGAACAAATAACATCAATGCTTTCGACCTGACCa      >  2:1005880/1‑36 (MQ=255)
                          aTAGAACAAATAACATCAATGCTTTCGACCTGACCa      >  1:2873623/1‑36 (MQ=255)
                           tAGAACAAATAACATCAATGCTTTCGACCTGACCAt     <  2:2139650/36‑1 (MQ=255)
                            aGAACAAATAACATCAATGCTTTCGACCTGACCATg    <  2:2285716/36‑1 (MQ=255)
                            aGAACAAATAACATCAATGCTTTCGACCTGACCATg    >  2:2032646/1‑36 (MQ=255)
                             gAACAAATAACATCAATGCTTTCGACCTGACCATgc   >  1:1607104/1‑36 (MQ=255)
                             gAACAAATAACATCAATGCTTTCGACCTGACCATgc   <  2:3144805/36‑1 (MQ=255)
                             gAACAAATAACATCAATGCTTTCGACCTGACCATgc   <  2:3558465/36‑1 (MQ=255)
                              aaCAAATAACATCAATGCTTTCGACCTGACCATgca  >  2:2852209/1‑36 (MQ=255)
                              aaCAAATAACATCAATGCTTTCGACCTGACCATgca  <  2:3725662/36‑1 (MQ=255)
                                 |                                
ACGGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCA  >  REL606/3288059‑3288124

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 25 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG