Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 904,176 +C 9.1% coding (249/369 nt) ECB_00850 → putative regulator of late gene expression

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606904,1761.→C 9.1% 7.7 ‑2.9 ‑0.033coding (249/369 nt)ECB_00850putative regulator of late gene expression
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (20/10):  total (20/13)
Fisher's exact test strand distribution p-value = 5.24e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.08e-01

TACACGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAGATGATGCGGAATGGACTATTACCAC  >  REL606/904141‑904211
                                    |                                   
tACACGGAAATGCCAGTGAAAGTCAGTGGCTTTaaa                                      <  1:3866211/36‑1 (MQ=255)
tACACGGAAATGCCAGTGAAAGTCAGTGGCTTTaaa                                      <  1:7222973/36‑1 (MQ=255)
tACACGGAAATGCCAGTGAAAGTCAGTGGCTTTaaa                                      >  1:3536014/1‑36 (MQ=255)
 acacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.c                                    <  1:1928277/36‑1 (MQ=255)
 acacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.c                                    <  1:6246112/36‑1 (MQ=255)
  cacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.Ca                                   <  1:7452014/36‑1 (MQ=255)
  cacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.Ca                                   <  1:6018067/36‑1 (MQ=255)
  cacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.Ca                                   <  1:6021003/36‑1 (MQ=255)
   acGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAg                                  <  1:4564316/36‑1 (MQ=255)
    cGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGc                                 <  1:1457950/36‑1 (MQ=255)
     ggAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGcc                                <  1:2876409/36‑1 (MQ=255)
     ggAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGcc                                >  1:4925524/1‑36 (MQ=255)
     ggAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGcc                                >  1:4839699/1‑36 (MQ=255)
     ggAAATGCCAATGAATCATAGTGGCTTTAAA.CAGcc                                <  1:388191/36‑1 (MQ=16)
      aaaaTGCCAGTGAAAGTCAGTGGCTTTAAA.CAGCCg                               <  1:2840482/35‑1 (MQ=255)
        aaTGCCAGTGAAAGTCAGTGGCTTTAAA.CAGCCGAt                             >  1:4236152/1‑36 (MQ=255)
        aaTGCCAGTGAAAGTCAGAGGATATAAA.CAGCCGag                             >  1:6293817/1‑35 (MQ=25)
        aaTGCCAGTGAAAGCCAGCGGCTTTAAA.CAGCAGAt                             >  1:4585991/1‑36 (MQ=21)
                tGAAAGTCAGTGGCTTTAAA.CAGCCGATAGatgatg                     >  1:2820758/1‑36 (MQ=255)
                tGAAAGGCAGTGGCTTTAAA.CAGACGATAGatgatg                     >  1:2489157/1‑36 (MQ=38)
                tGAAAGACAGTGGCTTTAAA.CAAACGATAGaggatg                     >  1:7695712/1‑36 (MQ=19)
                  aaaGTCAGTGGCTTTAAA.CAGCCGATAGATGATGCg                   >  1:6379192/1‑36 (MQ=255)
                   aaGTCAGTGGCTTTAAA.CAGCCGATAGATGATGCgg                  >  1:5391545/1‑36 (MQ=255)
                    aGTCAGTGGCTTTAAA.CAGCCGATAGATGATGCGGa                 >  1:6065521/1‑36 (MQ=255)
                      taaGTGGTTTAAAA.CAGCCGATAGATGATGCGGAAt               <  1:6461248/34‑1 (MQ=25)
                       cAGTGGCTTTAAA.CAGCCGATAGATGATGCGGAATg              >  1:7238670/1‑36 (MQ=255)
                       cAGTGGCTTTAAA.CAGCCGATAGATGATGCGGAATg              >  1:3642126/1‑36 (MQ=255)
                         gTGGTTTTAAC.CAGCCGATAGATGATGCGGAATGGa            <  1:195054/36‑1 (MQ=37)
                         gTGGCTTTAAA.CAGCCGATAGATGATGCGGAATGGa            >  1:4431909/1‑36 (MQ=255)
                           gggtttAAACCAGCCGATAGATGATGCGGAATGGAc           <  1:7574324/33‑1 (MQ=25)
                           ggCTTTAAA.CAGCCGATAGATGATGCGGAATGGACt          >  1:463371/1‑36 (MQ=255)
                           ggCTTTAAA.CAGCCGATAGATGATGCGGAATGGACt          >  1:3287702/1‑36 (MQ=255)
                            ggattAAACCAGCCGATAGATGATGCGGAATGGACt          <  1:2619035/33‑1 (MQ=25)
                            gCTTTAAA.CAGCCGATAGATGATGCGGAATGGACTa         >  1:2096016/1‑36 (MQ=255)
                            gCTTTAAA.CAGCCGATAGATGATGCGGAATGGACTa         >  1:2404257/1‑36 (MQ=255)
                              tttAAACCAGCCGATAGATGATGCGGAATGGACTAt        <  1:6548168/36‑1 (MQ=37)
                               ttAAA.CAGCCGATAGATGATGCGGTATGGACTATTa      >  1:425200/1‑36 (MQ=255)
                               ttAAA.CAGCCGATAGATGATGCGGAATGGACTATTa      >  1:7528369/1‑36 (MQ=255)
                                 aaa.CAGCCGATAGATGATGCGGAATGGACTATTacc    >  1:6702430/1‑36 (MQ=255)
                                 aaa.CAGCCGATAGATGATGCGGAATGGACTATTacc    >  1:5556136/1‑36 (MQ=255)
                                  aa.CAGCCGATAGATGATGCGGAATGGACTATTacca   >  1:6298142/1‑36 (MQ=255)
                                  aa.CAGCCGATAGATGATGCGGAATGGACTATTacca   >  1:1904090/1‑36 (MQ=255)
                                   a.CAGCCGATAGATGATGCGGAATGGACTATTaccac  >  1:5895496/1‑36 (MQ=255)
                                   a.CAGCCGATAGATGATGCGGAATGGACTATTaccac  >  1:4910920/1‑36 (MQ=255)
                                   a.CAGCCGATAGATGATGCGGAATGGACTATTaccac  >  1:1903860/1‑36 (MQ=255)
                                    |                                   
TACACGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAGATGATGCGGAATGGACTATTACCAC  >  REL606/904141‑904211

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG