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% 155.3 46R57C (CGT→TGT) arcBhybrid 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:483503/1‑36 (MQ=255)
aCGGCGTAATCAGCAAACCAAAGAAGATAGAACAaa                                <  1:1138743/36‑1 (MQ=255)
  ggCGTAATCAGCAAACCAAAGAAGATAGAACAAATa                              >  1:1505084/1‑36 (MQ=255)
  ggCGTAATCAGCAAACCAAAGAAGATAGAACAAATa                              >  1:6002096/1‑36 (MQ=255)
     gTAATCAGCAAACCAAAGAAGATAGAACAAATAACa                           >  1:2098780/1‑36 (MQ=255)
        aTCAGCAAACCAAAGAAGATAGAACAAATAACATCa                        >  1:6035536/1‑36 (MQ=255)
        aTCAGCAAACCAAAGAAGATAGAACAAATAACATCa                        >  1:6314485/1‑36 (MQ=255)
         tCAGCAAACCAAAGAAGATAGAACAAATAACATCaa                       >  1:4410805/1‑36 (MQ=255)
         tCAGCAAACCAAAGAAGATAGAACAAATAACATCaa                       >  1:3206236/1‑36 (MQ=255)
         tCAGCAAACCAAAGAAGATAGAACAAATAACATCaa                       >  1:92353/1‑36 (MQ=255)
          cAGCAAACCAAAGAAGATAGAACAAATAACATCAAt                      >  1:6652764/1‑36 (MQ=255)
          cAGCAAACCAAAGAAGATAGAACAAATAACATCAAt                      <  1:3507936/36‑1 (MQ=255)
           aGCAAACCGAAGAAGCTAGAACAAATAACATCAATg                     >  1:5992727/1‑36 (MQ=21)
           aGCAAACCAAAGAAGATAGAACAAATAACATCAATg                     <  1:5088755/36‑1 (MQ=255)
            gCAAACCAAAGAAGATAGAACAAATAACATCAATGc                    <  1:4584927/36‑1 (MQ=255)
              aaaCCAAAGAAGATAGAACAAATAACATCAATGCtt                  <  1:2099848/36‑1 (MQ=255)
               aaCCAAAGAAGATAGAACAAATAACATCAATGCttt                 >  1:6866864/1‑36 (MQ=255)
                aCCAAAGAAGATAGAACAAATAACATCAATGCTTTc                >  1:3408056/1‑36 (MQ=255)
                 ccAAAGAAGATAGAACAAATAACATCAATGCTTTCg               <  1:1364436/36‑1 (MQ=255)
                  cAAAGAAGATAGAACAAATAACATCAATGCTTTCGa              >  1:1900865/1‑36 (MQ=255)
                  cAAAGAAGATAGAACAAATAACATCAATGCTTTCGa              >  1:6653240/1‑36 (MQ=255)
                  cAAAGAAGATAGAACAAATAACATCAATGCTTTCGa              >  1:3640128/1‑36 (MQ=255)
                  cAAAGAAGATAGAACAAATAACATCAATGCTATcaa              >  1:1208116/1‑34 (MQ=21)
                   aaaGAAGATAGAACAAATAACATCAATGCTTTCGAc             >  1:2755693/1‑36 (MQ=255)
                    aagaagATAGAACAAATAACATCAATGCTTTCGAcc            >  1:85792/1‑36 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           >  1:928301/1‑36 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           <  1:4457757/36‑1 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           >  1:358566/1‑36 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           >  1:7012416/1‑36 (MQ=255)
                     agaagaTAGAACAAATAACATCAATGCTTTCGACCt           >  1:1864204/1‑36 (MQ=255)
                       aagaTAGAACAAATAACATCAATGCTTTCGACCTGa         <  1:4228681/36‑1 (MQ=255)
                        agaTTGAACAAATAACATCAATGCTTTCGACCTGAc        <  1:4048474/36‑1 (MQ=25)
                        agaTAGAACAAATAACATCAATGCTTTCGACCTGAc        >  1:536369/1‑36 (MQ=255)
                        agaTAGAACAAATAACATCAATGCTTTCGACCTGAc        >  1:6534279/1‑36 (MQ=255)
                         gaTAGAACAAATAACATCAATGCTTTCGACCTGAcc       >  1:1583155/1‑36 (MQ=255)
                         gaTAGAACAAATAACATCAATGCTTTCGACCTGAcc       >  1:6645481/1‑36 (MQ=255)
                          aTAGAACAAATAACATCAATGCTTTCGACCTGACCa      >  1:2873623/1‑36 (MQ=255)
                          aTAGAACAAATAACATCAATGCTTTCGACCTGACCa      >  1:4745090/1‑36 (MQ=255)
                           tAGAACAAATAACATCAATGCTTTCGACCTGACCAt     <  1:5878860/36‑1 (MQ=255)
                            aGAACAAATAACATCAATGCTTTCGACCTGACCATg    <  1:6024926/36‑1 (MQ=255)
                            aGAACAAATAACATCAATGCTTTCGACCTGACCATg    >  1:5771856/1‑36 (MQ=255)
                             gAACAAATAACATCAATGCTTTCGACCTGACCATgc   >  1:1607104/1‑36 (MQ=255)
                             gAACAAATAACATCAATGCTTTCGACCTGACCATgc   <  1:6884015/36‑1 (MQ=255)
                             gAACAAATAACATCAATGCTTTCGACCTGACCATgc   <  1:7297675/36‑1 (MQ=255)
                              aaCAAATAACATCAATGCTTTCGACCTGACCATgca  >  1:6591419/1‑36 (MQ=255)
                              aaCAAATAACATCAATGCTTTCGACCTGACCATgca  <  1:7464872/36‑1 (MQ=255)
                                 |                                
ACGGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCA  >  REL606/3288059‑3288124

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