from matlas.matches import DenovoAitac
motif_name = 'REST'
aitacdir = "/mnt/lab_data2/msharmin/oc-atlas/DanSkinData/fold_a_alex/{}".format(motif_name)
ob = DenovoAitac(aitacdir)
# ob.fetch_tomtom_matches(
# meme_db="/mnt/lab_data/kundaje/users/msharmin/annotations/HOCOMOCOv11_core_pwms_HUMAN_mono.renamed.nonredundant.annotated.meme",
# database_name="HOCOMOCO.nonredundant.annotated",
# save_report=True, tomtom_dir= "{0}/{1}_tomtomout".format(aitacdir, "HOCOMOCO.nonredundant.annotated"))
ob.load_matched_motifs(database_name="HOCOMOCO.nonredundant.annotated")
ob.get_motif_per_celltype(match_threshold=0.05, database_name="HOCOMOCO.nonredundant.annotated")
pattern_tab, pattern_dict = ob.visualize_pattern_table('Aitac')
tf_tab, tf_dict = ob.visualize_tf_table("Aitac")
from vdom.helpers import (b, summary, details)
from IPython.display import display
display(details(summary('Click here for ', b('Denovo Patterns'), ' by ', b('{}'.format('Aitac')),
' in ', b(motif_name),
": #{}".format(len(pattern_dict)),
), pattern_tab))
display(details(summary('Click here for ', b('Motifs'), ' by ', b('{}'.format('Aitac')),
' in ', b(motif_name),
": #{}".format(len(tf_dict)),
), tf_tab))