ó ˙fmSc@sŻddlZddlmZddlmZejdƒrMejdƒZndZdeZd„Zd„Z d „Z d „Z d „Z d „Z d „Zedd„ZdS(i˙˙˙˙N(tGenome(tenvirontCISTEMATIC_ROOTs /proj/genomes%s/C_elegans/celegans.genedbc CsŮg}tdd|ƒ}t|dƒ}|jƒ}x!|D]}|j|jƒƒq:Wtj|dƒ}t|ƒ} | dkr‰dGHnd|GHtdt|fd ƒ} | j |ƒ| j ƒ|j ||d ƒdS( NtceleganstdbFiletrtis#Problems reading sequence from fileswriting to file %ss%s%stwtfile( RtopentreadlinetappendtstriptstringtjointlentcisRoottwritetclosetaddChromosomeEntry( tdbtchromIDt chromPathtchromOuttseqArraytceGenometinFiletlinetseqtseqLentoutFile((se/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/celegans.pytloadChromosome+s       c Cstdd|ƒ}t|dƒ}g}xÓ|D]Ë}|ddkrJq.n|jdƒ}|ddkr|dd krq.n|d d krĽ|d d krĽq.n|d jdƒ}|djdƒ}|ddtjkr|dd } |dd} n|d} d} d|d| f} d| f} d} | dkrzyt| jƒƒd} Wqzd|dGHq.qzXnt|dƒd}t|dƒd}|d}|djƒ}|dkrŃd}nd}|j | ||||d | fƒq.Wdt |ƒGH|j |ƒdS(NRRRit#s itCoding_transcripttmiRNAit TranscripttmiRNA_primary_transcriptit"t.i˙˙˙˙tas%s.%si`s problem processing %s - skippingiiit+tFtRsAdding %d gene entries( RR tsplitR tletterstordtlowertintR R RtaddGeneEntryBatch(RtgffFileRtgeneFilet geneEntriesRtfieldtgidrevtgiddotstgidGenet gidLettertgidtgeneIDt gidVersiontstarttstoptsensetchrom((se/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/celegans.pytloadGeneEntries?sJ           &c Csśtdd|ƒ}t|dƒ}g}i}idd6dd6dd 6}xK|D]C}|d d krkqOn|jd ƒ}|d d!krqOn|d dkrź|djƒ|krźqOn|d d"krŐd} n||djƒ} |djdƒ} | d jdƒ} | d dtjkrE| d d } | d d} n| d } d} d| d | f}d|f}d }| dkržyt| jƒƒd}Wqžd| d GHqOqžXnt|dƒd }t|dƒd }|d}|d jƒ}|dkrd}nd}||kr4g||R?((se/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/celegans.pytloadFeatureEntriesps`   &           "$c CsĚg}t|dƒ}|jƒ}|jƒtdd|ƒ}xl|D]d}|jdƒ}yE|djƒ}d|d|df} |jd|f| fƒWqDqDXqDWd t|ƒGH|j|ƒdS( NRRRt,is%s %siisAdding %d annotations( R t readlinesRRR+R R RtaddAnnotationBatch( Rt geneIDPathtgeneAnnotationst geneIDFiletlinesRRR4R9R:((se/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/celegans.pytloadGeneAnnotations˛s   c Csötdddd|ƒ}t|dƒ}t|dƒ}t|dƒ}|jƒ}|jƒi} i} xT|D]L} | jdƒ} t| djƒƒdkrn| djƒ| | dn| jd ƒ}|d}|| krŒ| |}n|d d krŠ|d }n|djdƒ}|d }|d}|djdƒ}|}|d}t|ƒdkr"d|jƒ|df}ny!tj||d ddƒ}Wnd|GHd}nX|d}|| krg| || |j |ƒ|j d|f|||||||dfƒq>q>Wdt|ƒGH|j |ƒdS(NRtversiontWS200RRRNiiit!s i˙˙˙˙R'it ii t|is%s|%st'tps!could no map %s - using GOID onlyRi sAdding %d GO entries( RR RORR+RR R treplaceR taddGoInfoBatch(RtgoPatht goDefPathRQRRSt goDefFiletgoFileRTt geneIDmaptseenGORR4t goDefEntriestgoDefst goDefEntrytcolst goEntriestgoArraytfieldstnamet GOIDarraytGOIDtobjTypetobjNametgIDtisNottGOtermtevidence((se/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/celegans.pytloadGeneOntologyĹsd     )         !     6cCs)tdddd|ƒ}|j|ƒdS(NRRVRWR(Rt createGeneDB(RR((se/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/celegans.pyt createDBFilescCs&tdddd|ƒ}|jƒdS(NRRVRWR(Rt createIndices(RR((se/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/celegans.pytcreateDBindices sRcCs;|dkrdt}nd|}d|}d|}id|d6d|d 6d |d 6d |d 6d|d6d|d6}d|}d|GHt|ƒdGHt||ƒdGHt||ƒdGHt||ƒdGHt||||ƒx2|D]*}d|GHt||||d|ƒqńWdGHt|ƒd|GHdS(NRs %s/download/s%sgeneIDs.WS200s%sGO.terms_and_idss%sgene_association.wbs%sCHROMOSOME_I_softmasked.dnatIs%sCHROMOSOME_II_softmasked.dnatIIs%sCHROMOSOME_III_softmasked.dnatIIIs%sCHROMOSOME_IV_softmasked.dnatIVs%sCHROMOSOME_V_softmasked.dnatVs%sCHROMOSOME_X_softmasked.dnatXs%selegansWS200.gffsCreating database %ssAdding gene entriessAdding feature entriessAdding gene annotationssAdding gene ontologysLoading chromosome %ss/C_elegans/chr%s.binsCreating IndicessFinished creating database %s(RRwR@RMRURuRRy(Rt downloadRootRQR`R_tchromostgffPathR((se/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/celegans.pytbuildCelegansDBs8                  (R tcistematic.genomesRtosRtgetRtgeneDBRR@RMRURuRwRyRƒ(((se/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/celegans.pyts    1 B  @