Skip to content
Snippets Groups Projects
Commit 80db7b5b authored by Vermaat's avatar Vermaat
Browse files

Add chromAccession.py to SOAP tools

git-svn-id: https://humgenprojects.lumc.nl/svn/mutalyzer/trunk@574 eb6bd6ab-9ccd-42b9-aceb-e2899b4a52f1
parent 8300b90f
No related branches found
No related tags found
No related merge requests found
#!/usr/bin/env python
"""
Get chromosome accession number from a Mutalyzer installation.
Usage:
{command} chromosome [build]
chromosome: Chromosome to get accession number for, e.g. 'chr2'.
build: Human genome build, 'hg18' or 'hg19' (default: 'hg19').
The accession number is retrieved from the Mutalyzer SOAP webservice and
printed to standard output.
"""
from mutalyzer.util import monkey_patch_suds; monkey_patch_suds()
import sys
from suds.client import Client
from mutalyzer.util import format_usage
WSDL_LOCATION = 'http://localhost/mutalyzer/services/?wsdl'
def main(chromosome, build='hg19'):
"""
Get accession number and print it to standard output.
"""
service = Client(WSDL_LOCATION, cache=None).service
result = service.chromAccession(build=build, name=chromosome)
if result:
print 'Chromosome accession number: %s' % str(result)
if __name__ == '__main__':
if len(sys.argv) < 2:
print format_usage()
sys.exit(1)
main(*sys.argv[1:])
...@@ -22,8 +22,7 @@ from suds.client import Client ...@@ -22,8 +22,7 @@ from suds.client import Client
from mutalyzer.util import format_usage from mutalyzer.util import format_usage
#WSDL_LOCATION = 'http://localhost/mutalyzer/services/?wsdl' WSDL_LOCATION = 'http://localhost/mutalyzer/services/?wsdl'
WSDL_LOCATION = 'https://mutalyzer.nl/services/?wsdl'
def main(transcript, variant, build='hg19'): def main(transcript, variant, build='hg19'):
...@@ -34,7 +33,6 @@ def main(transcript, variant, build='hg19'): ...@@ -34,7 +33,6 @@ def main(transcript, variant, build='hg19'):
result = service.mappingInfo(LOVD_ver='3.0-beta-06', build=build, result = service.mappingInfo(LOVD_ver='3.0-beta-06', build=build,
accNo=transcript, variant=variant) accNo=transcript, variant=variant)
print result
if result: if result:
print 'Genomic start: %s' % result.start_g print 'Genomic start: %s' % result.start_g
print 'Genomic end: %s' % result.end_g print 'Genomic end: %s' % result.end_g
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment