美文网首页
gtf注释文件gene_id、gene_name、gene_bi

gtf注释文件gene_id、gene_name、gene_bi

作者: 啊吧啊吧鸭 | 来源:发表于2020-09-20 17:38 被阅读0次

    gtf注释文件中,第9列的内容并不是完全对齐的,所以用cut并不能很好的取出这三列内容。当然gene_id这列除外。

    图一

    这个时候可以考虑awk的匹配模式了

    cat test | awk 'BEGIN{FS=";"} {for(i=1;i<= NF;i++) if($i~/gene_name/) {print$i}}' > gene_name.txt
    
    

    这样分别取出gene_id、gene_name、gene_biotype并重定向到三个文件中,在paste即可。

    相关文章

      网友评论

          本文标题:gtf注释文件gene_id、gene_name、gene_bi

          本文链接:https://www.haomeiwen.com/subject/ucrpyktx.html