Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 330,478 C→G 26.3% A257P (GCC→CCC)  cynR ← DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606330,4780C→G26.3% 4.5 ‑4.5 0.022769cynRDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  new base (5/1):  ref base (3/13):  total (8/14)
Fisher's exact test strand distribution p-value = 1.09e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGG  >  REL606/330447‑330513
                               |                                   
aaaaCCGCCGTTCTCTCCAGTAGTGGCGGGGGAaga                                 >  1:1327204/1‑36 (MQ=255)
 aacccgccgTTCTCTCCAGTAGTGGCGGGGCAagag                                <  1:755274/33‑1 (MQ=255)
  caCCGCCGTTCGCTCCAGTAGTGGCGGGGCAAgaga                               <  1:730606/35‑1 (MQ=255)
  aaCCGCCGTTCTCTCCAGTAGTGGCGGGGGCAgaga                               >  1:5885622/1‑36 (MQ=37)
    ccgccgTTCTCTCCAGTAGTGGCGGGGGGAGAGgaa                             >  1:7116985/1‑36 (MQ=25)
        cgTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc                         <  1:1217390/36‑1 (MQ=255)
           tctctcCAGTAGTGGCGGGGCAAGAGAAATAGCttt                      <  1:6848098/36‑1 (MQ=255)
           tctctcCAGTAGTGGCGGGGCAAGAGAAATAGCttt                      <  1:496305/36‑1 (MQ=255)
            ctctcCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                     <  1:525115/36‑1 (MQ=255)
             tctcCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                    <  1:1202582/36‑1 (MQ=255)
               tcCAGTAGTGGCGGGGGGAGAGAAAAAACTTTaatt                  >  1:5936005/1‑34 (MQ=16)
                caagtggtGGCGGGGGAAGAGAAATAGCTTTAAGcc                 <  1:5459543/34‑1 (MQ=37)
                  agtagtGGCGGGGGAAGAGAAATAGATTTTAGCCCg               >  1:282881/1‑36 (MQ=37)
                       tGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg          <  1:3102018/36‑1 (MQ=255)
                        ggCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt         >  1:7172838/1‑36 (MQ=255)
                          cGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg       >  1:5296253/1‑36 (MQ=255)
                           ggggCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt      >  1:896117/1‑36 (MQ=255)
                           ggggCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt      <  1:5896781/36‑1 (MQ=255)
                           ggggCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt      <  1:1150146/36‑1 (MQ=255)
                            gggCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtg     <  1:1026322/36‑1 (MQ=255)
                             ggCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgt    <  1:2703953/36‑1 (MQ=255)
                             ggCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgt    <  1:834365/36‑1 (MQ=255)
                               cAAGAGAAATAGCTTTAAGCCCCTCATGTTGTGTgg  <  1:2189504/36‑1 (MQ=255)
                               |                                   
AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGG  >  REL606/330447‑330513

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