云之南

风声，雨声，读书声，声声入耳；家事，国事，天下事，事事关心

日志

关于我

云之南

专业背景：计算机科学研究方向与兴趣: JavaEE-Web软件开发, 生物信息学, 数据挖掘与机器学习, 智能信息系统目前工作: 基因组, 转录组, NGS高通量数据分析, 生物数据挖掘, 植物系统发育和比较进化基因组学

文章分类

如何在本地BLAST【转】

2011-01-12 08:42:05| 分类：生物信息学 | 标签： |举报 |字号大中小订阅

下载LOFTER 我的照片书 |

1. 下载Standalone BLAST executables
下载地址: http ftp
注意下载和你的系统对应的文件
2. 解压
以下步骤针对于linux操作系统,windows操作系统的方法请参考解压后的doc文件夹的index.html文件
3. 创建一个名为.ncbirc的文件
内容如下:
[NCBI]
Data=”path/data/”
其中，”path/data/”是Standalone BLAST程序data子目录的位置。如: Data=/root/blast/data
.ncbirc文件应保存在你调用Standalone BLAST程序的目录或者你的根目录
4. 创建BLAST database
下载或者收集BLAST搜索的数据库，通常是核酸序列或蛋白质序列，
通过Standalone BLAST程序的formatdb命令把这些序列格式化。如:
formatdb -i db -o F -p F -n “nr” -v 2000000000
各选项意义如下:
-i 输入文件，只能是一个文件，
-o Parse options （默认是F)
T – True: Parse SeqId and create indexes.
F – False: Do not parse SeqId. Do not create indexes.
Update: 如果需要用fastacmd通过SeqId获取序列，则必须使用-o选项;
如果只是BLAST，不需要获取BALST得到的序列，则不使用-o选项
-p 文件类型（默认是T)
T – protein
F – nucleotide [T/F] Optional
-n 数据库名称
不指出的话默认为输入文件名
-v 数据库每卷大小(单位是百万字母) （默认是0)
这个选项把大的FASTA文件分割成卷，每卷大小由-v后的参数指定，不超过2000000000，即20亿字母。
更多选项请参阅解压后的doc文件夹的formatdb.html文件或者man formatdb
关于错误解决方法，请参阅http://bioinformatics.ubc.ca/resources/tools/formatdb
5. 开始BLAST！
大家熟悉的”blastp”, “blastn”, “blastx”, “tblastn”, 或 “tblastx” 都通过Standalone BLAST程序的formatdb命令实现。
格式是:
blastall -p blastn -d nr -i QUERY -o out.QUERY -m 9
各选项意义如下:
-p 程序名
包括 blastp: 通过蛋白质序列搜索蛋白质序列数据库
blastn: 通过核酸序列搜索核酸数据库
blastx: 通过翻译后的核酸序列搜索蛋白质数据库
tblastn: 通过蛋白质序列搜索翻译后的核酸数据库
tblastx: 通过翻译后的核酸序列搜索翻译后的核酸数据库
-d 数据库名称
与formatdb中-n选项一致
-i 输入文件
不指明的话默认为STDIN
-o 输出文件
不指明的话默认为STDOUT
-m (blast2, blastall, blastall_old, blastcl3, blastpgp,impala, megablast, rpsblast)显示选项:
0 pairwise (default)
1 query-anchored showing identities
2 query-anchored, no identities
3 flat query-anchored, show identities
4 flat query-anchored, no identities
5 query-anchored, no identities and blunt ends
6 flat query-anchored, no identities and blunt ends
7 XML Blast output (not available for impala)
8 tabular (not available for impala)
9 tabular with comment lines (not available for impala)
10 ASN.1 text (not available for impala or rpsblast)
11 ASN.1 binary (not available for impala or rpsblast)
-T 输出html格式的文件
更多选项请参阅解压后的doc文件夹的formatdb.html文件，或者man blast
6. Standalone BLAST程序还带有其他命令，如
bl2seq: 比较2个序列相似性
megablast: 高度相似性核酸序列BLAST
等等
详情请参阅解压后的doc文件夹,或者man blast
2009.9.3 Update:
BENM同学也写了一系列BLAST教程，很好很详细，大家可以参考下：
Blast使用方法攻略(一)
Blast使用方法攻略(二)
Blast使用方法攻略(三)

评论这张

转发至微博

阅读(3069)| 评论(0)

历史上的今天

this.p={  m:2,
              b:2,
              loftPermalink:'',
              id:'fks_094071086080083070080087081095085094087074093087083071',
              blogTitle:'如何在本地BLAST【转】',
              blogAbstract:'<INS\><INS id=\"google_ads_frame2_anchor\" \></INS\></INS\><BR\>    <BR\>  <TABLE cellSpacing=\"0\"  cellPadding=\"0\" \>  <TBODY\>  <TR\>  <TD id=\"postmessage_7200\" \>1. 下载Standalone BLAST executables<BR\>下载地址: <A target=\"_blank\" rel=\"nofollow\" href=\"http://www.ncbi.nlm.nih.gov/blast/download.shtml\" \>http</A\> <A target=\"_blank\" rel=\"nofollow\" href=\"ftp://ftp.ncbi.nih.gov/blast/executables/LATEST/\" \>ftp</A\><BR\>注意下载和你的系统对应的文件<BR\>2. 解压 <BR\>以下步骤针对于linux操作系统,windows操作系统的方法请参考解压后的doc文件夹的index.html文件<BR\>3. 创建一个名为.ncbirc的文件<BR\>内容如下:<BR\>[NCBI]<BR\>Data=”path/data/”<BR\>其中，”path/data/”是Standalone BLAST程序data子目录的位置。如: Data=/root/blast/data<BR\>.ncbirc文件应保存在你调用Standalone BLAST程序的目录或者你的根目录<BR\>4. 创建BLAST database<BR\>下载或者收集BLAST搜索的数据库，通常是核酸序列或蛋白质序列，<BR\>通过Standalone BLAST程序的formatdb命令把这些序列格式化。如:</TD\></TR\></TABLE\>',
              blogTag:'',
              blogUrl:'blog/static/1869915420110128425843',
              isPublished:1,
              istop:false,
              type:0,
              modifyTime:1316772905740,
              publishTime:1294792925843,
              permalink:'blog/static/1869915420110128425843',
              commentCount:0,
              mainCommentCount:0,
              recommendCount:0,
              bsrk:-100,
              publisherId:0,
              recomBlogHome:false,
              currentRecomBlog:false,
              attachmentsFileIds:[],
              vote:{},
              groupInfo:{},
              friendstatus:'none',
              followstatus:'unFollow',
              pubSucc:'',
              visitorProvince:'',
              visitorCity:'',
              visitorNewUser:false,
              postAddInfo:{},
              mset:'000',
              mcon:'',
              srk:-100,
              remindgoodnightblog:false,
              isBlackVisitor:false,
              isShowYodaoAd:true,
              hostIntro:'专业背景：计算机科学                           \n\n研究方向与兴趣: JavaEE-Web软件开发,\n生物信息学, 数据挖掘与机器学习, 智能信息系统                                                                \n目前工作: 基因组, 转录组, NGS高通量数据分析, 生物数据挖掘, 植物系统发育和比较进化基因组学',
              hmcon:'1',
              selfRecomBlogCount:'0',
              lofter_single:'<iframe width="140" height="560" style="overflow:hidden;" src="http://www.lofter.com/mailEntry.do?blogad=1&blog" frameBorder="0"></iframe>'
            }

{list a as x}
    {if !!x}
    <div class="iblock nbw-fce nbw-f40">
      <a class="fc03 noul" target="_blank" hidefocus="true" href="http://blog.163.com/${x.visitorName}/">
      {if x.visitorName==visitor.userName}
      <img alt="${x.visitorNickname|escape}" onerror="this.src=location.f40" class="cwd bdwa bdc0" src="${fn1(x.visitorName)}&r=${visitor.imageUpdateTime}"/>
      {else}
      <img alt="${x.visitorNickname|escape}" onerror="this.src=location.f40" class="cwd bdwa bdc0" src="${fn1(x.visitorName)}"/>
      {/if}
      </a>
      <div class="cwd vname thide">
        {if x.moveFrom=='wap'}
          <a class="noul pnt" target="_blank" href="http://blog.163.com/services/wapblog.html?frompersonalbloghome"><span title="来自网易手机博客" class="iblock wapIcon"> </span></a>
        {elseif x.moveFrom=='iphone'}
          <a class="noul pnt" target="_blank"><span title="来自iPhone客户端" class="iblock iphoneIcon"> </span></a>
        {elseif x.moveFrom=='android'}
          <a class="noul pnt" target="_blank"><span title="来自Android客户端" class="iblock androidIcon"> </span></a>
        {elseif x.moveFrom=='mobile'}
          <a class="noul pnt" target="_blank" href="http://blog.163.com/services/emsblog.html?frompersonalbloghome"><span title="来自网易短信写博" class="iblock wapIcon"> </span></a>
        {/if}
        <a class="fc03 m2a"  target="_blank" hidefocus="true" href="http://blog.163.com/${x.visitorName}/">
          ${fn(x.visitorNickname,8)|escape}
        </a>
      </div>
    </div>
    {/if}
    {/list}

<#--最新日志，群博日志--> <#--推荐日志-->

<p class="fc06">推荐过这篇日志的人：</p>
    <div>
      {list a as x}
      {if !!x}
      <div class="iblock nbw-fce nbw-f40">
        <a class="fc03 noul" target="_blank" hidefocus="true" href="http://blog.163.com/${x.recommenderName}/">
        <img alt="${x.recommenderNickname|escape}" onerror="this.src=location.f40" class="cwd bdwa bdc0" src="${fn1(x.recommenderName)}"/>
        </a>
        <div class="cwd thide">
          <a class="fc03 m2a" target="_blank" hidefocus="true" href="http://blog.163.com/${x.recommenderName}/">
            ${fn(x.recommenderNickname,6)|escape}
          </a>
        </div>
      </div>
      {/if}
      {/list}
    </div>
    {if !!b&&b.length>0}
    <p  class="fc06">他们还推荐了：</p>
    <ul>
    {list b as y}
      {if !!y}
        <li class="rrb"><span class="iblock">·</span><a class="fc03 m2a" target="_blank" href="http://blog.163.com/${y.recommendBlogPermalink}/?from=blog/static/1869915420110128425843">${y.recommendBlogTitle|escape}</a></li>
      {/if}
    {/list}
    </ul>
    {/if}

<#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇，下一篇--> <#-- 热度 -->

{list a as x}
    {if !!x}
    <div class="hotItem iblock nbw-fce nbw-f40">
      <a class="fc03 noul" target="_blank" hidefocus="true" href="http://blog.163.com/${x.publisherUsername}/">
      {if x.publisherUsername==visitor.userName}
      <img alt="${x.publisherNickname|escape}" onerror="this.src=location.f40" class="cwd bdwa bdc0" src="${fn1(x.publisherUsername)}&r=${visitor.imageUpdateTime}"/>
      {else}
      <img alt="${x.publisherNickname|escape}" onerror="this.src=location.f40" class="cwd bdwa bdc0" src="${fn1(x.publisherUsername)}"/>
      {/if}
      </a>
      <div class="cwd vname thide">
        <a class="fc03 m2a"  target="_blank" hidefocus="true" href="http://blog.163.com/${x.publisherUsername}/">
          ${fn(x.publisherNickname,8)|escape}
        </a>
      </div>
      <a class="f-myLikeIcons hottype {if x.type==1} js-liketype{elseif x.type==2} js-reblogtype{elseif x.type==3} js-sharetype{else}{/if}" target="_blank" hidefocus="true" href="http://blog.163.com/${x.publisherUsername}/"> </a>
    </div>
    {/if}
    {/list}

<#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->

页脚

我的照片书 - 手机博客 - 下载LOFTER APP - 订阅此博客

云之南

导航

日志

如何在本地BLAST【转】

历史上的今天

最近读者

热度

评论

页脚