Please help me to rectify the error in mod9v2
My script
from modeller import * log.verbose() env = environ() #-- Prepare the input files #-- Read in the sequence database sdb = sequence_db(env) #sdb.read(seq_database_file='C:\Program Files\Modeller9v2\modlib\pdball.pir', seq_database_format='PIR',chains_list='ALL', minmax_db_seq_len=(30, 4000), clean_sequences=True) #-- Write the sequence database in binary form #sdb.write(seq_database_file='C:\Program Files\Modeller9v2\modlib\pdball.bin', seq_database_format='BINARY', chains_list='ALL') #-- Now, read in the binary database sdb.read(seq_database_file='C:\Program Files\Modeller9v2\modlib\pdball.bin', seq_database_format='BINARY', chains_list='ALL') #-- Read in the target sequence/alignment aln = alignment(env) aln.append(file='p17787.fas', alignment_format='FASTA', align_codes='p17787') #-- Convert the input sequence/alignment into # profile format prf = aln.to_profile() #-- Scan sequence database to pick up homologous sequences prf.build(sdb, matrix_offset=-450, rr_file='${LIB}/blosum62.sim.mat', gap_penalties_1d=(-500, -50), n_prof_iterations=1, check_profile=False, max_aln_evalue=0.01) #-- Write out the profile in text format prf.write(file='build_profile.prf', profile_format='TEXT') #-- Convert the profile back to alignment format aln = prf.to_alignment() #-- Write out the alignment file aln.write(file='build_profile.ali', alignment_format='PIR') aln.check()
error what i have got
'import site' failed; use -v for traceback Traceback (most recent call last): File "build_profile.py", line 21, in ? aln.append(file='p17787.fas', alignment_format='FASTA', align_codes='p17787' ) File "C:\Program Files\Modeller9v2\modlib\modeller\alignment.py", line 74, in append rewind_file, close_file) _modeller.error: read_al_373E> Protein specified in ALIGN_CODES(i) was not found in the alignment file; ALIGN_CODES( 1) = p17787
____________________________________________________________________________________ Be a better Globetrotter. Get better travel answers from someone who knows. Yahoo! Answers - Check it out. http://answers.yahoo.com/dir/?link=list&sid=396545469