2016-10-29 21:38 GMT+02:00 :
>
> Code:
>
> [...]
>
> for i in range (len(protein) & len(seq1)) :
>
> if protein[i] != mutantPRO[i] :
>print (protein[i] + str(i) + mutantPRO[i])
>A+= 1
> else:
> if seq1[i:i+3] != mutant[i:i+3]:
>
Disha,
Before you struggle to reinvent the wheel, you might want to check out the
Biopython package.
http://biopython.org/wiki/Biopython
I haven't used it for a few years, but the version that I did use was very
comprehensive.
--
https://mail.python.org/mailman/listinfo/python-list
On 2016-10-29 20:38, dishaachary...@gmail.com wrote:
Code:
A = 0
B= 0
i=0
j=0
# opening the files
infile1 = open("CDSsrebf1.txt")
infile2 = open("PROsrebf1.txt")
infile3 = open("mutant.txt")
print(" 1st line of WT SREBF1 (CDS):",infile1.readline())
print ("1st line of mutant protein of SREBF1: "
Code:
A = 0
B= 0
i=0
j=0
# opening the files
infile1 = open("CDSsrebf1.txt")
infile2 = open("PROsrebf1.txt")
infile3 = open("mutant.txt")
print(" 1st line of WT SREBF1 (CDS):",infile1.readline())
print ("1st line of mutant protein of SREBF1: ", infile3.readline())
print ("1st line of protein of S
