[Genome] short results when running my own blat server
Alex Holman
alexander.holman at tufts.edu
Tue Oct 10 08:32:55 PDT 2006
Hi,
I'm trying to run a local copy of the BLAT server so I can hit it
with a ton of search requests to the human genome. I downloaded the
compiled blatSuite for OSX, and have downloaded chr20.fa from the
latest alignment of the human genome for testing. However when I
take one line randomly from chr20 (50 bases), paste it into a fasta
file and use it to blat search against chr20.fa, I only turn up a 37
or 38 base hit. If I use the same sequence in the web blat at UCSC I
get the full 50 base hit.
My command is this:
./blat -t=dna -q=dna -out=pslx ../genome/chr20.fa testseq.fa output.psl
though it doesn't seem to make any difference if I specify the target
and query types or if I use the -fine option or set -minMatch and -
minScore to 0.
Searching with 100 bases of sequence I get a full length length hit.
What settings can I change to better replicate the sensitivity that I
see with the web based tool?
Thanks
-Alex Holman-
More information about the Genome
mailing list