Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,157,725 C→A 13.0% P759Q (CCA→CAA)  ebgA → cryptic beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,157,7250C→A13.0% 4.8 ‑4.5 ‑0.0512276ebgAcryptic beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  new base (0/7):  ref base (24/20):  total (24/27)
Fisher's exact test strand distribution p-value = 1.07e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.75e-01

GGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACC  >  REL606/3157689‑3157760
                                    |                                   
ggCAGGTGAATGGCG.ATTGCTTCTTACTCGCGGGcc                                     >  1:1605177/1‑36 (MQ=9)
  cAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCa                                    <  1:2481884/36‑1 (MQ=255)
  cAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCa                                    <  1:2795640/36‑1 (MQ=255)
   aGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCaa                                   >  1:4496561/1‑36 (MQ=255)
    ggTGAATGGCGAATCGCTGCTGACTCGCGAGCCaaa                                  <  1:637020/36‑1 (MQ=255)
    ggTGAATGGCGAATCGCTGCTGACTCGCGAGCCaaa                                  <  1:4821933/36‑1 (MQ=255)
    ggTGAATGGCGAATCGCTGCTGACTCGCGAGCCaaa                                  >  1:1939930/1‑36 (MQ=255)
    ggTGAATGGCGAATCGCTGCTGACTCGCGAGCCaaa                                  <  1:7615663/36‑1 (MQ=255)
     gTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAg                                 >  1:5027352/1‑36 (MQ=255)
     gTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAg                                 >  1:3883333/1‑36 (MQ=255)
      tGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGa                                <  1:6802429/36‑1 (MQ=255)
      tGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGa                                >  1:3264558/1‑36 (MQ=255)
      tGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGa                                >  1:2138718/1‑36 (MQ=255)
      tGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGa                                <  1:3095192/36‑1 (MQ=255)
        aaTGGCGAATCGCTGCTGACTCGCGAGCCAAAGagc                              >  1:8414669/1‑34 (MQ=255)
        aaTGGCGAATCGCTGCTGACTCGCGAGCCAAAGATc                              >  1:7827904/1‑36 (MQ=255)
            gCGAATCGCTGCTGACTCGCGAGCCAAAGATCAAct                          <  1:7402892/36‑1 (MQ=255)
             cGAATCGCTGCTGACTCGCGAGCCAAAGATCAActt                         <  1:5789617/36‑1 (MQ=255)
              gAATCGCTGCTGACTCGCGAGCCAAAGATCAActtc                        <  1:6837683/36‑1 (MQ=255)
              gAATCGCTGCTGACTCGCGAGCCAAAGATCAActtc                        <  1:4171610/36‑1 (MQ=255)
              gAATCGCTGCTGACTCGCGAGCCAAAGATCAActtc                        <  1:1682896/36‑1 (MQ=255)
               aaTCGCTGCTGACTCGCGAGCCAAAGATCAActtct                       >  1:593872/1‑36 (MQ=255)
               aaTCGCTGCTGACTCGCGAGCCAAAGATCAActtct                       <  1:3974825/36‑1 (MQ=255)
               aaTCGCTGCTGACTCGCGAGCCAAAGATCAActtct                       >  1:833775/1‑36 (MQ=255)
                aTCGCTGCTGACTCGCGAGCCAAAGATCAActtctt                      <  1:7993250/36‑1 (MQ=255)
                  cGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCa                    <  1:1713536/36‑1 (MQ=255)
                   gctgctGACTCGCGAGCCAAAGATCAACTTCTTCaa                   >  1:5345431/1‑36 (MQ=255)
                    ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg                  >  1:8214121/1‑36 (MQ=255)
                    ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg                  <  1:6716415/36‑1 (MQ=255)
                    ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg                  >  1:6787966/1‑36 (MQ=255)
                    ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg                  >  1:809999/1‑36 (MQ=255)
                    ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg                  <  1:8059951/36‑1 (MQ=255)
                    ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg                  <  1:3020884/36‑1 (MQ=255)
                      gctgACTCGCGAGCCAAAGATCAACTTCTTCAAGcc                >  1:1277112/1‑36 (MQ=255)
                       ctgACTCGCGAGCAAAAGATCAACTTCTTCAAGCCg               <  1:1304428/36‑1 (MQ=255)
                        tgACTCGCGAGCCAAAGATCAACTTCTTCAAGCCga              >  1:5055852/1‑36 (MQ=255)
                        tgACTCGCGAGCCAAAGATCAACTTCTTCAAGCCga              >  1:4779987/1‑36 (MQ=255)
                        tgAATCGCGAGCCAAAGATCAACTTCTTCAAGccta              <  1:968614/36‑3 (MQ=38)
                          aCTCGCGAGCCAAAGATCAACTTCTTCAAGCAgatg            >  1:3986893/1‑36 (MQ=255)
                          aCTCGCGAGAAAAAGATCAACTTCTTCAAGCCgatg            <  1:3339472/36‑1 (MQ=37)
                            tCGCGAGCCAAAGATCAACTTCTTCAAGCCgatgat          >  1:2321322/1‑36 (MQ=255)
                             cgcgAGCCAAACATAAACTTCTTCAAGCCGATGATc         <  1:7269077/36‑1 (MQ=37)
                             cgcgAGCAAAAGATCAACTTCTTCAAGCCGATGATc         <  1:1942250/36‑1 (MQ=255)
                             cgcgAGAAAAAGATCAACTTCTTCAAGCCGATGATc         <  1:8264184/36‑1 (MQ=37)
                             cgcgAGAAAAAGATAAACTTCTTCAAGCCGATGATc         <  1:3833658/36‑1 (MQ=25)
                               cgAGCCAAAGATCAACTTCTTCAAGCCGATGATCGa       >  1:7221201/1‑36 (MQ=255)
                               cgAGCAAAAGATCAACTTCTTCAAGCCGATGATCGa       <  1:721836/36‑1 (MQ=255)
                                gAGCCAAAGATCAACTTCTTCAATCCGATGATAGAc      >  1:8183398/1‑36 (MQ=37)
                                gAGCCAAAGATCAACTTCTTCAAGCCGATGATCGAc      >  1:3487193/1‑36 (MQ=255)
                                 aGCCAAAGATCAACTTCTTCAAGCCGATGATCGACa     >  1:740766/1‑36 (MQ=255)
                                 aGCCAAAGATCAACTTCTTCAAGCCGATGATCGACa     >  1:5958385/1‑36 (MQ=255)
                                 aGCAAAAGATCAACTTCTACAAGCCGATGATCGACa     <  1:3156827/36‑1 (MQ=37)
                                  gacAAAGATCAACTTCTTAAAGCCGATGATCGCCaa    <  1:6864121/34‑1 (MQ=37)
                                    cAAAGATCAACTTCTTCAAGCCGATGATCGACAAcc  >  1:1823399/1‑36 (MQ=255)
                                    cAAAGATCAACTTCTTCAAGCCGATGATCGACAAcc  >  1:2577235/1‑36 (MQ=255)
                                    cAAAGATCAACTTCTTCAAGCCGATGATCGACAAcc  >  1:1653345/1‑36 (MQ=255)
                                    |                                   
GGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACC  >  REL606/3157689‑3157760

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG