def checkCandidate(candidate, results, peptide): rna_complement = genome.transpose_rna(genome.reversecomplement(candidate)) rna = genome.transpose_rna(candidate) if genome.encodeProtein(rna) == peptide: results.append(candidate) elif genome.encodeProtein(rna_complement) == peptide: results.append(candidate)
import re import genome inFile=open('dataset_18_3.txt', 'r') pattern='AUGGCCAUGGCGCCCAGAACUGAGAUCAAUAGUACCCGUAUUAACGGGUGA' pattern=inFile.readline().strip() print genome.encodeProtein(pattern)