Skip to content

Instantly share code, notes, and snippets.

What would you like to do?
Script to get a list of all Schistosoma mansoni protein-coding genes from WormBase ParaSite
# script to retrieve a list of all protein-coding Schistosoma mansoni genes from wormbase parasite
# example script taken from
import requests, sys
server = ""
ext = "/rest-13/lookup/genome/schistosoma_mansoni_prjea36577?biotypes=protein_coding"
# took the PRJEA from
r = requests.get(server+ext, headers={ "Content-Type" : "application/json", "Accept" : ""})
if not r.ok:
decoded = r.json()
# go through the genes and print them out one at a time:
cnt = 0
for gene in decoded:
cnt += 1
print(cnt," ",gene["id"])
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
You can’t perform that action at this time.