Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 305,640 G→T 18.6% L207F (TTG→TTT yahA → predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606305,6400G→T18.6% 4.8 ‑3.2 0.024621yahApredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (8/11):  total (13/11)
Fisher's exact test strand distribution p-value = 4.11e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGT  >  REL606/305605‑305675
                                   |                                   
cATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTg                                     <  1:1302881/36‑1 (MQ=255)
 atcttGGCATCAACGTCTCGGCGGGTTGTTTTTTgg                                    <  1:3208969/33‑1 (MQ=255)
 atatTGGCATCAACGGCTCGGCGGGTTGTTTTTTgt                                    >  1:1030526/1‑35 (MQ=255)
     tGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGc                                <  1:3097441/36‑1 (MQ=255)
     tGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGc                                <  1:4760583/36‑1 (MQ=255)
        cATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCggg                             <  1:5582976/36‑1 (MQ=255)
         aTCAACGTCTCGGCGGGTTGTTTTTTTGGAGGgggg                            >  1:2712913/1‑35 (MQ=25)
         aTCAACGTCTCGGCGGGTTGTTTTTTTGGAGCgggg                            >  1:4940801/1‑35 (MQ=37)
           cAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGAtt                          <  1:3171435/36‑1 (MQ=255)
            aaCGTCTCGGCGGGTTGTTTTTTGGCAGCGGGAttt                         <  1:5697694/36‑1 (MQ=255)
            aaCGTCTCGGCGGGGTGTTTTTTTGGAGCGGGTttt                         >  1:4056265/1‑36 (MQ=19)
               ggctcGGCGGGTTGTTTTCTGGCAGCGTGATTTGaa                      <  1:5446196/34‑1 (MQ=25)
               gTCTCGGCGGGTTTTTTTTTTGCAGGGGGGttttta                      >  1:3856804/1‑33 (MQ=12)
                tctcGGCGGGTTGTTTTTTGGCAGCGGGATTTGaaa                     <  1:244180/36‑1 (MQ=255)
                tctcGGCGGGTTGTTTTTTGGCAGCGGGATTTGaaa                     <  1:3144024/36‑1 (MQ=255)
                tctcGGCGGGTTGGTTTTTTGCAGCGGGATTTGaaa                     >  1:105054/1‑36 (MQ=37)
                      cggGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAgt               <  1:4126488/36‑1 (MQ=255)
                      cggGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAgt               >  1:3904038/1‑36 (MQ=255)
                       gggTTGTTTTTTGGCAGCGGGATTTGAAAAagagcg              >  1:57222/1‑34 (MQ=255)
                       gggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAgtg              >  1:4276297/1‑36 (MQ=255)
                       gggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAgtg              >  1:5477640/1‑36 (MQ=255)
                           tgttTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTg          <  1:4204898/36‑1 (MQ=255)
                             ttttttGGCAGCGGGATTTGAAAAAGAGTGTCTGaa        >  1:1829072/1‑36 (MQ=255)
                               ttttGGCAGCGGGATTTGAAAAAGAGTGTCTGAAcc      <  1:851016/36‑1 (MQ=255)
                                 ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTg    >  1:815479/1‑36 (MQ=255)
                                 ttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTg    >  1:984729/1‑36 (MQ=255)
                                  tGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTgg   >  1:4013770/1‑36 (MQ=255)
                                   ggCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGt  <  1:4992471/36‑1 (MQ=255)
                                   ggCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGt  >  1:5052943/1‑36 (MQ=255)
                                   |                                   
CATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGT  >  REL606/305605‑305675

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG