click on the Biowiki logo to go to homepage



Research Teaching | Blog
Main | GBrowse | TWiki
Biowiki > GBrowse > DotplotTracks

Search

Advanced search...

Topics


Links

PageRank Checker

How about dotplot tracks for the genome browser?

No need to go through the TiledImage layer for this, either: we could generate the dotplot tiles independently, direct from sequence.

Simple tracks could just show dotplots a la dotlet. Or dotter (Erik Sonnhamer's cool X utility with nice, responsive contrast widget.... sadly don't think javascript is up to this!). More advanced tracks could show posterior probabilities for local direct repeats, inverted repeats and stem-loop structures, using various grammars.

-- IanHolmes - 28 Feb 2006

The program dottrack in dart/src/gfx now creates dotplot tracks, conveniently broken into tiles.

Type cd dart/src/gfx; make dottrack

Binary is placed in dart/bin. Type dottrack -h for help.

Basic usage is dottrack sequencefile.fasta but try also -log 5 for progress logging.

Currently does not do less than one base per pixel, i.e. super zoomed-in view (but can do more than one base/pixel, i.e. zoomed-out).

Should be pretty fast I think smile

For reference, at one-base-per-pixel, pixel intensity at (x,y) is M(S_{x-y/2}, S_{x+y/2}) where M(a,b)=\delta(a,b) for "forward track" and M(a,b)=\delta(\bar{a},b) (where \bar{a} is complement of a) for "reverse track".

Some contrast issues. Hard to spot repeats. Should be averaging over diagonals, not square neighborhoods. Fix this!

-- IanHolmes - 11 Mar 2006

Actions: Edit | Attach | New | Ref-By | Printable view | Raw view | Normal view | See diffs | Help | More...