Minsu Kim/python/Rosalind

Problem

basic problem

[1]   [2]   [3]

related with bioinformatics

[1]   [2]   [3]   [4]   [5]   [6]

[7]   [8]   [9]   [10]   [11]   [12]

Computing CG Content

f = open('rosalind_gc.txt','r')
seq_list = dict()
for line in f_fa:
    if line.startswith('>'):
        name = line.strip('>').strip()
        seq_list[name] = ''
    else:
        seq_list[name]= line.strip()
f.close()

seq_CG = dict()
for k in seq_list:
    count = 0
    s = seq_list[k]
    s = list(s)
    for i in range(0,len(s)):
        if s[i] == 'C' or s[i] == 'G':
            count = count + 1
    p = count / len(s)
    seq_CG[k] = p

inverse = [(value, key) for key, value in seq_CG.items()]
max(inverse)

import matplotlib.pyplot as plt

seq_list = []
f_fa = open('Ecoli_genome.fasta', 'r')
for line in f_fa:
    if not line.startswith('>'):
        seq_list.append(line.strip())
tmp_seq = ''.join(seq_list)

count_G = 0
count_C = 0
count_list = []
for tmp_n in tmp_seq:
    if tmp_n == 'G':
        count_G += 1
    if tmp_n == 'C':
        count_C += 1
    count_list.append(count_G-count_C)

pos_list = range(0,len(tmp_seq))

fig = plt.figure(figsize = (10,6))
ax1 = fig.add_subplot(1,1,1)
ax1.plot(pos_list, count_list, 'b-')
ax1.grid()
plt.show()
f_fa.close()

 

Dynamic Program

def maxvalue(a,b,c):
    a=a+5
    b=b-6
    c=c-6
    return max(a,b,c)
def maxvalue2(a,b,c):
    a=a-2
    b=b-6
    c=c-6
    return max(a,b,c)

def compare(a,b):
    import numpy as np
    scoring = np.zeros([len(a)+1,len(b)+1])
    for i in range(0,len(a)+1):
        scoring[i,0]=-6*i
    for j in range(0,len(b)+1):
        scoring[0,j]=-6*j

    for row in range(1,len(a)+1):
        for column in range(1,len(b)+1):
            if b[column-1]==a[row-1]:
                scoring[row,column] = maxvalue(scoring[row-1,column-1],scoring[row-1,column],scoring[row,column-1])
            if b[column-1]!=a[row-1]:
                scoring[row,column] = maxvalue2(scoring[row-1,column-1],scoring[row-1,column],scoring[row,column-1])
    print(scoring)

    list=[]
    for row in range(1,len(a)+1):
        for column in range(1,len(b)+1):
            new=scoring[row,column]
            list.append(new)
    mxv=max(list)
    print(mxv)