neolitec

package com.neolitec.examples;

import org.apache.commons.codec.binary.Base64;
import org.apache.commons.lang.StringUtils;
import org.slf4j.Logger;
import org.slf4j.LoggerFactory;

import javax.servlet.*;
import javax.servlet.http.HttpServletRequest;
import javax.servlet.http.HttpServletResponse;

brentp

UPSTREAM=400
INSTREAM=100
ORG=hg18

mysql --user genome --host genome-mysql.cse.ucsc.edu  -NAD $ORG -e \
    "select chrom, txStart, txEnd, X.geneSymbol, strand from knownGene as K, kgXref as X WHERE txStart != txEnd AND X.kgID = K.name" \
    | awk -v ups=$UPSTREAM -v ins=$INSTREAM 'BEGIN{OFS=FS="\t"}
           $5 == "-" { print $1,$3-ins,$3+ups,$4 } 
           $5 == "+" { print $1,$2-ins,$2+ups,$4 }' \
    | sort -k1,1 -k2,2n \