gtf注释文件gene_id、gene_name、gene_biotype的合并

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即可。

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。