本文介绍了一种使用Perl脚本进行基因注释文件转换的方法，包括从原始注释文件到人类基因组注释文件的转换过程。
cat id.v4_xloc | perl -ne '@t = split(/\t/, $_); chomp @t; $nId = "NONHSAG" . substr("000000", 0 , 6 - length($t[0]."") ) . $t[0]; print $t[1] . "\t" . $nId . "\n";' | perl /leofs/noncode/NONCODEv4/cmds/id_replace.pl - loci.combined.gtf > gtf.v4.human.lnc.1
cat id.v4.tcons | awk ' {print $2 "\t" $1; }' | perl /leofs/noncode/NONCODEv4/cmds/id_replace_tcons.pl - gtf.v4.human.lnc.1 | sed -e 's/oId .*//' > human.hg19.v4.lnc.gtf
scp human.hg19.v4.lnc.gtf.gz root@112.124.13.146:/data/var/www/html/NONCODE/datadownload/
cat id.v4.tcons | awk ' {print $2 "\t" $1; }' | perl /leofs/noncode/NONCODEv4/cmds/id_replace_tcons.pl - gtf.v4.human.lnc.1 | sed -e 's/oId .*//' > human.hg19.v4.lnc.gtf
scp human.hg19.v4.lnc.gtf.gz root@112.124.13.146:/data/var/www/html/NONCODE/datadownload/
扫一扫
216
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
