output should be sorted
default threshold for scanning
gimme genome: retrieve genome + annotation from UCSC and create index