Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,077,161 C→A 13.6% C364F (TGC→TTC)  hybB ← predicted hydrogenase 2 cytochrome b type component

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,077,1610C→A13.6% 7.0 ‑4.0 ‑0.0571091hybBpredicted hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  new base (1/7):  ref base (30/19):  total (31/26)
Fisher's exact test strand distribution p-value = 1.81e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.94e-01

GTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAA  >  REL606/3077126‑3077195
                                   |                                  
gTATCGGCAGTAGACGAATTAGTACGATGTAAgcgc                                    <  1:4867153/36‑1 (MQ=255)
gTATCGGCAGTAGACGAAAGAGTATGAAGTAAacgc                                    >  1:1354323/1‑36 (MQ=19)
 tATCGGCAGTAGACGAATGAGTACGATGTAAGCGCa                                   <  1:7194199/36‑1 (MQ=255)
 tATCGGCAGTAGACGAATGAGTACGATGTAAGCGCa                                   <  1:6168212/36‑1 (MQ=255)
 tATCGGCAGTAGACGAATGAGTACGATGTAAGCGCa                                   <  1:4435499/36‑1 (MQ=255)
 tATCGGCAGTAGACGAATGAGTACGATGTAAGCGCa                                   <  1:7118964/36‑1 (MQ=255)
 tATCGACAGTAGACGAATGAGTACGATGTAAGCGCa                                   <  1:7867679/36‑1 (MQ=255)
    cGGCAGTAGACGAATGAGTACGATGTAAGCGCAAAt                                <  1:607187/36‑1 (MQ=255)
    cGGCAGTAGACGAATGAGTACGATGTAAGCGCAAAt                                >  1:4108887/1‑36 (MQ=255)
     ggCAGTAGGCGAATGAGTACGATGTAAGCGCAAAtc                               <  1:6144067/36‑1 (MQ=255)
     ggCAGTAGACGAATGAGTACGATGTAAGCGCAAAtc                               <  1:5462379/36‑1 (MQ=255)
     ggCAGTAGACGAATGAGTACGATGTAAGCGCAAAtc                               <  1:2721082/36‑1 (MQ=255)
     ggCAGTAGACGAATGAGCACGATGTCAGCGAAAAtc                               >  1:4399474/1‑36 (MQ=25)
      gCAGTAGACGAATGAGTACGATGTAAGCGCAAAtct                              >  1:4771377/1‑36 (MQ=255)
       cAGTAGACGAATGAGTACGATGTAAGCGCAAAtctc                             <  1:1687953/36‑1 (MQ=255)
       cAGTAGACGAATGAGTACGATGTAAGCGCAAAtctc                             >  1:1894543/1‑36 (MQ=255)
       cAGTAGACGAATGAGTACGATGTAAGCGCAAAtctc                             >  1:6693189/1‑36 (MQ=255)
       cAGTAGACGAATGAGTACGATGTAAGCGCAAAtctc                             >  1:2838397/1‑36 (MQ=255)
        aGTAGACGAATGAGTACGATGTAAGCGCAAATCTCa                            >  1:5434842/1‑36 (MQ=255)
        aGTAGACGAATGAGTACGATGTAAGCGCAAATCTCa                            <  1:5371067/36‑1 (MQ=255)
         gTAGACGAATGAGTACGATGTAAGCGCAAATCTCaa                           >  1:6826455/1‑36 (MQ=255)
         gTAGACGAATGAGTACGATGTAAGCGCAAATCTCaa                           <  1:704344/36‑1 (MQ=255)
         gTAGACGAATGAGTACGATGTAAGCGCAAATCTCaa                           >  1:5086686/1‑36 (MQ=255)
         gTAGACGAATGAGTACGATGTAAGCGCAAATCTCaa                           <  1:3580272/36‑1 (MQ=255)
          tAGACGAATGAGTACGATGTAAGCGCAAATCTCAAt                          <  1:6637301/36‑1 (MQ=255)
           ataCGAATGAGTACGATGTAAGCGCAAATCTCAATa                         >  1:6035722/3‑36 (MQ=255)
           aGACGAATGATTACGATGTAAGCGCAAATCTCAATa                         <  1:140590/36‑1 (MQ=255)
            gACGAATGAGTACGATGTAAGCGCAAATCTCAATAg                        >  1:1106815/1‑36 (MQ=255)
            gACGAATGAGTACGATGTAAGCGCAAATCTCAATAg                        >  1:4082804/1‑36 (MQ=255)
            gACGAATGAGTACGATGTAAGCGCAAATCTCAATAg                        >  1:3872518/1‑36 (MQ=255)
             aCGAATGAGTACGATGTAAGCGCAAATCTCAATAGc                       >  1:7150771/1‑36 (MQ=255)
             aCGAATGAGTACGATGTAAGCGCAAATCTCAATAGc                       >  1:4429809/1‑36 (MQ=255)
              cGAATGAGTACGATGTAAGCGCAAATCTCAATAGcc                      >  1:6431944/1‑36 (MQ=255)
              cGAATGAGTACGATGTAAGCGCAAATCTCAATAGcc                      <  1:2799294/36‑1 (MQ=255)
               gAATGAGTACGATGTAAGCGCAAATCTCAATAGCca                     <  1:7786804/36‑1 (MQ=255)
               gAATGAGTACGATGTAAGCGCAAATCTCAATAGCca                     <  1:5403699/36‑1 (MQ=255)
                aaTGAGTACGATGTAAGCGCAAATCTCAATAGCcac                    >  1:4700097/1‑36 (MQ=255)
                 aTGAGTACGATGGAAGCGCAAATCTCAATAGCcaca                   >  1:5010373/1‑36 (MQ=255)
                  tGAGTACGATGTAAGCGCAAATCTCAATAGCCACaa                  >  1:6423207/1‑36 (MQ=255)
                  tGAGTACGATGTAAGCGCAAATCTCAATAGCCACaa                  >  1:2135075/1‑36 (MQ=255)
                  tGAGTACGATGTAAGCGCAAATCTCAATAGCCACaa                  >  1:6775519/1‑36 (MQ=255)
                   gAGTAAGATGTAAGCGCAAATCTCAATAGCCACaaa                 <  1:8198487/36‑1 (MQ=255)
                   gAGCACGATGTAAGCGCAAATATTAACAGCCACaaa                 >  1:48974/1‑36 (MQ=19)
                       acaaTGCAAGCGAAAATCAAAATAGCCACAAAACCa             <  1:4300993/33‑1 (MQ=12)
                       aagATGTAAGCAAAAATCTAAATAGCCACAAAACCa             <  1:8368222/34‑1 (MQ=25)
                       aCGATGTAAGCGAAAATCTAAATAGCCACAAAACCa             <  1:3341537/36‑1 (MQ=37)
                       aCGATGAAAGCGAAAATCTAAATAGCCACAAAACCa             <  1:7701739/36‑1 (MQ=25)
                       aCGATAAAAGCCAAAGTCTAAATAGCCACAAAACCa             <  1:4344097/36‑1 (MQ=12)
                        cGATGTAAGCGCAAATCTCAATAGCCACAAAACCaa            >  1:5747126/1‑36 (MQ=255)
                        cGATGTAAGCGCAAATCTCAATAGCCACAAAACCaa            >  1:486090/1‑36 (MQ=255)
                         gATGAAAGCAAAAATCTAAATAGCCACAAAACCAAt           <  1:1338019/36‑1 (MQ=16)
                          aTGTAAGCGCAAATCTCAATAGCCACAAAACCAATa          >  1:8072477/1‑36 (MQ=255)
                          aTGTAAGCGCAAATCTCAATAGCCACAAAACCAATa          >  1:758364/1‑36 (MQ=255)
                           tGTAAGCGCAAATCTCAATAGCCACAAAACCAATAg         >  1:6017626/1‑36 (MQ=255)
                           tGTAAGAGAAAATCTAAATAGCCACAAAACCAATAg         <  1:1561554/36‑1 (MQ=21)
                            gTAAGCGCAAATCTCAATAGCCACAAAACCAATAGa        >  1:6829845/1‑36 (MQ=255)
                             tAAGCGCAAATCTCAATAGCCACAAAACCAATAGaa       >  1:3758576/1‑36 (MQ=255)
                              aaGCGCAAATCTCAATAGCCACAAAACCAATAGaaa      >  1:4248201/1‑36 (MQ=255)
                              aaGCGCAAATCTCAATAGCCACAAAACCAATAGaaa      >  1:6981142/1‑36 (MQ=255)
                                gcgcAAATCTCAATAGCCACAAAACCAATAGAAATc    >  1:6755818/1‑36 (MQ=255)
                                  gcAAATCTCAATAGCCACAAAACCAATAGAAATCaa  >  1:7343851/1‑36 (MQ=255)
                                  gcAAATCTCAATAGCCACAAAACCAATAGAAATCaa  >  1:1596757/1‑36 (MQ=255)
                                  gcAAATCTCAATAGCCACAAAACCAATAGAAATCaa  >  1:6501160/1‑36 (MQ=255)
                                  gcAAATCTCAATAGCCACAAAACCAATAGAAATCaa  >  1:6656745/1‑36 (MQ=255)
                                  gcAAATCTCAATAGCCACAAAACCAATAGAAATCaa  >  1:6761956/1‑36 (MQ=255)
                                   |                                  
GTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAA  >  REL606/3077126‑3077195

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG