gtf转bed12
conda install -c bioconda ucsc-gtftogenepred
gtf=/path/to/gencode.v44.annotation.gtf
genePred=/path/to/gencode.v44.annotation.genePred
bed12=/path/to/gencode.v44.annotation.bed12
gtfToGenePred -genePredExt -geneNameAsName2 $gtf $genePred
awk '{print $2"\t"$4"\t"$5"\t"$1"\t0\t"$3"\t"$6"\t"$7"\t0\t"$8"\t"$9"\t"$10}' $genePred > $bed12
