Hi, is there a way of loading a ucsc genome browser track automatically from a script (preferably in python but any language would do) so that a user invoking the script can directly see the visualization with the loaded track?
Thanks a lot in advance
Check this I just wrote a tutorial on how to do it. Of course, you can automate to any level. I also included how to it within R. http://www.biostars.org/post/show/42844/visualizing-chip-seq-data-using-ucsc-bigwig/
Doing it in python is not particularly difficult, you can use directly urrlib or alternatively use something like twill to automate web scripting.
At the same time as an alternative you might want to look at Track Hubs