随笔分类 -  Python_script

如何得到给定序列的互补序列以及反向互补序列
摘要:用 python 实现如下: 1 #!/usr/bin/python 2 # Complementing DNA 3 4 my_dna = "ACTGATCGATTACGTATAGTATTTGCTATCATACATATATATCGATGCGTTCAT" 5 # 由于python 区分大小写,所以先... 阅读全文

posted @ 2015-09-17 20:39 OA_maque 阅读(3791) 评论(0) 推荐(0)

将给定序列翻译成蛋白质序列
摘要:利用 dictionary 可以将给定的cDNA序列翻译成蛋白序列 1 #!/bin/python 2 # Dictionary protein translation 3 4 my_dna = open("/home/maque/my_dna.txt").read().replace('\n',... 阅读全文

posted @ 2015-09-17 20:11 OA_maque 阅读(1032) 评论(0) 推荐(0)

统计给定序列的AT-content ,并设置小数点位数
摘要:习题来源: Python for Biologists: A complete programming course for beginner 1 #!/bin/python 2 # calculate the AT content of a DNA seq 3 4 def get_at_cont... 阅读全文

posted @ 2015-09-17 19:18 OA_maque 阅读(274) 评论(0) 推荐(0)

如何在一序列中寻找E-BOX 元件,并将其变为小写
摘要:给定一个序列,寻找所有的E-BOX motif (CNNTTG) , 并且将其变为小写, 用 python 实现如下: 1 #!/bin/python 2 # Date: 2015.8.01 3 # Author: 4 """ Search E-box motif(CNNTTG) in a DNA... 阅读全文

posted @ 2015-09-17 18:42 OA_maque 阅读(874) 评论(0) 推荐(0)

列出给出序列的crispr 引物
摘要:手动寻找cripsr 引物比较麻烦,而现在有些网站可以完成这一任务,但是,用python 去实现它也很简单。以下是脚本: 1 #!/usr/bin/python 2 # list all crispr-target(20 bp + NGG) 3 4 import re 5 from Bio.Seq... 阅读全文

posted @ 2015-09-14 23:19 OA_maque 阅读(685) 评论(0) 推荐(0)

统计任意两个正整数之间所有奇数的和
摘要:题目来源: Conditions and Loops 1 #!/usr/bin/python3 2 # sum all odds between two positive integer 3 4 a = int(input("Please input one positive integer > ... 阅读全文

posted @ 2015-09-12 11:10 OA_maque 阅读(674) 评论(0) 推荐(0)

去掉文件中重复的元件, 生成无重复的列表
摘要:有一个文件,内容如下 file_name1 file_name4 file_name3 file_name11 file_name7 file_name1 file_name3 file_name1 file_name8 文件里有重复的名字,现想去掉重复的,只保留一个,用 python 脚本实现,代... 阅读全文

posted @ 2015-09-11 15:43 OA_maque 阅读(284) 评论(0) 推荐(0)

导航