常用

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
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容