<p>所以这个python代码是给我的。这意味着它的工作,但显然我有问题。我只是想知道有没有人能告诉我出了什么问题?你知道吗</p>
<p>当我试着跑的时候,错误总是发生</p>
<pre><code> intron1_length = ((my_dna.find(‘gtcata’)) –
(my_dna.find(‘atcgat’)))
</code></pre>
<p>错误为“-”。只是不知道应该去哪里。谢谢!你知道吗</p>
<pre><code> # -*- coding: utf-8 -*-
# We’ll start with “data-aware” way to do this first.
# Since we have the sequence, we can see that exon 1 ends with
# ‘GACTA’. We can use that information, and that method, for
# finding the positions of the elements of this gDNA. This is NOT
# an elegant solution, and would not work if we didn’t know the
# DNA sequence ahead of time.
my_dna = "ATCGATCGATGGTCGAATGACTAgtcatagctatgcatgtagctactc
gatcgtattttattcgatcgatcgatCGATCGATCATGCTATCATCGATCGATATCGATGCATC
GACTACTATgtcatggctatgcatcgatcgtattttattcgatcgttcgatGATCGATCGATCGACTGACTTTGAA"
# here we introduce another useful operator for strings: len
gene_length = len(my_dna)
# we’ll use the starting position of useful substrings in the
# sequence to find the positions of the exon-intron boundaries.
# We’ll then use those to find the length of each segment.
exon1_length = (my_dna.find('Agtcata'))
intron1_length =((my_dna.find('gtcata')) - (my_dna.find('atcgat')))
exon2_length = ((my_dna.find('CGATCG'))- (my_dna.find('Tgtcatg')))
intron2_length = ((my_dna.find('gtcatg'))- (my_dna.find('tGATCGA')))
exon3_length = (gene_length(my_dna.find('GATCGA'))
print ("Gene length:" + str(gene_length))
print ("Exon1 length:" + str(exon1_length))
print ("Intron1 length:" + str(intron1_length))
print ("Exon2 length:" + str(exon2_length))
print ("Intron2 length:" + str(intron2_length))
print ("Exon2 length:" + str(exon3_length))
</code></pre>