In [1]:
# Filepaths and Hard-coded Defaults

proj_root = "/users/kcochran/projects/procap_models/"
sequence_path = proj_root + "genomes/GRCh38_no_alt_analysis_set_GCA_000001405.15.fasta"
chrom_sizes = proj_root + "genomes/hg38.chrom.sizes.withrRNA"

in_window = 2114
out_window = 1000
In [2]:
# stuff to get from config file

with open("2022-07-03_11-35-25_run1_modisco_config_CACO2_profile.txt") as config_f:
    config_dict = {line.split()[0] : line.strip().split()[1] for line in config_f}

modisco_out_path = config_dict["modisco_out_path"]
scoring_type = config_dict["scoring_type"]
score_center_size = int(config_dict["score_center_size"])
profile_display_center_size = int(config_dict["profile_display_center_size"])

# digest what's in config file

assay_type, model_type, cell, accession, modisco_dir_base = modisco_out_path.split("/")[-5:]
ts_part1, ts_part2, run_str, _ = modisco_dir_base.split("_")
timestamp = ts_part1 + "_" + ts_part2
run = int(run_str.replace("run", ""))
In [3]:
print(modisco_out_path)
print("cell_type:", cell, accession)
print("timestamp:", timestamp)
print("run:", run)
print("scoring_type:", scoring_type)
print("score_center_size:", score_center_size)
print("profile_display_center_size:", profile_display_center_size)

/users/kcochran//projects/procap_models/modisco_out/procap/bpnetlite_basic_v2/CACO2/ENCSR100LIJ/2022-07-03_11-35-25_run1_modisco
cell_type: CACO2 ENCSR100LIJ
timestamp: 2022-07-03_11-35-25
run: 1
scoring_type: profile
score_center_size: 1000
profile_display_center_size: 400
In [4]:
data_dir = proj_root + "/data/procap/processed/" + cell + "/" + accession + "/"
plus_bw_path = data_dir + "final.5prime.pos.bigWig"
minus_bw_path = data_dir + "final.5prime.neg.bigWig"
val_peak_path = data_dir + "peaks_uni_and_bi_train_and_val.bed.gz"

val_save_dir = proj_root + "model_out/" + assay_type + "/" + model_type + "/" + cell + "/" + accession + "/"
val_save_path = val_save_dir + timestamp + "_run" + str(run) + "_train_and_val"

attr_save_path = val_save_dir.replace("model_out", "deepshap_out") + timestamp + "_run" + str(run) + "_deepshap"

if not modisco_out_path.endswith("/"):
    modisco_out_path = modisco_out_path + "/"
In [5]:
# task-specific filepaths

import os

assert scoring_type in ["profile", "counts"], scoring_type

if scoring_type == "profile":
    scores_path = attr_save_path + "_prof.npy"
    onehot_scores_path = attr_save_path + "_prof_onehot.npy"
    modisco_obj_path = modisco_out_path + "results_allChroms_prof_slice" + str(score_center_size) + ".hdf5"
    seqlet_path = modisco_out_path + "seqlets_prof.txt"
else:
    scores_path = attr_save_path + "_count.npy"
    onehot_scores_path = attr_save_path + "_count_onehot.npy"
    modisco_obj_path = modisco_out_path + "results_allChroms_count_slice" + str(score_center_size) + ".hdf5"
    seqlet_path = modisco_out_path + "seqlets_count.txt"
    
assert(os.path.exists(scores_path)), scores_path
assert(os.path.exists(onehot_scores_path)), onehot_scores_path
In [6]:
# Imports, Plotting Defaults

import matplotlib.pyplot as plt
import matplotlib.font_manager as font_manager

plot_params = {
    "figure.titlesize": 22,
    "axes.titlesize": 22,
    "axes.labelsize": 20,
    "legend.fontsize": 18,
    "xtick.labelsize": 16,
    "ytick.labelsize": 16,
    "font.weight": "bold"
}
plt.rcParams.update(plot_params)

from IPython.display import display
import tqdm
tqdm.tqdm_notebook()

import numpy as np

from view_modisco_results_utils import *

/users/kcochran/miniconda3/envs/procap/lib/python3.7/site-packages/ipykernel_launcher.py:19: TqdmDeprecationWarning: This function will be removed in tqdm==5.0.0
Please use `tqdm.notebook.tqdm` instead of `tqdm.tqdm_notebook`
In [7]:
# Load in True Profiles and Sequences

import sys
sys.path.append('../1_train_models')

from data_loading import extract_peaks

one_hot_seqs, true_profs = extract_peaks(sequence_path, 
    plus_bw_path, minus_bw_path, val_peak_path, in_window, out_window,
    max_jitter=0, verbose=True)

one_hot_seqs = one_hot_seqs.swapaxes(1,2)
one_hot_seqs = one_hot_seqs[:, (in_window // 2 - score_center_size // 2):(in_window // 2 + score_center_size // 2), :]

Reading FASTA: 100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 24/24 [00:15<00:00,  1.53it/s]
Loading Peaks: 42149it [01:23, 504.23it/s]
In [8]:
# Load in Coordinates of Examples
    
coords = load_coords(val_peak_path, in_window)
In [9]:
# Import SHAP scores, predicted profiles

hyp_scores = np.load(scores_path).swapaxes(1,2)
hyp_scores = hyp_scores[:, (in_window // 2 - score_center_size // 2):(in_window // 2 + score_center_size // 2), :]
pred_profs = np.exp(np.load(val_save_path + ".profs.npy"))
In [10]:
# Load modisco results object
    
tfm_obj = import_tfmodisco_results(modisco_obj_path, hyp_scores, one_hot_seqs)
In [11]:
motif_pfms, motif_hcwms, motif_cwms, \
    motif_pfms_short, num_seqlets, \
    motif_seqlets, num_metaclusters = plot_all_metaclusters(tfm_obj, one_hot_seqs, hyp_scores,
                                                            true_profs, pred_profs, coords,
                                                            in_window, out_window,
                                                            score_center_size,
                                                            profile_display_center_size)

Metacluster 1/1

Pattern 1/43

16401 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 2/43

5452 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 3/43

4969 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 4/43

4868 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 5/43

3588 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 6/43

3413 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 7/43

2273 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 8/43

2123 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 9/43

1993 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 10/43

1957 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 11/43

1669 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 12/43

1213 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 13/43

1149 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 14/43

1124 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 15/43

897 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 16/43

826 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 17/43

759 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 18/43

559 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 19/43

541 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 20/43

540 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 21/43

471 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 22/43

434 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 23/43

366 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 24/43

360 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 25/43

353 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 26/43

288 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 27/43

259 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 28/43

212 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 29/43

199 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 30/43

190 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 31/43

185 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 32/43

185 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 33/43

66 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 34/43

56 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 35/43

52 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 36/43

51 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 37/43

51 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 38/43

49 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 39/43

48 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 40/43

47 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 41/43

35 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 42/43

26 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)

Pattern 43/43

25 seqlets

Sequence (PFM)
Hypothetical contributions (hCWM)
Actual contributions (CWM)
In [12]:
run_and_plot_tomtom(modisco_out_path, motif_pfms, motif_hcwms, motif_pfms_short, num_metaclusters)

Metacluster 1/1

The output directory '/users/kcochran//projects/procap_models/modisco_out/procap/bpnetlite_basic_v2/CACO2/ENCSR100LIJ/2022-07-03_11-35-25_run1_modisco/tomtom' already exists.
Its contents will be overwritten.
Processing query 1 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 2 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.972149
Processing query 3 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.998346
Processing query 4 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 5 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.956394
Processing query 6 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 7 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.961768
Processing query 8 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.988758
Processing query 9 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.967197
Processing query 10 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.986014
Processing query 11 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 12 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 13 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.983959
Processing query 14 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.963844
Processing query 15 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.967376
Processing query 16 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.985243
Processing query 17 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.973244
Processing query 18 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.986604
Processing query 19 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 20 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 21 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.99552
Processing query 22 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 23 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.993576
Processing query 24 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.97421
Processing query 25 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.998602
Processing query 26 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.955021
Processing query 27 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 28 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.974727
Processing query 29 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.963201
Processing query 30 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.961092
Processing query 31 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.968775
Processing query 32 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.964181
Processing query 33 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 34 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 35 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.962674
Processing query 36 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.978711
Processing query 37 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 38 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 39 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.982215
Processing query 40 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.979317
Processing query 41 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.995125
Processing query 42 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=1
Processing query 43 out of 43 
Estimating pi_0 from all 3912 observed p-values.
Estimating pi_0.
Estimated pi_0=0.988154
Can't locate HtmlMonolithWr.pm:   /root/lib/perl/HtmlMonolithWr.pm: Permission denied at /software/meme/4.11.2/bin/tomtom_xml_to_html line 48.
Warning: tomtom_xml_to_html exited abnormally and may have failed to create HTML output.

Motif 1/43

No TOMTOM matches passing threshold

Motif 2/43

Motif IDq-valPWM
MA0742.2-KLF120.000158461
MA0516.3-SP20.000158461
MA0685.2-SP40.000158461
MA1511.2-KLF100.000158461
MA0741.1-KLF160.000163241

Motif 3/43

No TOMTOM matches passing threshold

Motif 4/43

Motif IDq-valPWM
MA0060.3-NFYA1.03745e-05
MA1644.1-NFYC1.03745e-05
MA0314.2-HAP31.03745e-05
MA0316.1-HAP50.00902935
MA0502.2-NFYB0.00902935

Motif 5/43

Motif IDq-valPWM
MA1475.1-CREB3L40.00109599
MA1438.1-atf-70.00208936
MA0609.2-CREM0.00208936
MA1348.1-TGA90.00332866
MA1346.1-TGA100.00332866

Motif 6/43

No TOMTOM matches passing threshold

Motif 7/43

Motif IDq-valPWM
MA1832.1-Zm00001d0023649.8738e-09
MA1817.1-Zm00001d0202673.84147e-08
MA1821.1-Zm00001d0205957.56045e-08
MA1833.1-Zm00001d0493647.56045e-08
MA1819.1-Zm00001d0058921.87445e-07

Motif 8/43

Motif IDq-valPWM
MA0506.2-Nrf10.0036152
MA1650.1-ZBTB140.0675921
MA1880.1-Hey0.0735238
MA1816.1-O110.168551
MA1833.1-Zm00001d0493640.168551

Motif 9/43

Motif IDq-valPWM
MA1832.1-Zm00001d0023642.7108e-05
MA1821.1-Zm00001d0205952.7108e-05
MA1818.1-Zm00001d0522297.99929e-05
MA1819.1-Zm00001d0058920.000692292
MA1833.1-Zm00001d0493640.000899951

Motif 10/43

Motif IDq-valPWM
MA0750.2-ZBTB7A1.02456e-05
MA0076.2-ELK43.30479e-05
MA0759.2-ELK30.000203723
MA0156.3-FEV0.000203723
MA0763.1-ETV30.000203723

Motif 11/43

Motif IDq-valPWM
MA1628.1-Zic1::Zic20.455417
MA1726.1-ZNF3310.455417
MA1722.1-ZSCAN310.455417

Motif 12/43

Motif IDq-valPWM
MA0808.1-TEAD30.00301317
MA1121.1-TEAD20.00903948
MA0809.2-TEAD40.0476454
MA0090.3-TEAD10.0476454
MA0243.1-sd0.077189

Motif 13/43

Motif IDq-valPWM
MA1890.1-Klf151.41023e-07
MA1892.1-Klf5-like1.03988e-05
MA1893.1-Klf6-7-like1.03988e-05
MA1961.1-PATZ10.00067257
MA1653.1-ZNF1480.00109176

Motif 14/43

Motif IDq-valPWM
MA1833.1-Zm00001d0493646.41262e-08
MA2022.1-LOB1.33291e-06
MA1817.1-Zm00001d0202674.35124e-06
MA1820.1-Zm00001d0243249.81509e-06
MA1819.1-Zm00001d0058929.81509e-06

Motif 15/43

Motif IDq-valPWM
MA0139.1-CTCF3.92434e-09
MA1102.2-CTCFL8.11769e-07
MA1929.1-CTCF9.87725e-07
MA1930.1-CTCF5.28165e-06
MA0531.1-CTCF7.38485e-05

Motif 16/43

Motif IDq-valPWM
MA0489.2-Jun0.0019904
MA1988.1-Atf30.0019904
MA1448.1-fos-10.0019904
MA0099.3-FOS::JUN0.00251172
MA1928.1-BNC20.00251172

Motif 17/43

Motif IDq-valPWM
MA0153.2-HNF1B4.57794e-07
MA0046.2-HNF1A4.80425e-07
MA1198.1-HAT20.122134
MA1212.1-ATHB-130.122134
MA0951.1-ATHB-160.18482

Motif 18/43

Motif IDq-valPWM
MA1893.1-Klf6-7-like0.103622
MA0742.2-KLF120.103622
MA0079.5-SP10.103622
MA1511.2-KLF100.103622
MA0740.2-KLF140.108232

Motif 19/43

Motif IDq-valPWM
MA0527.1-ZBTB332.15706e-13

Motif 20/43

Motif IDq-valPWM
MA1573.2-Thap115.31448e-08
MA0088.2-ZNF1430.0199755
MA1716.1-ZNF760.0285848
MA1625.1-Stat5b0.315818
MA0519.1-Stat5a::Stat5b0.406725

Motif 21/43

Motif IDq-valPWM
MA0114.4-HNF4A0.00161612
MA0484.2-HNF4G0.00161612
MA0677.1-Nr2f60.00161612
MA0856.1-RXRG0.00161612
MA1550.1-PPARD0.00161612

Motif 22/43

No TOMTOM matches passing threshold

Motif 23/43

Motif IDq-valPWM
MA0833.2-ATF46.43339e-05
MA1636.1-CEBPG8.9149e-05
MA0025.2-NFIL30.0321598
MA1702.1-Pdp10.0321598
MA0043.3-HLF0.0927352

Motif 24/43

Motif IDq-valPWM
MA1961.1-PATZ10.000694171
MA1723.1-PRDM90.000694171
MA1890.1-Klf150.000699883
MA1653.1-ZNF1480.000725633
MA1892.1-Klf5-like0.000725633

Motif 25/43

Motif IDq-valPWM
MA1278.1-DOF3.40.150889
MA1279.1-DOF1.50.150889
MA1274.1-DOF3.60.150889
MA1281.1-DOF5.10.150889
MA1268.1-CDF50.150889

Motif 26/43

Motif IDq-valPWM
MA1820.1-Zm00001d0243240.00155765
MA0146.2-Zfx0.00155765
MA1817.1-Zm00001d0202670.00228389
MA1819.1-Zm00001d0058920.00326366
MA0026.1-Eip74EF0.0033888

Motif 27/43

Motif IDq-valPWM
MA0379.1-MOT20.062833

Motif 28/43

Motif IDq-valPWM
MA0753.2-ZNF7400.0062538
MA1819.1-Zm00001d0058920.0062538
MA1833.1-Zm00001d0493640.0062538
MA1713.1-ZNF6100.0062538
MA0079.5-SP10.0062538

Motif 29/43

Motif IDq-valPWM
MA0506.2-Nrf10.000416793
MA1880.1-Hey0.00349234
MA1650.1-ZBTB140.0182802
MA1833.1-Zm00001d0493640.0204018
MA0146.2-Zfx0.0315681

Motif 30/43

Motif IDq-valPWM
MA0748.2-YY20.000316945
MA1927.1-YY1-20.000316945
MA1833.1-Zm00001d0493640.000316945
MA1819.1-Zm00001d0058920.000316945
MA1651.1-ZFP420.00033601

Motif 31/43

Motif IDq-valPWM
MA1833.1-Zm00001d0493640.000705261
MA1961.1-PATZ10.000705261
MA1513.1-KLF150.000740664
MA1713.1-ZNF6100.00107834
MA1653.1-ZNF1480.00377721

Motif 32/43

Motif IDq-valPWM
MA0846.1-FOXC20.00303233
MA0032.2-FOXC10.00303233
MA1861.1-FoxB0.00303233
MA0546.1-pha-40.00303233
MA1683.1-FOXA30.00549028

Motif 33/43

No TOMTOM matches passing threshold

Motif 34/43

No TOMTOM matches passing threshold

Motif 35/43

Motif IDq-valPWM
MA1820.1-Zm00001d0243240.000627126
MA1817.1-Zm00001d0202670.000727272
MA1819.1-Zm00001d0058920.000831335
MA1832.1-Zm00001d0023640.00229992
MA1821.1-Zm00001d0205950.00275167

Motif 36/43

Motif IDq-valPWM
MA1051.1-RAP2-30.00375025
MA1053.1-ERF1090.00562944
MA1754.1-ERF0730.0056309
MA0997.1-ERF0690.0108585
MA0975.1-CRF20.0108585

Motif 37/43

Motif IDq-valPWM
MA1573.2-Thap118.37313e-06

Motif 38/43

Motif IDq-valPWM
MA0484.2-HNF4G0.124764
MA0809.2-TEAD40.124764
MA0114.4-HNF4A0.124764
MA1111.1-NR2F20.331627
MA1121.1-TEAD20.331627

Motif 39/43

Motif IDq-valPWM
MA1859.1-FoxA-a0.0210462
MA0846.1-FOXC20.0210462
MA0284.2-CIN50.0210462
MA1867.1-FoxI-b0.022439
MA1606.1-Foxf10.0263175

Motif 40/43

Motif IDq-valPWM
MA1653.1-ZNF1480.0424802
MA1723.1-PRDM90.0424802
MA0471.2-E2F60.0492121
MA1630.2-ZNF2810.0711584
MA1961.1-PATZ10.0711584

Motif 41/43

Motif IDq-valPWM
MA1448.1-fos-10.0825551
MA0580.1-DYT10.339943
MA0841.1-NFE20.339943
MA0655.1-JDP20.339943
MA1633.2-BACH10.34033

Motif 42/43

No TOMTOM matches passing threshold

Motif 43/43

Motif IDq-valPWM
MA0528.2-ZNF2630.0155904
MA1965.1-SP50.0155904
MA1522.1-MAZ0.0155904
MA1284.1-TCP10.0155904
MA1627.1-Wt10.0155904