data_name = "K562"
modisco_root = "/srv/scratch/msharmin/mouse_hem/k562_modisco"
tomtom_report_root = "http://mitra.stanford.edu/kundaje/msharmin/report/tomtom_outs/cells"
task_dir = "fold0"
#perf_file = "/srv/scratch/msharmin/mouse_hem/with_tfd/full_mouse50/fineFactorized/task_273-naivegw/NaiveauPRC.txt"
#homer_root = "/srv/scratch/msharmin/mouse_hem/with_tfd/full_mouse50/Naive_scans"
#reportfile = "/mnt/lab_data/kundaje/msharmin/annotations/filtering samples_MS2.xlsx"
#sheetname = "filter23"

from matlas.matches import DenovoModisco
from matlas.genome_data import *
# from matlas.matches import prepare_pattern_matches
ob = DenovoModisco("{}/{}".format(modisco_root, task_dir))
# ob.fetch_tomtom_matches(save_report=True,
#                             tomtom_dir= "{0}/{1}_tomtomout".format("{}/{}".format(modisco_root, task_dir), 
#                                                                    DEFAULT_DATABASE))
ob.load_matched_motifs()
ob.get_motif_per_celltype(match_threshold=0.05)
#

ob.display_individual_table()

Pattern Name	TF Name(s)	Modisco
metacluster_1/pattern_0 # seqlets: 3789 SequenceContrib ScoresHyp_Contrib Scores	Fosl2, Fosl1, Bach2, Fos, Jun, Mafk, Nfe2l2, Nfe2, Atf3, Bach1, Maf, Jund, Fosb, Mafg, Junb, Mafb, Maff, Batf, Mafa, Jdp2
metacluster_1/pattern_1 # seqlets: 2778 SequenceContrib ScoresHyp_Contrib Scores	Sp2, Sp3, Sp1, Klf3, Klf6, Sp4, Klf15, Klf1, Maz, Wt1, Sp5, Zbtb17, Egr1, Zfx, Klf7, Klf4, Klf5, E2f4, Egr2, Klf8, Klf12, Zfp281, Rela, Usf2, Nr1h4, Sall4, Ptf1a, E2f1, Srebf2, Fli1, Tfap2c, Egr3, Mxi1, E2f7
metacluster_1/pattern_2 # seqlets: 1843 SequenceContrib ScoresHyp_Contrib Scores	Gata6, Gata4, Tal1, Gata2, Gata1, Gata3, Gata5
metacluster_1/pattern_3 # seqlets: 1409 SequenceContrib ScoresHyp_Contrib Scores	Atf2, Atf7, Creb1, Atf1, Creb5, Crem, Nfil3, Rxrb, Tef
metacluster_1/pattern_4 # seqlets: 1006 SequenceContrib ScoresHyp_Contrib Scores
metacluster_1/pattern_5 # seqlets: 956 SequenceContrib ScoresHyp_Contrib Scores	Elk1, Gabpa, Elf1, Erg, Ehf, Ets1, Elf4, Etv6, Etv4, Elk4, Spi1, Etv2, Elf2, Elf3, Spib, Elf5, Etv3, XP_911724.4, Elk3, Etv5, Spdef, Etv1, Fev, Irf8, Irf4, Spic, Stat4, Irf2, Prdm1
metacluster_1/pattern_6 # seqlets: 889 SequenceContrib ScoresHyp_Contrib Scores	Nfyb, Foxi1, Nfyc, Nfya, Pbx3
metacluster_1/pattern_7 # seqlets: 818 SequenceContrib ScoresHyp_Contrib Scores	Usf1, Tfeb, Tfe3, Arntl, Tfec, Mitf, Srebf1, Npas2, Clock, Arnt, Bhlhe40, Creb3l2, Mlx, Tcfl5, Mycn, Max, Xbp1, Myc, Bhlhe41, Hes1, Mlxip
metacluster_1/pattern_8 # seqlets: 729 SequenceContrib ScoresHyp_Contrib Scores	Nrf1
metacluster_1/pattern_9 # seqlets: 435 SequenceContrib ScoresHyp_Contrib Scores
metacluster_1/pattern_10 # seqlets: 358 SequenceContrib ScoresHyp_Contrib Scores	Ctcfl, Ctcf
metacluster_1/pattern_11 # seqlets: 305 SequenceContrib ScoresHyp_Contrib Scores	Rest
metacluster_1/pattern_12 # seqlets: 257 SequenceContrib ScoresHyp_Contrib Scores
metacluster_1/pattern_13 # seqlets: 257 SequenceContrib ScoresHyp_Contrib Scores	Zfp143, Thap11, Tbx2
metacluster_1/pattern_14 # seqlets: 237 SequenceContrib ScoresHyp_Contrib Scores
metacluster_1/pattern_15 # seqlets: 201 SequenceContrib ScoresHyp_Contrib Scores	Rfx2, Rfx3, Rfx1, Rfx4, Rfx7, Rfx6
metacluster_1/pattern_16 # seqlets: 173 SequenceContrib ScoresHyp_Contrib Scores	Yy1, Taf1
metacluster_1/pattern_17 # seqlets: 121 SequenceContrib ScoresHyp_Contrib Scores
metacluster_1/pattern_18 # seqlets: 82 SequenceContrib ScoresHyp_Contrib Scores	Runx1, Runx2, Runx3, Cbfb
metacluster_1/pattern_19 # seqlets: 77 SequenceContrib ScoresHyp_Contrib Scores	Mef2d, Mef2a, Mef2c
metacluster_1/pattern_20 # seqlets: 67 SequenceContrib ScoresHyp_Contrib Scores
metacluster_1/pattern_21 # seqlets: 59 SequenceContrib ScoresHyp_Contrib Scores
metacluster_1/pattern_22 # seqlets: 55 SequenceContrib ScoresHyp_Contrib Scores	Rorc, Nfatc1, Nr1d2