Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,141,000 C→A 32.5% Q25K (CAA→AAA)  yijC → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,141,0000C→A32.5% 5.2 ‑2.7 0.02473yijCDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (6/10):  total (6/18)
Fisher's exact test strand distribution p-value = 6.64e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC.AAAAAGAAAAAACCCGCCGTTCGCTGGTGGA  >  REL606/4140965‑4141034
                                   |                                   
gCAAGGATTCAGACATCGTGATGGGCGTAAGAgcgc                                     >  1:4219217/1‑36 (MQ=255)
gCAAGGATTCAGACATCGTGATGGGCGTAAGAgcgc                                     >  1:1262241/1‑36 (MQ=255)
  aaGGATTCAGACATCGTGATGGGCGTAAGAGCGgaa                                   >  1:14964/1‑36 (MQ=255)
  aaGGATTCAGACATCGTGATGGGCGTAAGAGCGcaa                                   >  1:954894/1‑36 (MQ=255)
  aaGGATTCAGACATCGTGATGGGCGTAAGAGCGcaa                                   <  1:2589514/36‑1 (MQ=255)
   aGGATTCAGACATCGTGATGGGCGTAAGAGCGcaac                                  <  1:3679470/36‑1 (MQ=255)
   aGGATTCAGACATCGTGATGGGCGTAAGAGCGcaac                                  >  1:1220731/1‑36 (MQ=255)
      aTTCAGACATCGTGATGTGCGTAAGAGCGCAAC.aaa                              >  1:1453889/1‑36 (MQ=255)
      aTTCAGACATCGTGATGGGCGTAAGAGCGCAAC.aaa                              <  1:1109742/36‑1 (MQ=255)
       ttCAGGCNNNGTGAGGGGCGGAAGAGCGCAAC.aaaa                             <  1:2021034/36‑1 (MQ=19)
       ttCAGACATCGTGATGGGCGTAAGAGCGCAAC.aaaa                             <  1:2788684/36‑1 (MQ=255)
       ttCAGACATCGTGATGGGCGTAAGAGCGCAAC.aaaa                             >  1:1785625/1‑36 (MQ=255)
        tCAGACATCGTGATGGGCGTAAGAGCGCAAC.aaaaa                            <  1:5917773/36‑1 (MQ=255)
           gACATCGTGATGGGCGTAAGAGCGCAAC.AAAAAGaa                         >  1:1600010/1‑36 (MQ=255)
             cATCGTGATGGGCGTAAGAGCGCAAC.AAAAAGaaaa                       <  1:1726711/36‑1 (MQ=255)
             cATCGTGATGGGCGTAAGAGCGCAAC.AAAAAGaaaa                       <  1:1518439/36‑1 (MQ=255)
             cATCGTGATGGGCGTAAGAGCGCAAC.AAAAAGaaaa                       >  1:1486267/1‑36 (MQ=255)
              aTCGTGATGGGCGTAAGAGCGCAAC.AANNAGaaaaa                      >  1:1276078/1‑36 (MQ=255)
              aTCGTGATGGGCGTAAGAGCGCAAC.AAAAAGaaaaa                      <  1:3444039/36‑1 (MQ=255)
               tCGTGATGGGAGTAAGAGCGCAAC.AAAAAGaaaaaa                     <  1:4372940/36‑1 (MQ=255)
                 gTGATGGGCGTAAGAGCGCAAC.AAAAAGAAAAAAcc                   <  1:3541218/36‑1 (MQ=255)
                    aTGGGCGTAAGAGCGCAAA.AAAAAGAAAAAACccgc                <  1:16741/36‑1 (MQ=255)
                    aTGGGCGTAAGAGCGAAAC.AAAAAGAAAAAACccgc                <  1:4045768/36‑1 (MQ=255)
                      gggCGTAAGAGAGAAAC.AAAAAAAAAAAACccgccg              <  1:4921432/36‑1 (MQ=25)
                      gggCGTAAAAG.GAAACAAAAAAGATAAAACccgccg              <  1:5257769/36‑1 (MQ=2)
                            aaaagCGAAAA.AAAAAAAAAAAACCCGCCGTTCGCt        <  1:361882/33‑1 (MQ=12)
                            aaGAGCGCAAC.AAAAAGAAAAAACCCGCCGTTCGCt        >  1:4897980/1‑36 (MQ=255)
                            aaGAGCGAAAC.AAAAAGAAAAAACCCGCCGTTCGCt        <  1:1454098/36‑1 (MQ=255)
                            aaGAGCGAAAC.AAAAAAAAAAAACCCGCCGTTCGCt        <  1:1897846/36‑1 (MQ=37)
                             agagCGAATA.AAAAAGAAAAAACCCGCCGTTCGCtg       <  1:104865/36‑1 (MQ=25)
                              gagCGAAAA.AAAAAGAAAAAACCCGCCGTTCGCtgg      <  1:2809434/36‑1 (MQ=37)
                                  gcAAC.AAAAAGAAAAAACCCGCCGTTCGCTGGTGGa  >  1:2657228/1‑36 (MQ=255)
                                   |                                   
GCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAC.AAAAAGAAAAAACCCGCCGTTCGCTGGTGGA  >  REL606/4140965‑4141034

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 14 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG