生信 - 从repeatmasker传送门过来的 blast
以前有的是非完整时间写的博客,抽时间需要统一整理一下。
今天在重新装repeatmasker。
整个过程是这样的,有关联的事情有两个。
1. 装repeatmasker需要各种Prerequisites,其中就可能用到了blast,而之前一直找这个版本的blast,在ncbi硬是没有找到:
- For RMBlast ( NCBI Blast modified for use with RepeatMasker/RepeatModeler ) please go to our download page: http://www.repeatmasker.org/RMBlast.htm
blastdb_aliastool blastn deltablast makeblastdb rmblastn tblastn
blastdbcheck blastp dustmasker makembindex rpsblast tblastx
blastdbcmd blastx get_species_taxids.sh makeprofiledb rpstblastn update_blastdb.pl
blast_formatter convert2blastmask legacy_blast.pl psiblast segmasker windowmasker
2. 我有一些基因,需要找一下这些基因在各个物种的基因组上的位置,并注释出来。
然后我用了ncbi的版本:好吧,区别不到,少了 rmblastn
3. 安装repeatmasker
按照路径指示没有问题
4.安装repeatmodeler
perl ./configure
#
The following perl modules required by RepeatModeler are missing from
your system. Please install these first:
JSON
没看懂,去找,没找到怎么下的
$ wget https://centos.pkgs.org/7/centos-x86_64/perl-JSON-2.59-2.el7.noarch.rpm.html
$ rpm2cpio perl-JSON-2.90-1.of.el7.src.rpm |cpio -div
$ tar zxvf JSON-2.90.tar.gz
$ perl Makefile.PL
make
make install
完事
4.继续装repeatmodeler,
**REPEATMASKER INSTALLATION PATH** 将repeatmasker的路径喂给它
**RECON INSTALLATION PATH** 找找这个是啥
RECON - De Novo Repeat Finder, Bao Z. and Eddy S.R. Developed and tested with our patched version of RECON ( 1.08 ). The 1.08 version fixes problems with running RECON on 64 bit machines and supplies a workaround to a division by zero bug along with some buffer overrun fixes. The program is available at: http://www.repeatmasker.org/RECON-1.08.tar.gz. The original version is available at http://eddylab.org/software/recon/.
**RepeatScout INSTALLATION PATH**
RepeatScout - De Novo Repeat Finder, Price A.L., Jones N.C. and Pevzner P.A. Developed and tested with our multiple sequence version of RepeatScout ( 1.0.6 ). This version is available at http://www.repeatmasker.org/RepeatScout-1.0.6.tar.gz
http://bix.ucsd.edu/repeatscout/
**NSEG INSTALLATION PATH**
http://www.pdas.com/packages/nseg.zip
不对
ftp://ftp.ncbi.nih.gov/pub/seg/nseg/.
然后就可以了。
继续报错:
exist! RepeatModeler will not run correctly without these
然后去看贴吧
'm only guessing, but [n??] sounds like the extensions of files generated by formatdb when creating nucleotide database. Given that RepeatMasker.lib is a FASTA file, maybe you need to run formatdb on it.
./makeblastdb: /lib64/libc.so.6: version `GLIBC_2.14' not found (required by ./makeblastdb)
所以可能存在的问题有两个,一个是没有建库,一个是建库软件错误,无法建库。
不行!!!继续找别人的经验
不想折腾了,下载RepeatModeler-2.0.tar.gz
注意,如果不想折腾的话,不要选yes,我去下载LTR_Retriever,MAFFT
说实话,我想回滚了。
wget http://mafft.cbrc.jp/alignment/software/mafft-7.158-without-extensions-src.tgz (备份一个wget http://mafft.cbrc.jp/alignment/software/mafft-7.158-gcc_fc6.x86_64.rpm)
按照说明安装和编译,repeatmodeler说不想玩,下另一个,OK
https://github.com/traviswheelerlab/ninja/releases/tag/0.95-cluster_only
make all
应该通关了吧。
Congratulations! RepeatModeler is now ready to use.
今天准备将所有的常用软件备份到自己的网站undef.cn里面,但是这个网站还没有备案,所以可以先在cospure.cn上看到。
以上,repeatmodeler安装通关记录
abyss king
生信 - 从repeatmasker传送门过来的 blast的更多相关文章
- 生信工具汇总--OMICtools
各种生信工具: https://omictools.com/
- 生信基础知识【04】GO和pathway分析
非原创 参考资料: 一文掌握GO和pathway分析 - 生物信息学讨论版 -丁香园论坛http://www.dxy.cn/bbs/thread/34904124#34904124 GO富集 GO是G ...
- 生信软件的好帮手-bioconda--转载
http://mp.weixin.qq.com/s/nK1Kkf9lfZStoX25Y7SzHQ 这篇文章主要适用于Linux平台,当然MacOS也行,不过它有更好安装方法. 此外网上也会许多更好的关 ...
- 精心整理(含图版)|你要的全拿走!(R数据分析,可视化,生信实战)
本文首发于“生信补给站”公众号,https://mp.weixin.qq.com/s/ZEjaxDifNATeV8fO4krOIQ更多关于R语言,ggplot2绘图,生信分析的内容,敬请关注小号. 为 ...
- 生信-使用NCBI进行目的基因的引物设计
使用NCBI进行目的基因的引物设计 全文概述 利用生信工具进行目的基因的引物设计,使用了NCBI进行筛选与设计引物,使用 idtdna对筛选出的DNA进行检查.本文分享了如何筛选出高质量的基因引物,帮 ...
- 计科菜鸟玩生信(一)——Windows10下用docker安装GATK
1.官网下载GATK. (学校网速实在是太慢了,下载几个小时到自闭) 下载地址:https://software.broadinstitute.org/gatk/download/ 下载完成后文件夹中 ...
- 生信云实证Vol.12:王者带飞LeDock!开箱即用&一键定位分子库+全流程自动化,3.5小时完成20万分子对接
LeDock是苏黎世大学Zhao HongTao在博士期间开发的一款分子对接软件,专为快速准确地将小分子灵活对接到蛋白质而设计.LeDock优于大部分商业软件,在Astex多样性集合上实现了大于90% ...
- 生信分析常用脚本(二)--SOAPdenovo
1.SOAPDenovo配置文件示例 软件下载安装和使用:http://soap.genomics.org.cn/soapdenovo.html asm.cfg #maximal read lengt ...
- 生信-RPKM与FPKM学习[转载]
转自:https://baike.baidu.com/item/RPKM/1197657 均反应基因的表达水平 1.RPKM的计算公式 分母是总共比对到这个基因的reads的数目(条 为单位),分母是 ...
随机推荐
- Vue中的循环以及修改差值表达式
0828自我总结 一.Vue中的循环 v-for 常见的4总情况 #第一种 <div v-for="item in items"></div> #第二种 & ...
- AWVS安全渗透扫描
1.打开软件,点击New Scan 2.在website url中输入被扫描的网址,点击next 3.在scanning profile中选择测试的漏洞类型,默认选择default(默认) 在scan ...
- Java微服务(二):负载均衡、序列化、熔断
本文接着上一篇写的<Java微服务(二):服务消费者与提供者搭建>,上一篇文章主要讲述了消费者与服务者的搭建与简单的实现.其中重点需要注意配置文件中的几个坑. 本章节介绍一些零散的内容:服 ...
- MySQL make_set()的用法
MAKE_SET(bits,str1,str2,…)返回一个设定值(含子字符串分隔字符串","字符),在设置位的相应位的字符串.str1对应于位0,str2到第1位,依此类推.在s ...
- 毕业设计过程中的一些学习Android网站
安卓巴士:http://www.apkbus.com/CSDN:(下载资源)http://www.csdn.net/?ref=toolbar博客园:http://www.cnblogs.com/极客学 ...
- 利用python模拟菜刀反弹shell绕过限制
有的时候我们在获取到目标电脑时候如果对方电脑又python 编译环境时可以利用python 反弹shell 主要用到python os库和sokect库 这里的服务端在目标机上运行 from sock ...
- Docker卷
要解决的问题:在移除了现有容器或者添加了新容器时,之前容器的数据无法访问. 为避免上述问题,Docker提供了以下策略来持久化数据 tmpfs挂载 绑定挂载 卷 1.tmpfs挂载 2.绑定挂载 将D ...
- ElasticSearch业务逻辑案例
ElasticSearch业务逻辑案例 一.业务难题 我们有一个索引: myindex/mytype(为了方便,我们下文以a/b表示) 索引类型中的一个字段group之前是a.b.c(历史遗留问题), ...
- Js极客之路 - 优化操作(性能优化)
1.因为每次For循环都会计算一次arr.length,所以有必要存储数组长度以减少计算.针对这篇文章(http://www.crimx.com/2015/04/21/should-array-len ...
- JavaScript 变量作用域和声明提升
一.变量作用域 说到这个概念,不有自主的想到this,scope 这两个关键字. JavaScript的this总是指向一个明确的对象,这个对象是在执行的时候动态绑定的.通俗的说就是谁调用我,我的th ...