This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Table Browser refGene identifiers question

Hi,

I want to put a list of genes into table browser (identifiers tab) and get a BED file with whole gene intervals however when I input 10K gene names into Genes and Gene Predictions I get a file out with 17K intervals. Is there any straightforward way to keep it to the 10K genes only and not include the predictions as well? I am assuming thats where the additional numbers come from. Thanks.

Rob.

table browser ucsc
• 1,579 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by rbronste rbronste

what kind of identifiers is it ? ucsc id ? ensembl-id ? common names ?

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Its common names I believe, those spit out by cuffdiff.

• 0 views • link
written by rbronste rbronste

1 answer

use mysql ucsc, something like:

$ join -t $'\t' -1 1 -2 4 <(cat geneids.txt | LC_ALL=C sort | uniq ) <(mysql --user=genome --host=genome-mysql.cse.ucsc.edu -A -D hg19 -N -e 'select chrom,txStart,txEnd,name2  from wgEncodeGencodeBasicV7 limit 10' | LC_ALL=C sort -t $'\t' -k4,4 | uniq )
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this