def test_peptide_round_trip():
pm = ProtMapper()
pos = '187'
motif, site_pos = pm.motif_from_position('P28482', '187')
ms = pm.map_peptide_to_human_ref('MAPK1', 'hgnc', motif, site_pos)
assert isinstance(ms, MappedSite)
assert ms == MappedSite(up_id='P28482', error_code=None, valid=True,
orig_res=None, orig_pos=None, mapped_id='P28482',
mapped_res='Y', mapped_pos='187',
description=None, gene_name='MAPK1')
def test_motif_from_position():
sm = ProtMapper()
# Try with string position
(motif, respos) = sm.motif_from_position('O95218-2', '120', window=7)
assert motif == 'EYIEREESDGEYDEF'
assert respos == 8
# Also try with int
motif, respos = sm.motif_from_position('O95218-2', '120')
assert motif == 'EYIEREESDGEYDEF'
assert respos == 8
# Change the window size
motif, respos = sm.motif_from_position('O95218-2', '120', window=5)
assert motif == 'IEREESDGEYD'
assert respos == 6
# Try with a residue that is close to the end of the protein
motif, respos = sm.motif_from_position('Q04637-1', '1596', window=7)
assert motif == 'LREAEEESDHN'
assert respos == 8
# Try with a residue that is close to the start of the protein
motif, respos = sm.motif_from_position('Q04637-1', '3', window=7)
assert motif == 'MNKAPQSTGP'
assert respos == 3
(motif, respos) = sm.motif_from_position('P28482', '187')
assert motif == 'HTGFLTEYVATRWYR'
assert respos == 8
