Integrated Genome Viewer "Motif Finder" like-tool for the command line?

Question

[ccc]

I recently asked about help finding promoters in particular, and got a helpful pointer for that... however, it wasn't exactly what I wanted I found out.

I then found something much closer to what I wanted: the Integrated Genome Viewer has a "Motif Finder" tool, and you can even enter in ambiguous nucleotides to it. The only catch (for me at least) is... is there a command line tool like it? Do I have to interact with IGV through the GUI? It seems there is some command line support for IGV, but it seems like it is geared more for pre-processing. Sorry if I missed something though, would appreciate any pointers!

Answer 1

[dariober]

If any useful, I wrote fastaRegexFinder for this.

[ccc]

This sounds exactly what I was looking for! Much appreciated :)

[danvoronov]

Does it work with IUPAC codes, i.e. does it find sequences with a G or T if I input K into regex? For instance, my fasta would be ATGGTACAATTGC, would using -r KKTAC return a GGTAC from the sequence? It would be useful since most genomic sequences only use ATGC while most transcription factor motifs use IUPAC ambiguity codes.