ó ÷VêNc@sšddlZddlmZddlmZejdƒrMejdƒZndZdeZd„Zd„Z d „Z d „Z d „Z ed „Z dS( iÿÿÿÿN(tGenome(tenvirontCISTEMATIC_ROOTs /proj/genomes%s/C_brenneri/cbrenneri.genedbcCs«g}d}tdd|ƒ}t|dƒ}|jƒ}xa|dkrœg}d}|jƒdjƒ}|d} |jƒ} xU| dkr×| ddkr×| jƒ} |t| ƒ7}|j| ƒ|jƒ} qƒWtj|dƒ} |dkr7d | GH|j d| f| d t |ƒƒ|j | | d ƒn\d || f} td t | fdƒ}|j | ƒ|jƒd| GH|j | | dƒ| }q<W|jƒdS(Nit cbrenneritdbFiletrtit>i †sAdded contig %s to databaset chromosometdbs%s%s.bins%s%stws Added contig file %s to databasetfile(Rtopentreadlinetstriptsplittlentappendtstringtjoint addSequencetstrtaddChromosomeEntrytcisRoottwritetclose(R t chromPatht chromOutPathtseqArraytseqLentcbGenometinFiletheadert chromHeadertchromIDt currentLinetlineSeqtseqt outFileNametoutFile((sW/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM/cistematic/genomes/cbrenneri.pytloadChromosome+s:       "    c Cs$tdd|ƒ}t|dƒ}i}i}i}i}g}x^|D]V} | ddkrbqFn| ddkrxqFn| d jdƒ} | d d kr±| d d kr±qFn| d jd ƒ} | d} d| f} | d}| djƒ|| <|dkrd|| t gidVersiontstarttstopR?tchrom((sW/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM/cistematic/genomes/cbrenneri.pytloadFeatureEntriesƒs2      &cCs)tdddd|ƒ}|j|ƒdS(NRtversiont PB2801_001R(Rt createGeneDB(R R((sW/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM/cistematic/genomes/cbrenneri.pyt createDBFile§scCs&tdddd|ƒ}|jƒdS(NRRJRKR(Rt createIndices(R R((sW/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM/cistematic/genomes/cbrenneri.pytcreateDBindices¬scCs‚dt}dt}d}d|GHt|ƒdGHt||ƒdGHt||ƒdGHt|||ƒdGHt|ƒd |GHdS( Ns %s/download/PB2801_2007feb09.gffs"%s/download/PB2801_supercontigs.fas /C_brenneri/sCreating database %ssAdding gene entriessAdding feature entriessLoading genomic sequencesCreating IndicessFinished creating database %s(RRMR@RIR(RO(R tgffPatht chromoPatht chromoOutPath((sW/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM/cistematic/genomes/cbrenneri.pytbuildCbrenneriDB±s       (Rtcistematic.genomesRtosRtgetRtgeneDBR(R@RIRMRORS(((sW/woldlab/castor/data00/home/georgi/code/erange-4.0a-BAM/cistematic/genomes/cbrenneri.pyts   % 3 $