def getPathsReacts(project): kegg = Kegg(project) # Get the pathway - reaction links paths = {} for pR in kegg.getPathReacts(): if pR.path_id in ['path:rn01100','path:rn01110','path:rn01120']: continue if pR.path_id not in paths: paths[pR.path_id] = [] paths[pR.path_id].append(pR.re_id) return paths