python基础操作-Deliaccy

#打印10以内的偶数#
n=0;
while n<=10:
    if n%2==0:
        print(n,end="")
    n += 1

#计算gc含量#
seq='ATGGATGTAGATGATCCCGTCGCAACCCCCTATATAATCGATATATATTGGGGGTTGCGTTCATG' 
GC=(seq.count('G')+seq.count('C'))/len(seq)
print ("GC content is: %.2f%%" % (GC*100))

#生成随机数及其随机数最大值最小值#
import random
lb=[]#建立一个空列表
for i in range(10):
        lb.append(random.randint(0,100))
print(lb)
lb.sort() #随机数排序
print(lb)
print("MAX=",max(lb))
print("MIN=",min(lb))

#函数递归调用斐波那契数列#
def fib(n):
    if n==1:
        return 0
    if n==2:
        return 1
    return fib(n-1)+fib(n-2)#递归调用

print(fib(10))#输出第10个斐波那契数列

#检测orf#
seq='CAATGAAATCTTCAATCTACAGAATATTGATTGGAAAGGTGCGTCAAGCTTCGCCGTTGACACATGTTTGTTGACTGGTATCACAAACAAAGGCTTTGTCACAGACGCATTACATTGAGGGCCTTAACATGGATAACTTCGGTAGCATGATGACAAAATGCTAAGAAAACGCCGCCGGTATGCTCTGCTCTCCATGATATGGACTATCGAACCCATAGTCTAAAGGTACATCAACATTCAACACGCGCCCTGCTCCTTGGGGAATGAAAGGAGGCCAGCTCAGACCGCAGCAGCCACAAATCTCGCCATGTCCTCCACTCGTCCGCCGCCACCCTGGCCTCATGCTCAACCTAACCAAGCAATGCCAGCACCACTGGCCCACGCAAATCGACCCGGACAAGCGCTCCGCATGGCTCGCCCAGCCGACGCACGAGAGCTAATCTTCAAGTTCGGCCTCGAGTTTCTCGAAGCGGGCAGCCC';
for n,n1,n2 in seq:
    while n1,n2,n3=="ATG"
    print(seq)

#dna浓度计算
n=2.02;
if n>=0.1 and n<=1.5:
    
    print("对应DNA浓度=",n*50,"ug/ml");
    
else: print("超出光度计吸光值，无法精确计算");

#ORF（开放阅读框）是指一段ATG开头终止密码子（TAA/TGA/TAG）结尾的基因组片段，一般长度达到50bp以上ORF才是有效的。现在已将某基因片段存储与变量$seq中，请编写一个最简单的ORF预测程序,查找$seq中所有有效ORF（长度超过50bp），并输出到显示屏。#
import re
import sys

dna = """CAATGAAATCTTCAATCTACAGAATATTGATTGGAAAGGTGCGTCAAGCTT
CGCCGTTGACACATGTTTGTTGACTGGTATCACAAACAAAGGCTTTGTCACAGACGCATTACATTGAGGGC
CTTAACATGGATAACTTCGGTAGCATGATGACAAAATGCTAAGAAAACGCCGCCGGTATGCTCTGCTCTCCATGA
TATGGACTATCGAACCCATAGTCTAAAGGTACATCAACATTCAACACGCGCCCTGCTCCTTGGGGAATGAAAGGAGGC
CAGCTCAGACCGCAGCAGCCACAAATCTCGCCATGTCCTCCACTCGTCCGCCGCCACCCTGGCCTCATGCTCAACCTAACCAAG
CAATGCCAGCACCACTGGCCCACGCAAATCGACCCGGACAAGCGCTCCGCATGGCTCGCCCAGCCGACGCACGAGAGC
TAATCTTCAAGTTCGGCCTCGAGTTTCTCGAAGCGGGCAGCCC"""
def find_orf(dna):
    ret = []
    dna_len = int(len(dna) / 3) * 3
    i = 0
    for i in range(0, dna_len, 3):
        codon = dna[i:i + 3]
        if codon != 'ATG': continue
        base = codon
        for j in range(i + 3, dna_len, 3):
            codon = dna[j:j + 3]
            base += codon
            if ((codon == 'TAA') | (codon == 'TAG') | (codon == 'TGA')):
                ret.append(base)
                continue
    return ret
for i in range(0,len(find_orf(dna))):
    print(len(find_orf(dna)[i]))
print(find_orf(dna))