Witam!
jako ze nie jestem za dobry w programowaniu, a pythona zoboaczylem na oczy 2 dni temu, mam mala prozbe.
Czy moglby mi ktos pomoc w zadaniu z bioinformatyki?
Mam program ktory konwertuje zbiory RSF na format FASTA.
Kod programu:
<div class='codetop'>CODE</div><div class='codemain' style='height:200px;white-space
re;overflow:auto'>#!/usr/bin/python
import re
ENTRIES = re.compile('name\s+?(\S+).+?sequence\s+?([^\}]+)',re.S)
FASTA_ALN = True
def convert(options,linewidth=60):
try: f = open(options.inputfile,'r')
except IOError: raise " File '%s' not found." % options.inputfile
data = f.read()
f.close()
entries = re.findall(ENTRIES,data)
f = open(options.outputfile,'w')
for name, sequence in entries:
f.write('>%s\n' % name)
if FASTA_ALN:
cleaned = re.subn('\s','',sequence)[0].upper()
cleaned = re.subn('\.','-',cleaned)[0].upper() + '\n'
else: cleaned = re.subn('[\s\.]','',sequence)[0].upper()
f.write('\n'.join([cleaned[i:i+linewidth] \
for i in range(0,len(sequence),linewidth)]))
f.close()
def main():
from optparse import OptionParser
usage = """ %prog [options]
program for convertion biological sequence data from RST to Fasta format """
parser = OptionParser(usage=usage)
parser.add_option('-i', '--input', help='input RST filename',
action='store', type='string', dest='inputfile')
parser.add_option('-o', '--output', help='output FASTA filename',
action='store', type='string', dest='outputfile',
default='sequence.fasta')
options, args = parser.parse_args()
if not options.inputfile:
parser.error(" You MUST provide input filename. "+\
"Run program with '-h' for help.")
convert(options)
if __name__ == '__main__': main()
</div>
a teraz zadania z tym zwiazene:
-jakie moduły ze standardowej biblioteki Pythona są wykorzystane w programie: czy jest to modol "re" i "optparse"
-Wyjaśnić działanie ostatniej konstrukcji warunkowej.
-Napisać funkcję readFastaAlignment, która będzie wczytywać plik z dopasowaniem w formacie FASTA a następnie będzie zwracać w postaci listy, której każdym z elementów będzie dwuelementowa krotka zawierająca opis sekwencji oraz samą sekwencje.
- Napisać funkcję, która jako argument będzie pobierać listę zwracaną przez funkcję readFastaAlignment i w kolejnych liniach będzie drukować oznaczenie sekwencji oraz jej długość.
Z gory dziekuje za pomoc
jako ze nie jestem za dobry w programowaniu, a pythona zoboaczylem na oczy 2 dni temu, mam mala prozbe.
Czy moglby mi ktos pomoc w zadaniu z bioinformatyki?
Mam program ktory konwertuje zbiory RSF na format FASTA.
Kod programu:
<div class='codetop'>CODE</div><div class='codemain' style='height:200px;white-space
re;overflow:auto'>#!/usr/bin/pythonimport re
ENTRIES = re.compile('name\s+?(\S+).+?sequence\s+?([^\}]+)',re.S)
FASTA_ALN = True
def convert(options,linewidth=60):
try: f = open(options.inputfile,'r')
except IOError: raise " File '%s' not found." % options.inputfile
data = f.read()
f.close()
entries = re.findall(ENTRIES,data)
f = open(options.outputfile,'w')
for name, sequence in entries:
f.write('>%s\n' % name)
if FASTA_ALN:
cleaned = re.subn('\s','',sequence)[0].upper()
cleaned = re.subn('\.','-',cleaned)[0].upper() + '\n'
else: cleaned = re.subn('[\s\.]','',sequence)[0].upper()
f.write('\n'.join([cleaned[i:i+linewidth] \
for i in range(0,len(sequence),linewidth)]))
f.close()
def main():
from optparse import OptionParser
usage = """ %prog [options]
program for convertion biological sequence data from RST to Fasta format """
parser = OptionParser(usage=usage)
parser.add_option('-i', '--input', help='input RST filename',
action='store', type='string', dest='inputfile')
parser.add_option('-o', '--output', help='output FASTA filename',
action='store', type='string', dest='outputfile',
default='sequence.fasta')
options, args = parser.parse_args()
if not options.inputfile:
parser.error(" You MUST provide input filename. "+\
"Run program with '-h' for help.")
convert(options)
if __name__ == '__main__': main()
</div>
a teraz zadania z tym zwiazene:
-jakie moduły ze standardowej biblioteki Pythona są wykorzystane w programie: czy jest to modol "re" i "optparse"
-Wyjaśnić działanie ostatniej konstrukcji warunkowej.
-Napisać funkcję readFastaAlignment, która będzie wczytywać plik z dopasowaniem w formacie FASTA a następnie będzie zwracać w postaci listy, której każdym z elementów będzie dwuelementowa krotka zawierająca opis sekwencji oraz samą sekwencje.
- Napisać funkcję, która jako argument będzie pobierać listę zwracaną przez funkcję readFastaAlignment i w kolejnych liniach będzie drukować oznaczenie sekwencji oraz jej długość.
Z gory dziekuje za pomoc