创建gtf下载

最新推荐文章于 2024-01-09 12:08:52 发布

NeoBilly

最新推荐文章于 2024-01-09 12:08:52 发布

阅读量1.2k

点赞数

分类专栏：随手记

本文链接：https://blog.csdn.net/NeoBilly/article/details/17715199

版权

随手记专栏收录该内容

11 篇文章 0 订阅

订阅专栏

搞错俩，麻利的改回来

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 loci.tracking | perl -ne '@t = split(/\t/, $_); $id1 = $t[0]; $t[4] =~ /.*?\|(.*?)\|/; $id2 = $1; print $id1 . "\t" . $id2 . "\n";' | perl /leofs/noncode/NONCODEv4/cmds/leftJoin.pl - 2 id.v4_tcons 2 | awk '{print $3 "\t" $1 ;}' | perl -ne '@t = split(/\t/, $_); chomp @t; $nId = "NONHSAT" . substr("000000", 0 , 6 - length($t[0]."") ) . $t[0]; print $t[1] . "\t" . $nId . "\n";' | perl /leofs/noncode/NONCODEv4/cmds/id_replace_tcons.pl - gtf.v4.human.lnc.1 | sed -e 's/oId .*//' > human.hg19.v4.lnc.gtf