分割 fasta 文件
文章目录
- 一、分割 fasta 文件
一、分割 fasta 文件
效果:对 fasta 文件进行分割,每个 > 的序列作为一个文件保存,并以 > 后的字符作为文件名。
# -*- coding: utf-8 -*-
"""
Created on Fri Jul 17 09:00:59 2020@author: dujidan
"""
import sysfile_fasta = sys.argv[1] #fna 文件
out_path = sys.argv[2] # 输出路径f_fa = open (file_fasta) flag = True
while flag is True:line_fa = f_fa.readline()if line_fa.startswith('>'):#name = line_fa.replace('\n','').replace('>','').replace('|','-').replace('/','_').replace('/','_').replace(' ','_')name = line_fa.replace('>','')f_out = open(out_path + "'" + name.split(' ')[0] + ".fna'",'w')f_out.write(line_fa)else:f_out.write(line_fa)if not line_fa :flag = Falsef_out.close()
f_fa.close()
分割 fasta 文件相关推荐
- 如何提取fasta文件的部分信息
需要将下面的文件整理成三列 即整理成三列,第一列为蛋白ID,第二列为整个一行,第三列是gene ID,中间用制表符分割 #打开文件 f1 = open('X101SC19051138-Z01-bos_ ...
- 使用python读取和分析fasta文件
fasta文件格式 在生物信息学中,FASTA格式(又称为Pearson格式)是一种基于文本的.用于表示核苷酸序列或氨基酸序列的格式.FASTA文件以序列标识和序列作为一个基本单元,每个基本单元分为两 ...
- [笔记]python对FASTA文件的处理
这学期选了生信的选修课-perl/python在生物信息学中的应用 把结课作业的代码整理出来主要是python对FASTA文件的读取和数据处理 FASTA文件数据处理 FASTA文件读取: 只含一个基 ...
- java拆分excel_Jexcel实现按一定规则分割excel文件的方法
本文实例讲述了Jexcel实现按一定规则分割excel文件的方法.分享给大家供大家参考.具体如下: 现有一个excel文档,需要读取它并按照一定的规则,分割之,分割出来的每一段记录需要单独创建一个ex ...
- linux下压缩并分割稀疏文件
稀疏文件是指一个文件中大部分内容都是空字符的文件. 如虚拟机创建的100G的磁盘文件是,但系统不会马上分割出100G的空间,而是仅标识出虚拟硬盘文件是稀疏文件,待到数据存储时,再分配空间,这样可以大大 ...
- linux下使用split 来分割大文件
平常都是使用ssh来进行远程操作和传输文件 , 这个工具可以很方便的实现我们所需要的功能,一直使用正常,没遇到什么问题,今天因为需要传一个412M的文件 , 无论如何也传不过来, 一直系统报错.无奈之 ...
- pupload 文件分块 php,基于Plupload实现Base64分割的文件上传方案
标题:基于Plupload实现Base64分割的文件上传方案 关键词:文件上传.Base64.Plupload.Blob.分割上传 领域:Web前端 作者:孙振强 日期:2018-04-13 目录 背 ...
- linux 分隔大文件,linux系统下分割大文件的方法
本文介绍下在linux系统中,分割大文件的方法,比如一个5gb日志文件切割为很多小块. 在linux中分割大文件,比如一个5gb日志文件,需要把它分成多个小文件,分割后以利于普通的文本编辑器读取. 有 ...
- python分割合并文件
python分割合并文件 python分割文件 将notebook.txt文件分割放到yyy文件夹中: python合并文件 将yyy文件夹中的文件合并到以日期开头的日志文件中: 运用案例:仿真出问题 ...
最新文章
- C#生成唯一的字符串或者数字
- 教你打造一道超级防御的电脑防火墙
- 数控g71编程实例带图_数车编程G71实例讲解
- 电力笔记-30个行业专业词汇(Ⅰ期)
- 直入灵魂的Python教学:《看动漫学Python》让学习不再枯燥
- JavaScript获取图片的原始尺寸
- 柯达新款Wi-Fi数码相机闪亮登场
- 构建项目AppFuse+QuickStart
- 千军万马跑协程goroutine
- 微信聊天小程序——(四、聊天页面)
- javascript常用校验器
- 获取京东JDCookie
- 华为网吧服务器型号,网吧服务器 主机配置
- lte网络是什么意思
- apicloud传递数据
- 人工智能三大学派:符号主义,连接主义,行为主
- 物联卡中心:物联卡能永久使用吗?放手机会不会被封卡
- 如何将EXCEL中同一个工作薄多张工作表合并
- php类的定义与实例化方法
- Rebar3使用介绍(一)基础用法
热门文章
- esp8266--04wifi开关
- 联想服务器系统机房,联想10月1日完成收购IBM x86服务器
- (详细)解决hive报错FAILED: SemanticException Cartesian products are disabled for safety的问题
- Maven:解决IDEA无法下载源码
- logit回归怎么看显著性_SPSS:Logistic回归(Logistic regression)概述
- 倍福CX9020(WINCE 7)使用技巧汇总
- cf卡格式化后还能数据恢复么?相机cf卡误删怎么恢复
- 使用Galen进行Responsive Web测试四部曲
- Three.js 3D 动画场景搭建
- ajax 报错 Content-Type is not allowed by Access-Control-Allow-Headers in preflight respon