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

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

友情链接更多精彩内容