Hi,

I am trying to use a software require full dbsnp annotation from UCSC with gtf format.

The uscs table browser offer a online server, so that we can export its with gtf format, but Its not work when file size is large

I also try to download dbsnp from UCSC, eg: https://hgdownload.cse.ucsc.edu/goldenpath/hg38/database/snp151.txt.gz, and convert it by genePredToGtf, but its also not work because of

invalid unsigned integer: "A"

So, how can I get the gtf format dbsnp data from UCSC?

genePredToGtf cannot work, as it convertes the genePred format to GTF so the input must be in genePred format. The dbSNP tables are not in genePred format as they are not gene predictions. GTF is a bad match for this type of data, REDItools made a weird choice, as GTF has fields for exons and frames that don't make sense for SNPs. Either way you can use awk to convert the file to a very bare-bones GTF file which is probably good enough for this software:

curl -s https://hgdownload.cse.ucsc.edu/goldenpath/hg38/database/snp151.txt.gz | zcat | awk '{OFS="\t"; print ($2,"ucsc", "snp", $3,$4, "1", "+", "0", "")}' > snp151.gtf