Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,285,943 G→T 15.6% A773E (GCA→GAA)  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,285,9430G→T15.6% 2.3 ‑4.2 ‑0.1322318arcBhybrid sensory histidine kinase in two‑component regulatory system with ArcA
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (10/17):  total (15/17)
Fisher's exact test strand distribution p-value = 1.49e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.49e-01

GGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGC  >  REL606/3285911‑3285978
                                |                                   
ggTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGcca                                  <  1:3238215/36‑1 (MQ=255)
  tCTAGCCGGGGTCATTTTTTAGTGGCTTTTGtcacc                                >  1:1241372/1‑36 (MQ=255)
  tCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccacc                                <  1:990461/36‑1 (MQ=255)
   cTAGCCGGGGTCATTTTTTAGTGGCTTTTTtcaccc                               >  1:2829206/1‑36 (MQ=37)
     aGCCGGAGTCATTTTTTAGTGGCTTTTGccacccac                             <  1:4440226/36‑1 (MQ=255)
      gCCGGGGTCATTTTTTAGTGGCTTTTTCCACCCACg                            >  1:2485524/1‑36 (MQ=255)
      gCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                            >  1:1929573/1‑36 (MQ=255)
        cGGGGTCATTTTTTAGTGTCTTTTGCCACCCACGCt                          >  1:2291174/1‑36 (MQ=255)
         ggggTCATTTTTTAGTGGCTTTTTCCACCCCCGCtt                         >  1:5743524/1‑36 (MQ=37)
          gggTCATTTTTTAGTGGCTTTTGCCACCCACGCttt                        <  1:3340404/36‑1 (MQ=255)
            gTCATTTTTTAGTGGCTTTTTCCACCCACGCTTTCa                      >  1:7279541/1‑36 (MQ=255)
             tCATTTTTTAGTGGCTTTTTCCACCCACGCTTTCAg                     >  1:7910456/1‑36 (MQ=255)
             tCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                     <  1:434227/36‑1 (MQ=255)
              cATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGc                    <  1:3679107/36‑1 (MQ=255)
                 tttttAGTGGCTTTTGCCACCCACGCTTTCAGCACt                 <  1:1515741/36‑1 (MQ=255)
                  ttttAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                <  1:5329531/36‑1 (MQ=255)
                   tttAGTGGCTTTTGCCACCCACGCTTTCAGCACTTc               <  1:3690205/36‑1 (MQ=255)
                    ttAGTGGCTTTTGCCACCCACGCTTTCCGCAATTCt              >  1:3521372/1‑36 (MQ=37)
                    ttAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCt              <  1:7885/36‑1 (MQ=255)
                    ttAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCt              <  1:4113721/36‑1 (MQ=255)
                     tAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTa             <  1:3138654/36‑1 (MQ=255)
                      aGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAc            <  1:7199610/36‑1 (MQ=255)
                       gTGGCTTTTGCCACCCATGCTTTGAGCACTTCTAcg           >  1:6165873/1‑36 (MQ=25)
                       gTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcg           >  1:5920289/1‑36 (MQ=255)
                       gTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcg           >  1:392448/1‑36 (MQ=255)
                         ggCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtc         <  1:3596761/36‑1 (MQ=255)
                          gCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg        <  1:1658629/36‑1 (MQ=255)
                           cTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcgt       <  1:911851/36‑1 (MQ=255)
                              ttGCCACCCACGCTTTCAGCACTTCTACGTCGTGAc    <  1:5105332/36‑1 (MQ=255)
                              ttGCCACCCACGCTTTCAGCACTTCTACGTCGTGAc    >  1:6995541/1‑36 (MQ=255)
                               tGCCACCCACGCTTTCAGCACTTCTACGTCGTGaag   >  1:1092205/1‑34 (MQ=255)
                               tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACg   >  1:1344396/1‑36 (MQ=255)
                                gCCACCCACGCTTTCACCACTTCTACGTCGTGACGc  >  1:6570829/1‑36 (MQ=255)
                                |                                   
GGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGC  >  REL606/3285911‑3285978

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG