Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 306,033 G→T 20.6% L338F (TTG→TTT yahA → predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606306,0330G→T20.6% 2.5 ‑4.0 0.0351014yahApredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  new base (7/1):  ref base (10/17):  total (17/18)
Fisher's exact test strand distribution p-value = 1.78e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTATTTGTACTCTCCGCCAGTACCGGGTAA  >  REL606/306000‑306068
                                 |                                   
ggATTTAATGATCGGTCAATGCGTTCACTTTTTGCa                                   <  1:3389521/36‑1 (MQ=37)
  aTTTAATGATCGGTAAAGGCGTTCACTTTTTGCAgg                                 <  1:211713/36‑1 (MQ=255)
   tttAATGATCGGTAAAGGTGTTCACTTTTTGCAtgg                                >  1:1335955/1‑36 (MQ=37)
   tttAATGATCGGTAAAGGCGTTCACTTTTTTCAgtg                                >  1:882694/1‑34 (MQ=37)
   tttAATGATCGGTAAAGGCGTTCACTTTTTGCAggg                                <  1:421671/36‑1 (MQ=255)
    ttAATGATCGGTAAAGGCGTTCACTTTTTTCAgggg                               >  1:4135391/1‑35 (MQ=255)
    ttAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGc                               >  1:571281/1‑36 (MQ=255)
     tAATGATCGGTAAAGGCGTTCGCTTTTTGCAGGGCt                              <  1:4831399/36‑1 (MQ=255)
      aaTGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTa                             >  1:4536835/1‑36 (MQ=255)
        tGATCGGTAAAGGCGTTTACTTTTTTAAGGGCTAtt                           >  1:3272404/1‑36 (MQ=25)
         gATCGGTAAAGGCGTTCACTTTTTGCAGGGCTAttt                          <  1:843962/36‑1 (MQ=255)
          aTCGGTAAAGGCGTTCCCTTTTTTCAGGGCTCtttt                         >  1:2350305/1‑35 (MQ=21)
          aTCGGTAAAGGCGTTCACTTTTTGCAGGGCTATTTg                         <  1:2549207/36‑1 (MQ=255)
           tCGGTAAAGGCGTTCACTTTTTGCAGGGCTATTTGt                        <  1:2177242/36‑1 (MQ=255)
            ctgTAAAGCCGTTCCCTTTTTTCAGGGCTATTTGTa                       <  1:1410150/34‑1 (MQ=21)
            cGGTAAAGGCGTTCACTTTTTTCAGGGCTATTTTTa                       >  1:3406383/1‑36 (MQ=37)
             ggTAAAGGCGTTCACTTTTTTCTTGGCTATTTGTAc                      >  1:1615872/1‑36 (MQ=25)
              gTAAAGGCGTTCACTTTTTGCAGGGCTATTTTTAct                     >  1:3022754/1‑36 (MQ=255)
              gTAAAGGCGTTCACTTTTTGCAGGGCTATTTGTAct                     <  1:628262/36‑1 (MQ=255)
                aaaGGCGTTCACTTTTTGCAGGGCTATTTGTActct                   <  1:3629814/36‑1 (MQ=255)
                aaaGGCGTTCACTTTTTGCAGGGCTATTTGTActct                   <  1:3236007/36‑1 (MQ=255)
                 aaGGCGTTCACTTTTTGCAGGGCTATTTGTACtctc                  >  1:2831941/1‑36 (MQ=255)
                  aGGCGTTCACTTTTTTCAGGGCTATTTGTActctac                 >  1:5724830/1‑34 (MQ=38)
                  aGGCGTTCACTTTTTGCAGGGCTATTTGTACTCTcc                 <  1:2236733/36‑1 (MQ=255)
                    gCGTTCACTTTTTGCAGGGCTATTTGTACTCTCCGc               >  1:3758807/1‑36 (MQ=255)
                     cGTTCACTTTTTGCAGGGCTATTTGTACTCTCCGcc              <  1:4222126/36‑1 (MQ=255)
                      gTTCACTTTTTGCAGGGCTATTTGTACTCTCCGCCa             >  1:501088/1‑36 (MQ=255)
                       ttCACTTTTTGCAGGGCTATTTGTACTCTCCGCCAg            >  1:5466457/1‑36 (MQ=255)
                        tCACTTTTTGCAGGGCTATTTTTACTCTCCGCCAGt           >  1:3585792/1‑36 (MQ=255)
                          aCTTTTTGCAGGGCTATTTGTACTCTCCGCCAGTAc         <  1:5229084/36‑1 (MQ=255)
                          aCTTTTTGCAGGGCTATTTGTACTCTCCGCCAGTAc         >  1:1507925/1‑36 (MQ=255)
                           cTTTTTGCAGGGCTATTTGTACTCTCCGCCAGTAcc        <  1:5371917/36‑1 (MQ=255)
                           cTTTTTGCAGGGCTATTTGTACTCTCCGCCAGTAcc        <  1:2637290/36‑1 (MQ=255)
                             ttttGCAGGGCTATTTGTACTCTCCGCCAGTACCgg      <  1:3272743/36‑1 (MQ=255)
                               ttGCAGGGCTATTTGTACTCTCCGCCAGTACCGGGt    <  1:4053085/36‑1 (MQ=255)
                                 gCAGGGCTATTTGTACTCTCCGCCAGTACCGGGtaa  <  1:3452665/36‑1 (MQ=255)
                                 |                                   
GGATTTAATGATCGGTAAAGGCGTTCACTTTTTGCAGGGCTATTTGTACTCTCCGCCAGTACCGGGTAA  >  REL606/306000‑306068

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG