def analyse(uniprot):
print('***************** ANALYSIS *******************************')
p = ProteinAnalyser(taxid='9606', uniprot=uniprot).load()
p.mutation = f'{p.sequence[65]}66W'
p.predict_effect()
p.analyse_structure()
print(p)
import json
json.dump(p.asdict(), open('test.json', 'w'))
def test_ProteinAnalyser():
p = ProteinAnalyser(uniprot='Q86V25').load()
print(p)
p.mutation = Mutation('p.N127W')
p.analyse_structure()
print(p.get_features_near_position())
print(p.get_gnomAD_near_position())
print(p.model.get_structure_neighbours())
print(p.get_superficiality())