擅长:python、mysql、java
<p>以下是我对代码的建议,非常直截了当,简短易懂,除了<code>re</code>之外没有任何其他导入包:</p>
<pre><code>import re
def findGaps(aSeq):
# Get and print the list of gaps present into the sequence
gaps = re.findall('[-]+', aSeq)
print('Number of gaps = {0} \n'.format(len(gaps)))
# Get and print start index, end index and length for each gap
for i,gap in enumerate(gaps,1):
startIndex = aSeq.index(gap)
endIndex = startIndex + len(gap) - 1
print('Index Position of Gap region {0} = {1} to {2}'.format(i, startIndex, endIndex))
print('Length of Gap region {0} = {1} \n'.format(i, len(gap)))
aSeq = aSeq.replace(gap,'*' * len(gap), 1)
findGaps("ATC----GCTGTA--A-----T")
</code></pre>