Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 904,176 +C 6.8% 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 6.8% 3.5 ‑3.6 0.044coding (249/369 nt)ECB_00850putative regulator of late gene expression
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (30/11):  total (30/14)
Fisher's exact test strand distribution p-value = 2.75e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TACACGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAGATGATGCGGAATGGACTATTACCAC  >  REL606/904141‑904211
                                    |                                   
tACACGGAAATGCCAGTGAAAGTCAGTGGCTTTaaa                                      <  1:2903234/36‑1 (MQ=255)
tACACGGAAATGCCAGTGAAAGTCAGGGGCTTTaaa                                      >  1:2242138/1‑36 (MQ=255)
 acacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.c                                    >  1:5527300/1‑36 (MQ=255)
 acacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.c                                    >  1:3891497/1‑36 (MQ=255)
 acacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.c                                    <  1:597402/36‑1 (MQ=255)
 acacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.c                                    >  1:409521/1‑36 (MQ=255)
 acacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.c                                    >  1:7409423/1‑36 (MQ=255)
 acacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.c                                    <  1:7175841/36‑1 (MQ=255)
  cacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.Ca                                   >  1:3768563/1‑36 (MQ=255)
  cacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.Ca                                   >  1:3522170/1‑36 (MQ=255)
  cacGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.Ca                                   >  1:2674939/1‑36 (MQ=255)
   acGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAg                                  >  1:2686663/1‑36 (MQ=255)
   acGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAg                                  >  1:922185/1‑36 (MQ=255)
    cGGAAATGCCAGTGAAAGTTAGTCGCTTTTAA.GAGc                                 >  1:4021069/1‑36 (MQ=19)
    cGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGc                                 >  1:641717/1‑36 (MQ=255)
    cGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGc                                 >  1:8129577/1‑36 (MQ=255)
    cGGAAAGGCCAGTGAAAGGCAGTGGCTTTAAA.CAGc                                 >  1:4473456/1‑36 (MQ=255)
     ggAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGcc                                >  1:7854083/1‑36 (MQ=255)
     ggAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGcc                                >  1:8090409/1‑36 (MQ=255)
       aaaTGCCAGTGAAAGTCAGTGGCTTTAAA.CAGCCGa                              <  1:6080189/36‑1 (MQ=255)
         aTGCCAGTGAAAGTCAGTGGCTTTAAA.CAGCAGATa                            >  1:2688592/1‑36 (MQ=255)
          tGCCAGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAg                           >  1:5765483/1‑36 (MQ=255)
           gCCAGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAga                          >  1:5821500/1‑36 (MQ=255)
           gCCAGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAga                          >  1:6158384/1‑36 (MQ=255)
            ccAGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAgat                         >  1:1394130/1‑36 (MQ=255)
            ccAGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAgat                         >  1:8133516/1‑36 (MQ=255)
             cAGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAgatg                        >  1:3094827/1‑36 (MQ=255)
             cAGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAgatg                        >  1:477891/1‑36 (MQ=255)
              aGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAgatga                       >  1:939369/1‑36 (MQ=255)
               aaaaaaGTAAGTTGCTTAAAC.CAGCCGATAgatgat                      <  1:1327789/33‑1 (MQ=16)
                 aaaaGACAGTGGCTTAAAA.CAGCCGATAGATGATGc                    <  1:1350662/35‑1 (MQ=37)
                  aaaGTCAGTGTTTTAAAC.CAGCCGATAGATGATGCg                   <  1:6869764/36‑1 (MQ=19)
                   aaGTCTGTCGCTAAAAC.CAGCCGATAGATGATGCgg                  <  1:4471076/36‑1 (MQ=16)
                   aaGTAAGTGGTTTAAAA.CAGCCGATAGATGATGCgg                  <  1:1803936/36‑1 (MQ=25)
                   aaGACAGTGGTTTAAAC.CAGCCGATAGATGATGCgg                  <  1:2846322/36‑1 (MQ=19)
                    aGTCAGTGGCTTTAAA.CAGCCGATAGATGATGCGGa                 >  1:8408416/1‑36 (MQ=255)
                    aGTCAGTGGCTTTAAA.CAGCCGATAGATGATGCGGa                 >  1:1924869/1‑36 (MQ=255)
                    aGTAAGGGGTTTAAAA.CAGCCGATAGATGAAGCGGa                 <  1:7538888/36‑1 (MQ=12)
                      aaaGTGGATTTAAA.CAGCCGATAGATGATGCGGAAt               <  1:8421005/34‑1 (MQ=38)
                        ggTGGCTAAAAC.CAGCCGATAGATGATGCGGAATgg             <  1:1987312/35‑1 (MQ=25)
                         gTGGCTTTAAA.CAGCCGATAGATGATGCGGAATGGa            >  1:1341691/1‑36 (MQ=255)
                         gTGGCTTTAAA.CAGCCGATAGATGATGCGGAATGGa            >  1:7517280/1‑36 (MQ=255)
                           ggCTTTAAA.CAGCCGATAGATGATGCGGAATGGACt          >  1:6412094/1‑36 (MQ=255)
                           ggCTTTAAA.CAGCCGATAGATGATGCGGAATGGACt          >  1:2554675/1‑36 (MQ=255)
                             ttttAAACCAGCCGATAGATGATGCGGAATGGACTa         <  1:6515169/35‑1 (MQ=37)
                             atttAAACCAGCCGATAGATGATGCGGAATGGACTa         <  1:3747890/35‑1 (MQ=37)
                              tttAAACCAGCCGATAGATGATGCGGAATGGACTAt        <  1:168952/36‑1 (MQ=37)
                              tttAAA.CAGCCGATAGATGATGCGGAATGGACTAtt       >  1:2483539/1‑36 (MQ=255)
                               ttAGA.CACCCTATAGATGATGCTGAATGGACTATTa      <  1:7313555/36‑1 (MQ=19)
                               ttAAA.CAGCCGATAGATGATGCGGAATGGACTATTa      >  1:8293549/1‑36 (MQ=255)
                               ttAAA.CAGCCGATAGATGATGCGGAATGGACTATTa      >  1:625141/1‑36 (MQ=255)
                               taaaa.CAGCCGATAGATGATGCGGAATGGACTATTa      <  1:3097564/34‑1 (MQ=255)
                                tAAA.AAGCCGATAGATGATGCGGAATGGACTATTac     >  1:7617853/1‑36 (MQ=255)
                                 aaa.CAGCCGATAGATGATGCGGAATGGACTATTacc    >  1:7647714/1‑36 (MQ=255)
                                   a.CAGCCGGTAGATGATGCGGAATGGACTATTgccac  <  1:446641/36‑1 (MQ=37)
                                    |                                   
TACACGGAAATGCCAGTGAAAGTCAGTGGCTTTAAA.CAGCCGATAGATGATGCGGAATGGACTATTACCAC  >  REL606/904141‑904211

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 18 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG