创建gtf下载

最新推荐文章于 2024-12-12 21:02:47 发布

原创最新推荐文章于 2024-12-12 21:02:47 发布 · 1.3k 阅读

0 ·

CC 4.0 BY-SA版权

随手记专栏收录该内容

11 篇文章

订阅专栏

本文介绍了一段Perl脚本用于修改基因注释文件中的ID，包括替换和整合不同来源的数据，最终输出格式化的GTFFILE。通过脚本的执行，实现了从原始数据到标准化基因ID的转换过程。

搞错俩，麻利的改回来

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