custom Alt_PAS_AAUGAA 50 bases upstrand from poly(A) site

[CSV file]

Density of motif

t-test: p=6.3525e-06

Neutral set average density: 0.00266257

Different set average density: 0.00406439

Presence of motif

           motif_present
gene_set    absent present
  different    798     196
  neutral     1703     250

Fisher's exact test: p=1.2084e-06

Binomial GLM: p=8.64095e-07

Call:
glm(formula = motif_present ~ different, family = binomial, data = daf)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-0.6628  -0.6628  -0.5234  -0.5234   2.0276  

Coefficients:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept)   -1.91869    0.06773  -28.33  < 2e-16 ***
differentTRUE  0.51469    0.10461    4.92 8.64e-07 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 2505.1  on 2946  degrees of freedom
Residual deviance: 2481.4  on 2945  degrees of freedom
AIC: 2485.4

Number of Fisher Scoring iterations: 4