Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,141,009 G→A 23.6% E28K (GAA→AAA)  yijC → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,141,0090G→A23.6% 15.1 ‑0.7 ‑0.05382yijCDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  new base (7/6):  ref base (27/13):  total (34/19)
Fisher's exact test strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.88e-01

CAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTT  >  REL606/4140974‑4141044
                                   |                                   
cAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg                                     >  1:7802241/1‑36 (MQ=255)
 aGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa                                    >  1:5572699/1‑36 (MQ=255)
  gACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaa                                   >  1:2816648/1‑36 (MQ=255)
    cATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaa                                 <  1:7454934/36‑1 (MQ=255)
     aTCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaa                                <  1:5579115/36‑1 (MQ=255)
        gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACc                             <  1:4787648‑M2/36‑2 (MQ=255)
        gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACc                             >  1:7341501‑M2/1‑35 (MQ=255)
        gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACc                             >  1:1767827‑M2/1‑35 (MQ=255)
         tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACcc                            <  1:656620‑M2/36‑3 (MQ=255)
         tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACcc                            <  1:2923902‑M2/36‑3 (MQ=255)
         tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACcc                            <  1:6852849‑M2/36‑3 (MQ=255)
            tGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgcc                         <  1:5571892‑M2/36‑6 (MQ=255)
            tGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgcc                         >  1:7433392‑M2/1‑31 (MQ=255)
                cGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTc                     >  1:819558/1‑36 (MQ=255)
                agTAAGAGCAAAACAAAAAGAAAAAACCCGCCGTTc                     <  1:7841874/35‑1 (MQ=25)
                 gTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCg                    >  1:3447303/1‑36 (MQ=255)
                   aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCt                  >  1:1930025/1‑36 (MQ=255)
                   aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCt                  >  1:8334920/1‑36 (MQ=255)
                    agagCGAAAAAAAAAGAAAAACCCCGCCGTTCGCtg                 <  1:6090784/36‑1 (MQ=25)
                     gagCGCAACAAAAAGAAAAAACCCGCCGTTCGCtgg                >  1:4956375/1‑36 (MQ=255)
                      agCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggt               >  1:6638661/1‑36 (MQ=255)
                      agCGAAAAAAAAAAAAAAAACCCGCCGTTCGCtggt               <  1:5434703/36‑1 (MQ=18)
                       gcgcAACAAAAAGAAAAAACCCGCCGTTCGCtggtg              >  1:621075/1‑36 (MQ=255)
                        cgcAACAAAAAGAAAAAACCCGCCGTTCGCtggtgg             >  1:809768/1‑36 (MQ=255)
                        cgcAACAAAAAGAAAAAACCCGCCGTTCGCtggtgg             >  1:2747246/1‑36 (MQ=255)
                         gaaacaaAAAGAAAAACCCCGCCGTTCGCTGGTGGa            <  1:873195/34‑1 (MQ=38)
                          caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGaa           >  1:8032097/1‑36 (MQ=255)
                           aaaaaaaaGAAAAAACCCGCCGTTCGCTGGTGGAAg          <  1:1817160/33‑1 (MQ=255)
                            acaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc         >  1:4757174/1‑36 (MQ=255)
                            acaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc         >  1:2777313/1‑36 (MQ=255)
                            acaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc         >  1:3377767/1‑36 (MQ=255)
                            aaaaaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGc         <  1:849823/34‑1 (MQ=255)
                            aaaaaaaAAAAAAACCCGCCGTTCGCTGGTGGAAGc         <  1:8007084/34‑1 (MQ=37)
                            aaaaaaaAAAAAAACCCGCCGTTCGCTGGTGGAAGc         <  1:5901898/34‑1 (MQ=37)
                             caaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGcc        >  1:5843321/1‑36 (MQ=255)
                             caaAAAGAAAAAACCCGCCGCTCGCTGGTGGAAGcc        >  1:8301022/1‑36 (MQ=255)
                             caaAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGcc        <  1:5913146/36‑1 (MQ=255)
                             aaaaaaGAAAAACCCCGCCGTTCGCTGGTGGAAGcc        <  1:3153254/35‑1 (MQ=255)
                             aaaaaaAAAAAAACCCGCCGTTCGCTGGTGGAAGcc        <  1:6589166/35‑1 (MQ=255)
                              aaaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCg       >  1:4411566/1‑36 (MQ=255)
                              aaaaaAAAAAAACccgccagtccacgtaataacgtt       >  1:6205237‑M2/1‑13 (MQ=255)
                              aaaaaAAAAAAACccgccagtccacgtaataacgtt       >  1:3203353‑M2/1‑13 (MQ=255)
                              aaaaaAAAAAAACCCGCCGTTCGCTGGTGGAAGCCg       <  1:2574093/36‑1 (MQ=255)
                               gaaaaAAAAAACccgccgaaaactgtaccgagaaca      >  1:1713578‑M2/2‑12 (MQ=255)
                               aaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGc      >  1:4498743/1‑36 (MQ=255)
                               aaaaGAAAAAACCCGACGTTCGCTGGTGGAAGCCGc      >  1:7363172/1‑36 (MQ=255)
                               aaaaAAAAAAACccgccgtaccagattaaggtattc      >  1:2278046‑M2/1‑12 (MQ=255)
                               aaaaAAAAAAACccgccgaaaactgtaccgagaaca      >  1:2798616‑M2/1‑12 (MQ=255)
                               aaaaAAAAAAACccgccagtccacgtaataacgttt      >  1:1408672‑M2/1‑12 (MQ=255)
                               aaaaAAAAAAACccgccagtccacgtaataacgttt      >  1:3633073‑M2/1‑12 (MQ=255)
                                aaaGAAAAAACCCGCCGTGCGCTGGTGGAAGCCGCa     >  1:2614074/1‑36 (MQ=255)
                                aaaAAAAAAACccgccgaaaactgtaccgagaacaa     >  1:4701667‑M2/1‑11 (MQ=255)
                                 aaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt    >  1:5964286/1‑36 (MQ=255)
                                 aaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt    >  1:4585948/1‑36 (MQ=255)
                                  ggaAAGAACccgcctctgcgctggctcatgcggtac   >  1:2980193‑M2/2‑9 (MQ=255)
                                  ggaAAAAAAccgcctctgcgctggctcatgcggtac   >  1:1686106‑M2/2‑9 (MQ=255)
                                  aGAAAAAACccgcctctgcgctggctcatgcggtac   >  1:3421879‑M2/1‑9 (MQ=255)
                                   gAAAAAACccgcctctgcgctggctcatgcggtaca  >  1:2043722‑M2/1‑8 (MQ=255)
                                   gAAAAAACCCGCCGTTCGCTGGGGGAAGCCGCAttt  >  1:3239008/1‑36 (MQ=255)
                                   |                                   
CAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTT  >  REL606/4140974‑4141044

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG