python学习——读取染色体长度(四:获取最长染色体的编号)

# 读取fasta

# 解析每条序列的长度

chr_len = {‘chr1‘:10,‘chr2‘:20,‘chr3‘:30,‘chr4‘:40,‘chr5‘:15}

# 求和

total_len = sum(chr_len.values())

# 求最长染色体编号

max_chr =‘‘

max_len = 0

for [chr,len] in chr_len.items():

if len > max_len:

print(‘current chr > max chr‘)

max_chr = chr

max_len = len

elif len == max_len:

print(‘current chr == max chr‘)

else:

print(‘current chr < max chr‘)

# 输出结果

print(total_len)

print(max_len)

print(max_chr)

原文地址:https://www.cnblogs.com/caicai2019/p/10784047.html

时间: 2024-10-15 13:18:40

python学习——读取染色体长度(四:获取最长染色体的编号)的相关文章

python学习之最简单的获取本机ip信息的小程序

文章是从我的个人博客粘贴过来的,大家可以直接访问我的个人博客哦 http://www.iwangzheng.com 获取本机ip信息的命令ifconfig总是在用,这次拿到pyhton代码里,感觉python好麻烦啊,每次都要import. 可以在程序里直接read出来,也可以先在命令行执行ifconfig > xx.txt,这样可以直接访问倒入后的文本信息. 这个小程序有个问题,与朋友讨论了还没解决彻底,如果用findall是可以匹配所有的,但是怎么解决用match只匹配一次的问题呢,当然如果

王亟亟的Python学习之路(四)-循环,条件,Range,list和tuple

转载请注明出处:王亟亟的大牛之路 这是本次Python学习历程的第4篇,之前的文章:传送门 这篇主要讲以下几个知识点 - 循环 - 条件语句 - range函数 - list - tuple 循环 Python中的循环语句有 for 和 while.同样需要注意冒号和缩进.另外,在Python中没有do..while循环. 需要注意冒号和缩进,因为没有{}来帮助我们把循环里面的内容封起来 while循环 while 判断逻辑: 业务内容 或者for-in循环 names = ['12',12,

python学习——读取染色体长度(三、用循环或者函数求总长并获取最长染色体长度)

# 读取fasta # 解析每条序列的长度 chr_len = [10,20,30,40,50] # 求和 # 方法一:通过循环 total_len = 0 #定义total_len的初始长度 for len in chr_len: # 从列表chr_len中每次取一个值交给len total_len += len # 或者total_len = total_len + len # 方法二:通过函数sum() total_len = sum(chr_len) # 求最长染色体编号 max_len

python学习——读取染色体长度(七:读取fasta文件)

读取fasta文件genome_test.fa,并计算染色体总长,同时输出最长染色体编号.序列以及长度 fasta文件genom_test.fa的内容如下: >chr1ATATATATAT>chr2ATATATATATCGCGCGCGCG>chr3ATATATATATCGCGCGCGCGATATATATAT>chr4ATATATATATCGCGCGCGCGATATATATATCGCGCGCGCG>chr5ATATATATATCGCGCGCGCGATATATATATCGCGCG

python学习——读取染色体长度(二、向前一步:通过染色体序列获得长度信息)

# 读取fasta # 解析每条序列的长度 chr1_seq = 'ATATATATAT' chr2_seq = 'ATATATATATCGCGCGCGCG' chr3_seq = 'ATATATATATCGCGCGCGCGATATATATAT' chr4_seq = 'ATATATATATCGCGCGCGCGATATATATATCGCGCGCGCG' chr5_seq = 'ATATATATATCGCGCGCGCGATATATATATCGCGCGCGCGATATATATAT' # 求和 tot

python学习笔记(基础四:模块初识、pyc和PyCodeObject是什么)

一.模块初识(一) 模块,也叫库.库有标准库第三方库. 注意事项:文件名不能和导入的模块名相同 1. sys模块 import sys print(sys.path) #打印环境变量 print(sys.argv) #打印相对路径,在pycharm输出结果是绝对路径,因为在pycharm中调用 print(sys.argv[2]) #可以从数据列表中取值,[2]指取第三位.        标准库存放位置:C:\Python35\Lib 第三方库存放位置C:\Python35\Lib\site-p

python学习之路(四)

继续昨天的学习,学到了数组. 首先有两个数组,name1和name2.我们可以将两个数组合并 name1=[1,2,3,4] name2=[5,6,7,8] names=name1.extend(name2) print(names) 最后输出的是[1,2,3,4,5,6,7,8],这就是在name1后面接上name2 然后来到了一个小重点,copy(复制) 这里分为深copy和浅copy,我在学习之路(二)中说过一个复制的问题,就是在赋值的时候,值是存在内存里面的,而不是随着变量变化的,这里的

python学习之第十四天补充

本节内容 学员管理系统练习 Django ORM操作进阶 用户认证 Django练习小项目:学员管理系统设计开发 带着项目需求学习是最有趣和效率最高的,今天就来基于下面的需求来继续学习Django 项目需求: 1.分讲师\学员\课程顾问角色,2.学员可以属于多个班级,学员成绩按课程分别统计3.每个班级至少包含一个或多个讲师4.一个学员要有状态转化的过程 ,比如未报名前,报名后,毕业老学员5.客户要有咨询纪录, 后续的定期跟踪纪录也要保存6.每个学员的所有上课出勤情况\学习成绩都要保存7.学校可以

python学习第七十四天:单表查询

单表查询 创建表 创建模型 在models.py中定义类,一个类就是一张表 from django.db import models class Book(models.Model): id = models.AutoField(primary_key=True) name = models.CharField(max_length=64) pub_data = models.DateField() price = models.DecimalField(max_digits=5, decima