Changes
no edit summary
<h4>I solve the problem relative related to the bioinformatics in web site : rosalind</h4>
<ph1> Problem</ph1>
<h4><a href="http://biolecture.org/index.php?title=Minsu_Kim/python/Rosalind / function"p>basic function</a>problem</h4p>
<h4p><a href="http://biolecture.org/index.php?title=Minsu_Kim/python/Rosalind / b1">[1]</a> <a href="http://biolecture.org/index.php?title=Minsu_Kim/python/Rosalind / b2">[12]</a> <a href="http://biolecture.org/index.php?title=Minsu_Kim/python/Rosalind / b3">[3]</a></h4p>
<h2p>functionrelated with bioinformatics</h2p>
<p>def patterncount(text, pattern)<a href="http://biolecture.org/index.php?title=Minsu_Kim/python/Rosalind / problem1">[1]<br /a> count<a href=0"http://biolecture.org/index.php?title=Minsu_Kim/python/Rosalind / problem2">[2]<br /a> for i in range(0,len(text))<a href="http://biolecture.org/index.php?title=Minsu_Kim/python/Rosalind / problem3">[3]<br /a> <a href="http://biolecture.org/index.php?title=Minsu_Kim/python/Rosalind /problem4">[4]</a> if(text<a href="http://biolecture.org/index.php?title=Minsu_Kim/python/Rosalind /problem5">[i:len(pattern)+i5]==pattern):<br /a> count +<a href="http://biolecture.org/index.php?title= 1Minsu_Kim/python/Rosalind /problem6">[6]<br /a> return count</p>
<p>def complementing(sequence):<br /> sequence a href= sequence["http::-1]<br /> text /biolecture.org/index.php?title= list(sequence)<br Minsu_Kim/python/> for i in range(0,len(text)):<br Rosalind /problem7"> if text[i7] == 'A':<br /a> text[i] = 'T'<br /> elif text[i] =a href= 'G'"http:<br /> text[i] /biolecture.org/index.php?title= 'C'<br Minsu_Kim/python/Rosalind /problem8"> elif text[i8] == 'C':<br /a> text[i] = 'G'<br /> elif text[i] a href== 'T'"http:<br /> text[i] = 'A'<br /> <br biolecture.org/> s index.php?title= ''.join(text)<br Minsu_Kim/> print(s)<br python/> return 0<Rosalind /pproblem9"> <h1>rosalind transcript basic[9]</h1> <pa>f = open('nbsp;rosalind_rna.txt')<br />a href= f"http://biolecture.readline()<br org/>findex.close<php?title=Minsu_Kim/p> <p>text = list(a)<python/p> <p>for i in range(0,len(text)):<br Rosalind /problem10"> if text[i10] == 'T':<br /a> text[i] = 'U'<br />a href= ''.join(text)</p> <p>print(a)<"http:/p> <p>f = open('rosalind_revc2.txt','r')<br />b = fbiolecture.readline()<br org/>findex.close</p> <p>b php?title= b[::-1]<br Minsu_Kim/>text = list(b)<python/p> <p>for i in range(0,len(text)):<br Rosalind /problem11"> if text[i11] == 'A':<br /a> text[i] = 'T'<br /> elif text[i] =a href= 'G'"http:<br /> text[i] = 'C'<br /> elif text[i] biolecture.org/index.php?title== 'C':<br Minsu_Kim/> text[i] = 'G'<br python/> elif text[i] == 'T':<br Rosalind /problem12"> text[i12] = 'A'<br /a> <br />b = ''.join(text)</p> <p>print(b)</p>
<h2>Computing CG Content</h2>
<p>f_fa f = open('rosalind_gc.txt','r')<br />
seq_list = dict()<br />
for line in f_fa:<br />
else:<br />
seq_list[name]= line.strip()<br />
<p>seq_CG = dict()<br />
<p> </p>
<h2>Dynamic Program</h2>
<p>def maxvalue(a,b,c):<br />
a=a+5<br />
b=b-6<br />
c=c-6<br />
return max(a,b,c)<br />
def maxvalue2(a,b,c):<br />
a=a-2<br />
b=b-6<br />
c=c-6<br />
return max(a,b,c)</p>
<p>def compare(a,b):<br />
import numpy as np<br />
scoring = np.zeros([len(a)+1,len(b)+1])<br />
for i in range(0,len(a)+1):<br />
scoring[i,0]=-6*i<br />
for j in range(0,len(b)+1):<br />
scoring[0,j]=-6*j</p>
<p> for row in range(1,len(a)+1):<br />
for column in range(1,len(b)+1):<br />
if b[column-1]==a[row-1]:<br />
scoring[row,column] = maxvalue(scoring[row-1,column-1],scoring[row-1,column],scoring[row,column-1])<br />
if b[column-1]!=a[row-1]:<br />
scoring[row,column] = maxvalue2(scoring[row-1,column-1],scoring[row-1,column],scoring[row,column-1])<br />
print(scoring)</p>
<p> list=[]<br />
for row in range(1,len(a)+1):<br />
for column in range(1,len(b)+1):<br />
new=scoring[row,column]<br />
list.append(new)<br />
mxv=max(list)<br />
print(mxv)</p>
<p><img alt="" src="/ckfinder/userfiles/images/bandicam%202017-12-12%2019-14-32-803.jpg" style="height:530px; width:640px" /></p>