2021-09-29

gffcompare 这个报错真让我 yue 了!

报错

$ gffcompare -r genome/genes.gtf -o stringtie -i stringtie/stringtie_merged.gtf
Error: cannot locate input file: 
1      StringTie       transcript      998962  1000172 1000    -       .       gene_id "MSTRG.1"; transcript_id "ENST00000428771"; gene_name "HES4"; ref_gene_id "ENSG00000188290";

不对,怎么报错了?我前面也用了这个命令,能成功运行啊?

经过

那有错也还得解决啊,还能咋办?

把“老三套”独立解决报错的办法搬出来:

(1)搜索引擎检索

基本上啥也检索不到,唯一相关的链接上,只有问题,没有答案。
相关链接:https://github.com/gpertea/gffcompare/issues/3

(2)软件的帮助文档

找不到有用信息。

(3)重装一遍

没有用

反思

以上就是我觉得 rue 了的原因,但谁又能想到小丑竟是我自己!

上面的方法全部试个遍后,仍然没有找到解决办法,那就不能再相信前面成功运行过的命令了,找找别人是怎么用 gffcompare 的。

然后我就发现,-i 选项原来不是那样子用的。

正在我想写文章分享,梳理流程时,发现原来“老三套”中早就提示过了 。

(1)搜索引擎检索

相关链接:https://github.com/gpertea/gffcompare/issues/3
相关链接上是有写的:

(2)软件的帮助文档

帮助文档也是有写的:

没看见啊,没看见!

解决办法

解决办法:

  • 单个gtf 文件:不需要用 -i,直接接 gtf 文件;
  • 多个 gtf 文件:修改为 -i gtf_list.txt;
gffcompare -r genome/genes.gtf -o stringtie stringtie/stringtie_merged.gtf
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

  • 今天尝试录制读书类产品,感觉流量数据不太好,但先做吧,至少迈出了第一步,尽量字数控制在,300字左右,正好是录制1...
    努力生活呀阅读 1,247评论 0 0
  • 2021年9月28日,正向分享302天 感激生活的赐予,让爱融入生活 收尾工作断断续续的完成,节前各项检查都如期进...
    钟楚君阅读 1,171评论 0 0
  • 阅读蒙顶茶-百度百科[https://baike.baidu.com/item/%E8%92%99%E9%A1%B...
    蔚空阅读 1,670评论 0 0
  • 16宿命:用概率思维提高你的胜算 以前的我是风险厌恶者,不喜欢去冒险,但是人生放弃了冒险,也就放弃了无数的可能。 ...
    yichen大刀阅读 11,521评论 0 4
  • 公元:2019年11月28日19时42分农历:二零一九年 十一月 初三日 戌时干支:己亥乙亥己巳甲戌当月节气:立冬...
    石放阅读 11,899评论 0 2

友情链接更多精彩内容