因缺少模块儿rnammer安装失败,转使用barrnap预测基因组上的核糖体RNA

一、预安装bedtools

因为barrnap依赖于bedtool和nhmmer,所以额就傻到先去安装bedtools了,结果从github上下载最新版本2.29,报错,因为额的系统版本低,软件版本高,不兼容。


这是我的报错提示

所以额就退而求其次下载2.25版本了,手动安装,还挺顺利的。(代码从github上获取)

wgethttps://github.com/arq5x/bedtools2/archive/v2.25.0.tar.gz

tar xzvf v2.25.0

cd bedtools2-2.25.0/

make

cd bin/

export PATH=$PWD:$PATH        #直接这样输路径好像不管用,并没有将路径输进去,也有可能是我操作问题

echo $PATH    #先查看一下有没有路径,后来再来一遍方便对照

vim ~/.bashrc    #我是打算直接永远生效的,所以vim修改

在最后一行加上:补充上这个软件的安装路径到bin为止

source ~/.bashrc    #source之后,直接生效

echo $PATH          #此时在返回的最前面可以看到路径

输入bedtools(或者加上-h),如果返回有提示,而不是出现command not found,那么就算是安装好了。

插一句话:这两种export的格式应该是一样的吧?

export PATH=$PATH:~/software/eggnog-mapper/bin

export PATH=~/transcriptome/soft:~/conda/bin:$PATH

参考这篇路径的配置方式linux可执行文件添加到PATH环境变量的方法https://www.cnblogs.com/henrylinux/p/9746818.html

二、conda直接安装(大转折!额真是太意外了,毕竟要买智商税!)

当我满心欢喜的要去安装我们的主题barrnap的时候,我上去github上,结果此网页下面说可以conda安装。嗯,是的,明明一行代码就能解决的事情,我为什么要手动呢?!!主要是,当我看到barrnap下面要求bedtools的版本要>=2.7的时候,我哭聊。

conda install -c bioconda -c conda-forge barrnap

不过好在conda好就好在,你拿它安装,你欠缺的辅助软件只要你yes,他会帮你自动安装,我的就是自动安装的,帮我安装上了bedtools-2.28.0版本,省的我哭了;还有就是不用管配置,直接帮我配置到我家目录下的conda的bin下

barrnap -h        #有提示,说明我安装成功了

which bedtools        #查看一下我的路径

  ~/conda/bin/bedtools        #显示在这里

还有之前说的nhmmer,应该也是随着一块儿下了,

nhmmer -h    #有提示

which nhmmer        #返回的路径~/conda/bin/bedtools 

脑残的我之前以为nhmmer是hmmer,而且之前我确实检索服务器,确实没有nhmmer(当时竟然蠢到不行手动安装,但是我百度也没百度到)。

三、barrnap使用

参数的话可以 --help查看

--kingdom参数指定物种类型,bac代表细菌,arc代表古菌,euk代表真核生物,mito代表后生动物线粒体;--threads指定并行的线程数,预测结果以GFF3格式保存。

(关于线程数的浅解:线程数,就是核心数,跟人的脑子一样,核心数2就说明CPU有两个脑子。脑子越多解决问题速度越快。CPU的核心数越高处理速度就越高。核心数2通俗地说就是双核CPU了。但自超线程技术问世后,一个核心可以同时2个线程了。使CPU性能上升百分之40。----摘自百度百科

但是也不是线程数越高越好,要看服务器的配置)

barrnap --kingdom bac --threads 8 --quiet SXXXXXXXX011.fas > 011rRNA.gff3 

#1)其实我在纠结输入的文件格式fas是不是对的;2)在所要分析的文件目录下运行

 barrnap --kingdom bac --threads 8 --quiet S115.fas > 115rRNA.gff3

 barrnap --kingdom bac  --threads 8 --quiet S114.fas > 114rRNA.gff3

 barrnap --kingdom bac  --threads 8 --quiet S108.fas > 108rRNA.gff3

 barrnap --kingdom bac  --threads 8 --quiet S111.fas > 111rRNA.gff3

按道理该得到三个文件,但是我没有得到5s的文件,所以不知道为什么。可能是我的序列的问题,也可能是做个饭软件的问题吧(毕竟我搜的时候这个软件的相关信息就很很少,感觉很小流的样子)。

参考:

1、barrnap:预测基因组上的核糖体RNAhttps://www.jianshu.com/p/315a7bcc93dc

2、bedtools 的安装与使用庐州月光https://www.cnblogs.com/xudongliang/p/5051503.html

关于这块儿基本也就这些了,后期的分析还没做到,再补一期分析。

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

推荐阅读更多精彩内容