I'm receiving the following error on linux with Modeller 9.2: _modeller.error: read_al_373E> Protein specified in ALIGN_CODES(i) was not found in the alignment file; ALIGN_CODES( 1) = mod2hg42
Here's the alignment line from my python file: aln = alignment(env, file='mod2hg42.fa', align_codes=('mod2hg42'))
Here's the contents of mod2hg42.fa (located in the same directory): >mod2hg42 EEKLRRYLKRTVTELDSVTARLREVEHRAGEPIAIVGXACRFPGDVDSPESFWEFVSGGGDAIAEAPADRGWEPD PDARLGGXLAAAGDFDAGFFGISPREALNXDPQQRIXLEISWEALERAGHDPVSLRGSATGVYTGSSTNDYGFLA VSDPDEAHPYAITGTASSVASGRVAYCLGLEGPAXTVDTACSSSLTALHLAXESLRRDECGLALAGGVNVLSTPM VTLGFDSQGGLAADGRIKPFSKAADGYTRSEGAGVLVLQRLSAARREGRPVLAVLRGSAVNQDGASNGLIAPSGP AQQRVIRRALENAGVRAGDVDYVEAHGTGTRLGDPIEVHALLSTYGAERDPDDPLWIGAVKTNVGHLESAAGVAG VXKAVLALRHGEXPRTLHFDEPSPQIEWAVSVVSQARSWPAGERPRRAGVSSFGFGGTNAHVIVEEAPEADEPV
I've tried dozens of combinations, but I can't get around this error. Any advice?
Dave Roe wrote: > I'm receiving the following error on linux with Modeller 9.2: > _modeller.error: read_al_373E> Protein specified in ALIGN_CODES(i) was not found in the alignment file; ALIGN_CODES( 1) = mod2hg42 > > Here's the alignment line from my python file: > aln = alignment(env, file='mod2hg42.fa', align_codes=('mod2hg42'))
Your alignment file looks like a FASTA file, and you haven't specified 'alignment_format' when reading it in, so Modeller will treat it as a PIR file by default. Since there is no valid PIR sequence in that file, it won't find your sequence.
Add alignment_format='FASTA' when reading your alignment, or (better) switch your alignment to PIR format, and everything should work.
Ben Webb, Modeller Caretaker
participants (2)
-
Dave Roe -
Modeller Caretaker