Download short match track as bed file?

[Froberg, John]

Hello,

I am interested in downloading the output of the Short Match tool run for a short sequence against mm9 as a bedfile, which will provide a list of all instances of the motif of interest in the mm9 genome assembly. Is there any way to download the output of the Short Match tool as a bedfile? One previous answer to a similar question (https://groups.google.com/a/soe.ucsc.edu/forum/#!searchin/genome/short$20match/genome/85TponZ_RjM/7EJI3qg-Xc4J) suggested downloading the findMotif tool to obtain this information, and an early answer (https://groups.google.com/a/soe.ucsc.edu/forum/#!searchin/genome/short$20match/genome/rf7PJNcIjyg/ntq3-Je3SQoJ) implied that there is a way to download Short Match alignment coordinates from the browser, however that response links to a third response which appears to have been deleted. 

So has there been any progress obtaining Short Match alignment coordinates from the browser, without downloading and running findMotif?

Thank you very much,

John Froberg

Graduate Student, Lee Lab

Harvard University

[Brian Lee]

Dear John,

Thank you for using the UCSC Genome Browser and your question about downloading the short match track as a bed file.

The current path available will require running a command line utility, such as tacg or findMotif.

http://www.ncbi.nlm.nih.gov/pmc/articles/PMC99049/

http://tacg.sourceforge.net/

If you wish to use our findMotif utility, you can download it from here:

http://hgdownload.cse.ucsc.edu/admin/exe/

For example, if you are working in a macOSX.x86_64 environment:

Step1:   Download the appropriate utility for your environment

http://hgdownload.cse.ucsc.edu/admin/exe/macOSX.x86_64/findMotif

Step2: Make it executable 

chmod +x findMotif

Step3:  Run it on a 2bit file for a sequence, such as "aagaca":

./findMotif -motif=aagaca /path/to/a/sequenceSuchAS.GL339449.2bit

Get results in a bed format:

GL339449.1 938 944 1 1000 -

GL339449.1 2249 2255 2 1000 +

You can obtain the hg19.2bit file with the following command:

rsync -aP rsync://hgdownload.cse.ucsc.edu/goldenPath/hg19/bigZips/hg19.2bit .

You may want to read more about the findMotif utility in our mailing list archives.

https://groups.google.com/a/soe.ucsc.edu/forum/?hl=en&fromgroups#!search/findMotif

Thank you again for your inquiry and using the UCSC Genome Browser. If you have any further questions, please reply to gen...@soe.ucsc.edu. All messages sent to that address are archived on a publicly-accessible forum. If your question includes sensitive data, you may send it instead to genom...@soe.ucsc.edu.

All the best,

Brian Lee

UCSC Genome Bioinformatics Group

--