#!/usr/bin/python
import sys
from SOAPpy import WSDL

accession = sys.argv[1]

# create SOAP services entry points
GenomeAtlasWSDL = 'http://www.cbs.dtu.dk/ws/GenomeAtlas/GenomeAtlas_3_3.wsdl'

AtlasProxy = WSDL.Proxy(GenomeAtlasWSDL)

# input parameters for retrieving the sequence information
input = {'accession' : accession }

# method which will retrieve the sequences
sequences = AtlasProxy.getSeq(input)

print ">" + sequences.contig.sequencedata.sequence.id
for i in range(0,len(sequences.contig.sequencedata.sequence.seq),60):
	print sequences.contig.sequencedata.sequence.seq[i:i+60]