利用毒株的accession号,批量下载fasta格式的文件
一、安装好Entrez Driect
见:http://www.cnblogs.com/lmt921108/p/8087474.html
二、在NCBI官网找到需要下载的文件的accession号
二.创建bash脚本文件
1.下载全基因组fasta序列(get_comseq.sh)
1 #!/bin/bash 2 3 cat $1 | while read line 4 do 5 efetch -db nucleotide -id $line -format fasta > $line.fasta 6 done 7 8 #利用bash的for循环实现fasta文件的批量下载 9 #局限,需要知道毒株的acession号 10 #代码解释,第三行:逐行读取文件,得到变量line--保存了每行文件的内容。
运行:bash get_comseq.sh accession.txt
2.下载各个区段的碱基序列(get_genecds.sh)
1 #!/bin/bash 2 3 cat $1 | while read line 4 do 5 efetch -db nucleotide -id $line -format gene_fasta > gene$line.fasta 6 done
运行:bash get_genecds.sh accession.txt
3.下载各个区段的氨基酸序列(get_protcds.sh)
1 #!/bin/bash 2 3 cat $1 | while read line 4 do 5 efetch -db nucleotide -id $line -format fasta_cds_aa > cds$line.fasta 6 done
运行:bash get_protcds.sh accession.txt
四、accession文件格式
AB007451.1 AB022288.1 AB022289.1 AB548618.1 AB548619.1
一行只能有一个号
转载于:https://www.cnblogs.com/lmt921108/p/8136417.html
利用毒株的accession号,批量下载fasta格式的文件相关推荐
- 通俗易懂【Springboot】 单文件下载和批量下载(多个文件合成一个压缩包下载)
文章目录 一.单文件下载 1.简单理解文件下载 2.单文件下载的具体代码实现 3.测试 4.单文件下载整体代码 二.多文件批量下载(多个文件合成一个压缩包下载) 1.多文件下载的实现方式,这里使用了Z ...
- 【数据下载】使用DownThemAll批量下载网站上的文件
1 下载DownThemAll扩展程序 网址:DownThemAll! 2 批量下载网站上的文件 2.1 全球人工不透水面积GAIA数据 以全球人工不透水面积GAIA数据为例,网页界面如下: 右键后, ...
- element-ui实现下载动态进度条,可实现批量下载,附Blob文件类型对照表
element-ui实现下载动态进度条,可实现批量下载,附Blob文件类型对照表 ①页面结构使用element-ui的进度条组件,示例为圆形,具体样式根据需求确定 <el-button>下 ...
- 【模板】批量下载视频到指定文件脚本
批量下载视频到指定文件脚本 以下代码是在学习工作中总结出来的方法所构成的模板,仅适合在某些特定情况下使用,如下模板需要拥有视频的URL地址,相关三方库,稳定的网络等 需要库: you-get os f ...
- 下载torrent格式的文件
下载torrent格式的文件 下载uTorrent 这个可以从网上下载就好啦. (总感觉这个软件有点奇怪) 打开torrent文件 现在我们打开torrent文件就可以进行下载了,因为系统已经进行识别 ...
- 如何优雅的批量下载m3u8 格式视频
我们去很多网站上 去下载视频,解析出的地址 是m3u8格式视频, 使用 普通的下载方式下载就是一个m3u8文件.里面装的是一段一段的视频.一方面,不做特殊处理播放器可能播放不了,另一方面该格式可能是加 ...
- 利用wget批量下载http目录下文件
因为网络特殊,连不上互联网,只好自己制作一个puppet本地源.因为自己的电脑是是可以连互联网的,所以就直接将puppetlabs.repo导入进来了 rpm -Uvh http://yum.pupp ...
- 前端接收bolb格式下载文件,前端下载blob格式的文件
最近遇到一个下载文件的需求,不是请求后端接口直接下载文件,而是后端返回blob格式的文件,前端自行下载文件 下面的具体的操作代码 1.修改request请求(请求是封装好的axios方法),增加响应参 ...
- 【InSAR 笔记1】ASF网站哨兵一号批量下载
加入下载筐 -> 打开下载筐 -> 下载.py文件 -> 运行.py文件 -> 输入用户名和密码 ASF 下载网站 python download-all-2020-05-12 ...
最新文章
- ubuntu 18 mysql5.7_Ubuntu18.04 下 MySQL5.7 的安装
- 利用Object.defineProperty挂接set钩子,监控对象属性的修改事件
- linux ping监控脚本,WINDOWS和LINUX下带时间的PING包监控脚本
- android获取子线程id,Android 开发 知晓各种id信息 获取线程ID、activityID、内核ID
- :Android网络编程--XML之解析方式:SAX
- MySQL 安装手册(5.7.20)
- 蓝桥杯java 大纲,2019 第十届蓝桥杯Java省赛B组个人总结
- 靠,竟然有如此沙雕的代码注释!
- jsdroid 教程_南方Cass专题,全系列教程+插件汇总打包
- php bi系统,bi系统是什么系统
- 串联和并联拓扑PI控制器
- proteus三输入与非门名字_proteus元件名称中文名说明
- 创业公司,老板说等公司做大了,给5%股权,建议你不要轻易相信
- 复合赋值运算符“+=、-=、*=、/=、%=”详解
- Android使用Github Actions持续集成并自动上传apk到蒲公英App内测分发平台(含证书密码脱敏)
- Git的稀疏检出功能
- 常用URL分享,实用地址
- windows下wav格式歌曲文件信息修改
- 【genius_platform软件平台开发】第二十八讲:NEON指令集优化(附实例)
- 公司绝对不会告诉你的潜规则 - 看到职场如战场了吧
热门文章
- http服务(nginx、apache)停用不安全的SSL协议、TLS1.0和TLS1.1协议/启用TLS1.3
- 内部系统自己搭建还是花钱购买?这是一个问题
- Android 自动换行添加控件
- Golang-Beego(go mod)
- linux复制整个目录_如何在Linux中复制整个目录?
- php安全开发正则表达式,动态网页制作PHP常用的正则表达式
- 伦敦交通局设置较低的速度限制
- 论文笔记:HIE-SQL:History Information Enhanced Network for Context-Dependent Text-to-SQL Semantic Parsing
- python_turtle库丘比特之箭
- 移动直播技术秒开优化经验