ó ÿfmSc@s£ddlZddlmZddlmZejdƒrMejdƒZndZdeZd„Zd„Z d „Z d „Z d „Z d „Z ed „ZdS(iÿÿÿÿN(tGenome(tenvirontCISTEMATIC_ROOTs /proj/genomes %s/M_domestica/mdomestica.genedbcCs›g}d}tdd|ƒ}t|dƒ}|jƒ}xQ|dkrŒg}d}|jƒd}|jƒ} xU| dkrÇ| ddkrÇ| jƒ} |t| ƒ7}|j| ƒ|jƒ} qsWtj|dƒ} |dkr'd |GH|jd|f| d t |ƒƒ|j ||d ƒn\d ||f} td t | fdƒ} | j | ƒ| j ƒd| GH|j || dƒ| }q<W|j ƒdS(Nit mdomesticatdbFiletrtit>i ¡sAdded contig %s to databaset chromosometdbs%s%s.bins%s%stws Added contig file %s to databasetfile(Rtopentreadlinetstriptlentappendtstringtjoint addSequencetstrtaddChromosomeEntrytcisRoottwritetclose(R t chromPatht chromOutPathtseqArraytseqLentmdGenometinFiletheadertchromIDt currentLinetlineSeqtseqt outFileNametoutFile((sg/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/mdomestica.pytloadChromosome+s8      "    c Csôg}tdd|ƒ}t|dƒ}xª|D]¢}|jdƒ}|d}t|dƒ}t|dƒ} |d} |d } | d kr–d } nd } d|f} d } |j| | || | d | fƒq.Wdt|ƒGH|j|ƒdS(s< FIXME - NEED TO DEAL WITH ALTERNATIVE SPLICING ENTRIES RRRs iiiiit+tFtRtgenesAdding %d gene entriesN(RR tsplittintRRtaddGeneEntryBatch(R tgFilet geneEntriesRtgeneFiletlinetcolstgidtstarttstoptsensetchromtgeneIDt gidVersion((sg/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/mdomestica.pytloadGeneEntriesNs$       &c CsÌg}t|dƒ}tdd|ƒ}x‚|D]z}yj|jdƒ}|d}|d}t|ƒdkr|jd|ftj|jƒddƒfƒnWq.q.Xq.Wd t|ƒGH|j|ƒdS( NRRRs iit'tpsAdding %d annotations( R RR+RRRtreplaceRtaddAnnotationBatch( R t annotPathtgeneAnnotationst annotFileRR1R2tlocIDtgeneDesc((sg/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/mdomestica.pytloadGeneAnnotationshs   5c Cs@t|dƒ}idd6dd6dd6}g}g}xË|D]Ã}|jdƒ}t|dƒ}|d jd ƒ} |d jd ƒ} |d } ||d } t|dƒd } t|dƒd }|d}yd|f}Wn q=nXd}||kr d}n |j|ƒxæt|ƒD]Ø}t| |ƒd }t| |ƒd }||kr|| kr|j||| ||| dfƒq$||krÛ| dkr°d}nd}|j|d | ||| |fƒq$|| kr'| dkrüd}nd}|j|d | ||| |fƒq$|| kr™| dkrHd}nd}|j|d | || | dfƒ|j|d | | d || |fƒq$| dkr®d}nd}|j|d | ||| dfƒ|j|d | ||d | |fƒq$Wq=W|jƒtdd|ƒ}dt|ƒGH|j|ƒdS(NRR(R'R)t-t.s iit,i iiiiiRt1t2tCDSt5UTRt3UTRRsAdding %d features( R R+R,RtrangeRRRtaddFeatureEntryBatch(R tgfileR0t senseArrayt seenArrayt insertArraytgeneLinet geneFieldstexonNumt exonStartst exonStopsR7R6tgstoptgstarttgeneidR8R9tindextestarttestoptfTypeR((sg/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/mdomestica.pytloadGeneFeatureszsl        %   %   %   ")  ". cCs#tdd|ƒ}|j|ƒdS(NRR(Rt createGeneDB(R R((sg/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/mdomestica.pyt createDBFileÂscCs tdd|ƒ}|jƒdS(NRR(Rt createIndices(R R((sg/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/mdomestica.pytcreateDBindicesÇscCs‚dt}dt}d}d|GHt|ƒdGHt||ƒdGHt||ƒdGHt|||ƒdGHt|ƒd |GHdS( Ns%s/download/mondom/genscan.txts%s/download/mondom/softMask.fas /M_domestica/sCreating database %ssAdding gene entriessAdding gene featuressLoading chromosomessCreating IndicessFinished creating database %s(RRaR:R_R&Rc(R tgenePatht chromoPatht chromoOutPath((sg/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/mdomestica.pytbuildMdomesticaDBÌs       (Rtcistematic.genomesRtosRtgetRtgeneDBR&R:RDR_RaRcRg(((sg/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/mdomestica.pyts   #   H