blat command line output completely empty but webserver output correct

5 views
Skip to first unread message

Federica Luppino

unread,
Jun 11, 2025, 12:53:08 PMJun 11
to genome...@soe.ucsc.edu
Hello,

I have an important topic concerning the comparison between blat command line and webserver version at (Human BLAT Search). 

I am using BLAT v38 (and I have reproduced the error as well with version v39). As suggested in the FAQ, I am using the following command line arguments to reproduce the webserver results:

blat -stepSize=5 -repMatch=2253 -minScore=20 -minIdentity=0 database.2bit query.fa output.psl

The command I am using is:
blat ../reference/hg38.2bit try_chr7.fa -stepSize=5 -repMatch=2253 -minScore=20 -minIdentity=0 out.psl

Some examples on the fasta sequence of the try_chr7.fa:
>chr7:38290616-38291666/rc_sliding:6-25
GAAAGAAGATGAGGTGGTGT
>chr7:38290616-38291666/rc_sliding:10-29
GAAGATGAGGTGGTGTACCA
>chr7:38290616-38291666/rc_sliding:11-30
AAGATGAGGTGGTGTACCAC
>chr7:38290616-38291666/rc_sliding:40-59
TGGATTAGGCACCACAGTGT

The output file out.psl is completely empty, while if I run those sequences on the webserver I get the correct results. In addition, the problem is only related to sequences on chromosome7.

Would you be able to assist me with that? 

Thank you a lot,

Dr. Federica Luppino

Reply all
Reply to author
Forward
0 new messages