def map_peptide(target_up_id, peptide, pos):
seq = uniprot_client.get_sequence(target_up_id)
peptide_start = seq.find(peptide)
# Peptide not found in sequence, return None
if peptide_start == -1:
return None
target_pos = peptide_start + pos
return target_pos
def test_get_sequence():
seq = uniprot_client.get_sequence('P00533')
assert len(seq) > 1000
def motif_from_position(up_id, pos, window=7):
seq = uniprot_client.get_sequence(up_id)
return ProtMapper.motif_from_position_seq(seq, pos, window)