gffread 常用使用自己速记

gffread gff3文件 -g 基因组文件 -x cds文件 -y pep 文件

gffread gff3文件 -T -o gtf文件
gffread gtf文件 -o gff文件

大家常见的-o- 是指定输出打印到屏幕，所以也可以写成
gffread gtf文件 -o- > gff文件

今天gffread 翻译后的蛋白投入 orthofinder 发现 orthofinder 在调用diamond报错，错误内容是无法识别 “.“字符，这是因为gffread 把stop codon翻译为”.”(dot)，而不是我们使用mega时候的""(asterisk)。
gffread 可以在使用-y pep 时候 -S "" 指定终止密码子为"*" ,但是作者在github写So ugly. 人们为什么看到星号不是通配符？他就觉得很丑。

我呢懒得重新翻译了，所以就写了一个小脚本更改一下：

stop_codon_replace.py

import sysdef replace(infile):with open(infile) as fi:for line in fi:line = line.strip()if line.startswith(">"):print(line)else:   if "." in line:line = line.replace(".","*")print(line)else:print(line)if __name__ == "__main__":replace(sys.argv[1])

然后 for循环就好了
for i in *pep; do python stop_codon_replace.py $i> $(basename $i).change ; done

gffread 常用使用自己速记相关推荐

jQuery速记Ajax方法简介
如果您从未听说过Ajax这个词,请举手. 我敢打赌,几乎所有人的手臂都向下垂,靠近他们的身体. Ajax 最初代表异步JavaScript和XML ,它是最常用的客户端方法之一,可帮助创建异步网站和W ...
Windows轻松实现linux shell环境：gitforwindows
Windows缺少shell命令支持用过Linux服务器分析数据的小伙伴,一定对Linux强大Shell命令所折服,经常会感觉windows缺少这些命令而感觉不方便. 还有想学习Linux Shel ...
Windows不用虚拟机或双系统，轻松实现linux shell环境：gitforwindows
windows缺少shell命令支持用过Linux服务器分析数据的小伙伴,一定对Linux强大Shell命令所折服,经常会感觉windows缺少这些命令而感觉不方便. 还有想学习Linux Shel ...
服务器windows模拟linux环境,科学网—Windows不用虚拟机或双系统，轻松实现shell环境：gitforwindows - 刘永鑫的博文...
windows缺少shell命令支持用过Linux服务器分析数据的小伙伴,一定对Linux强大Shell命令所折服,经常会感觉windows缺少这些命令而感觉不方便. 还有想学习Linux Shel ...
Scala 2.8馆藏图书馆是“历史上最长的遗书”吗？ [关闭]
我刚刚开始研究即将发布的2.8版本中的Scala集合库重新实现 . 熟悉2.7中的库的人会注意到,从使用角度来看,库几乎没有变化. 例如... > List("Paris", ...
学习正则表达式的简单方法
What is Regular Expression? 正则表达式是一组字符或符号,用于从文本中查找特定的模式. 正则表达式是一种从左到右与主题字符串匹配的模式. "正则表达式"一 ...
TensorFlow2中的常用API速查速记
这篇博客为<深度学习与TensorFlow 2入门实战>的部分笔记,为上一篇"TensorFlow2中的Tensor的常用操作速查速记"的后续,今天也正好看见龙龙老师的 ...
Matlab 常用语法速记 1
===快捷帮助=== 查找帮助:doc 函数名强制终止: ctrl + c 注释: % 清屏 clc,清临时变量窗口 clear NaN: not a number inf: infinite == ...
linux常用命令速记
一.命令提示符说明 1. [root@localhost ~]# root: 当前登录用户 localhost: 主机名 ~: 当前所在目录 #: 超级用户提示符($: 普通用户) 2. -rwxr- ...

gffread 常用使用自己速记

gffread 常用使用自己速记相关推荐

最新文章

热门文章