[BiO BB] Getting sequences by base pair locations
MEC at Stowers-Institute.org
Fri Jul 28 12:19:42 EDT 2006
There are many options.
But, it looks like you only have start end coordinates! Where do you
know which chromosome/contig the hit was on?
Assuming you have this, if you did the blat with a local copy of the
blat program and a the genome, then in addition to the blat command, you
have the twoBitToFa command which can extract the hits from the blat
index (see http://genome.ucsc.edu/goldenPath/help/blatSpec.html)
Or did you do the blat at ucsc?
Database Applications Manager, Bioinformatics
Stowers Institute for Medical Research
bio_bulletin_board-bounces+mec=stowers-institute.org at bioinformatics.org
[mailto:bio_bulletin_board-bounces+mec=stowers-institute.org at bioinformat
ics.org] On Behalf Of Yuval Itan
Sent: Thursday, July 27, 2006 11:23 AM
To: bio_bulletin_board at bioinformatics.org
Subject: [BiO BB] Getting sequences by base pair locations
I was BLATing a few hundred human genes against the chimp
genome, and kept the best chimp hits for every human gene.
I have the base pair start and end location for every chimp hit,
and I need to get the sequence for each of these chimp hits. Here is an
example for a few chimp hits bp locations:
I have one chimp genome file (about 3GB) including all
chromosomes, but I could also get one file per chromosome if that would
make things easier. Does anyone have a script or a link for an interface
that can do the job?
Thank you very much,
More information about the BBB