python学习——读取染色体长度(七:for循环对染色体序列进行反向互补)
对fasta文件genome_test.fa中的染色体序列进行反向互补,并输出到文件genome_test_RC.fa
genome_test.fa
>chr1
ATATATATAT
>chr2
ATATATATATCGCGCGCGCG
>chr3
ATATATATATCGCGCGCGCGATATATATAT
>chr4
ATATATATATCGCGCGCGCGATATATATATCGCGCGCGCG
>chr5
ATATATATATCGCGCGCGCGATATATATATCGCGCGCGCGATATATATAT
新疆Reverse_Complement.py文件,并输入如下python脚本
python脚本
1 # import os # 导入模块os 2 import sys # 导入模块sys 3 f_fasta = sys.argv[1] # 从命令行获取文件名 4 f = open(f_fasta) # 打开文件 5 f_RC = open("genome_test_RC.fa","w+") 6 # 逐行读取 7 lines = f.readlines() 8 for line in lines: 9 line = line.strip() # 去掉行尾的换行符 10 if (line.startswith(">")): 11 chr_id = line + '_RC' 12 else: 13 chr_seq = line[::-1].replace('A','t').replace('T','a').replace('C','g').replace('G','c').upper() 14 # 输出结果 15 print(chr_id) 16 print(chr_seq) 17 18 f_RC.write(chr_id + '\n') 19 f_RC.write(chr_seq + '\n') 20 f.close() 21 f_RC.close()
从cmd终端命令行输入参数,调用上述python脚本,并对genome_test.fa进行处理
1 E:\15_python\DEBUG>python Reverse_Complement.py genome_test.fa
结果
genome_test_RC.fa
>chr1_RC
ATATATATAT
>chr2_RC
CGCGCGCGCGATATATATAT
>chr3_RC
ATATATATATCGCGCGCGCGATATATATAT
>chr4_RC
CGCGCGCGCGATATATATATCGCGCGCGCGATATATATAT
>chr5_RC
ATATATATATCGCGCGCGCGATATATATATCGCGCGCGCGATATATATAT
转载于:https://www.cnblogs.com/caicai2019/p/10799760.html
python学习——读取染色体长度(七:for循环对染色体序列进行反向互补)相关推荐
- python学习[第十三篇] 条件和循环
python学习[第十三篇] 条件和循环 if语句 单一if 语句 if语句有三个部分构成,关键字if本身,判断结果真假的条件表达式,以及表达式为真或非0是执行的代码 if expression: e ...
- Python学习笔记D2(条件与循环)
Python学习笔记D2(条件与循环) 条件语句 if expression:(可以利用布尔操作符and,or,not) *****语句- elif:(else if) ***** 语句- else: ...
- Python学习笔记:第七站 夫妻站
Python学习笔记 文章目录 Python学习笔记 第七站 夫妻站 1. 什么是字典 2. 字典的创建 3. 字典的查询操作 4. 字典元素的增.删.改操作 5. 字典推导式 6. 本章作业 课程笔 ...
- Python学习笔记三之编程练习:循环、迭代器与函数
Python学习笔记三之编程练习 1. 编程第一步 # 求解斐波纳契数列 #/user/bin/python3#Fibonacci series:斐波那契数列 #两个元素的总和确定了下一个数 a,b= ...
- 王亟亟的Python学习之路(七)-date,continue,迭代对象,生成器
转载请注明出处:王亟亟的大牛之路 大多数的语言都提供了时间操作的相应类,诸如java的java.util.Date,java.util.Calendar等等, Python中是time 和 calen ...
- 王亟亟的Python学习之路(四)-循环,条件,Range,list和tuple
转载请注明出处:王亟亟的大牛之路 这是本次Python学习历程的第4篇,之前的文章:传送门 这篇主要讲以下几个知识点 - 循环 - 条件语句 - range函数 - list - tuple 循环 P ...
- [Python学习] 专题二.条件语句和循环语句的基础知识
前面讲述了"专题一.函数的基础知识",而这篇文章讲述的Python的条件语句和循环语句的基础知识.主要内容包括: 1.条件语句:包括单分支.双分支和 ...
- python学习笔记(四)、条件、循环及其他语句
1 再谈print和import 1.1 打印多个参数 print 能够同时打印多个表达式,并且能自定义分隔符.如下: print('a','b','c') --> a b c print('a ...
- Python学习:命令行运行,循环结构
一.安装配置和运行方法 1.安装OpenCV 3.1: 假设安装目录为"C:\Python34" 2.配置环境变量: 方法一:直接配置:打开"控制面板",搜索& ...
最新文章
- oracle一条sql的执行过程,请问执行一条sql的整个过程是怎样的,谢谢!
- 这家自动驾驶公司融资2.7亿美元,京东美团齐出手
- 2014届华为校园招聘机试题
- Redis中的Sentinel 配置
- css3中的 @Keyframes
- 34. 二分查找左右边界
- python输出函数使用_python基本输入输出函数
- 特斯拉上海工厂启动试生产 首批车被曝将销售给工厂员工
- apk逆向思路_安卓逆向和手游辅助学习路线
- java需要下载哪些,全网最新
- 页面三个txt加载联动省市县的代码,类似淘宝的收货地址的布局
- 联想小新电脑摄像头黑屏、检测不到设备、指示灯不亮解决方案
- php 京东 联盟 链接,京东联盟推广链接在微信上如何直接跳转接入到京东小程序下单购买?...
- python提取一句话中的数字_从Python的字符串列表中提取数字
- Google的C++编程规范总结
- MySQL-存储表情字符
- android分屏后布局变化问题,Android分屏多窗口实践
- Linux关闭系统更新
- 七周七并发模型与七周七语言
- [C++]<numeric>头文件介绍