Motif presence at genomic location about basset HOT 4 OPEN

Yea this has been on my radar for awhile. I suppose I've hesitated because tomtom is so often unable to identify what proteins are most relevant to the motif. I'll post here if I code up a solution.

from basset.

Yeah, but I can definitely say that this is something that is absolutely crucial for this to catch real traction, need to connect to previous people's stuff. I started working on this and it should just be feeding the model the seqeunce of a specific location right?

from basset.

Just to make sure that we're on the same page- you want to annotate the motifs that pop up with loss scores in the saturated mutagenesis heat maps, right? For example, the sequence in Figure 5 in the paper would be annotated as CTCF. That's how I've been thinking about it.

Obviously, there are already tools for annotating sequences with motifs, like Tomtom. Basset can help suggest which motifs are actually relevant. So maybe the easiest approach would be to add an option to basset_sat.py to query the sequence for significant motifs with Tomtom, but filter the list for only those that overlap a nucleotide with a loss score above some threshold.

from basset.

Yeah that's exactly right (CTCF in figure 5).
And although the method you point out would show the identified ones, it would only show the identified ones? I think that there is definitely a TON of value in showing the Motifs that TomTom can't find. But that would mean starting w/ an all inclusive list of motifs which cause a significant loss score and just labeling those which TomTom is able to label, not discarding any (we know that current motifs suck).

from basset.