ó ÿfmSc@sˆddlZddlmZddlmZejdƒrMejdƒZndZdeZd„Zd„Z d „Z ed „Z dS( iÿÿÿÿN(tGenome(tenvirontCISTEMATIC_ROOTs /proj/genomes"%s/S_purpuratus/spurpuratus.genedbcCsÙg}d}tdd|ƒ}t|dƒ}|jƒ}x|dkrÊg}d}|jƒ}d|kry|d} n|djƒ} | jd ƒ} | d} |jƒ} xU| dkr| dd kr| jƒ} |t| ƒ7}|j| ƒ|jƒ} q±Wtj|dƒ}|d kred | GH|j d| f|d t |ƒƒ|j | | dƒn\d|| f}tdt |fdƒ}|j |ƒ|jƒd|GH|j | |dƒ| }q<W|jƒdS(Nit spurpuratustdbFiletrtt purpuratusiiÿÿÿÿt_t>iPÃsAdded contig %s to databaset chromosometdbs%s%s.bins%s%stws Added contig file %s to databasetfile(Rtopentreadlinetsplittstriptlentappendtstringtjoint addSequencetstrtaddChromosomeEntrytcisRoottwritetclose(R t chromPatht chromOutPathtseqArraytseqLentspGenometinFiletheadertfieldstidparttpartstchromIDt currentLinetlineSeqtseqt outFileNametoutFile((sh/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/spurpuratus.pytloadChromosome+sB          "    cCs)tdddd|ƒ}|j|ƒdS(NRtversions2.1R(Rt createGeneDB(R R ((sh/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/spurpuratus.pyt createDBFileUscCs&tdddd|ƒ}|jƒdS(NRR-s2.1R(Rt createIndices(R R ((sh/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/spurpuratus.pytcreateDBindicesZscCsTdt}d}d|GHt|ƒdGHt|||ƒdGHt|ƒd|GHdS(Ns%s/download/Spur2.1_Nmasked.txts/S_purpuratus/sCreating database %ssLoading genomic sequencesCreating IndicessFinished creating database %s(RR/R,R1(R t chromoPatht chromoOutPath((sh/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/spurpuratus.pytbuildSpurpuratusDB_s    ( Rtcistematic.genomesRtosRtgetRtgeneDBR,R/R1R4(((sh/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM-2014-05-09-fix/cistematic/genomes/spurpuratus.pyts   *