Esempio n. 1
0
def lanceModifHeader():
	
	liStrains = ['NCYC543','62-1041','62-196','CBS2861','CBS4104','CBS6545','CBS10367','CBS10368','DBVPG3108','dd281a','55-86_1','67-588','77-1003','CBS4568','CBS5828','CBS6547','CBS6626','DBVPG4002',"CBS3082b","CBS3082a","NRBC1892","NRBC1811","NRBC10955","NRBC10572","NRBC101999","DBVPG3452","CBS6546","CBS10369","68917-2"]
	liStrains = ['55-86_1','67-588','77-1003','CBS4568','CBS5828','CBS6547','CBS6626','DBVPG4002']
		
	liKmer = [75,63,83,57]
        for strain in liStrains:
                print "traite strain %s" % strain
                for kmer in liKmer:
			ace = "ace%s_%smerPrep.scafSeq" % (strain,kmer)
			fileOut =  "ace%s_%smer.fasta" % (strain,kmer)
			suff = "%s_%s" % (strain,kmer)
			if os.path.isfile(ace):
				fasta.ajoutSuffHeader(ace,suff,fileOut)
Esempio n. 2
0
def lanceALaChaineModifHeader():

	liStrains = ['NCYC543','62-1041','62-196','CBS2861','CBS4104','CBS6545','CBS10367','CBS10368','DBVPG3108','dd281a','55-86_1','67-588','77-1003','CBS4568','CBS5828','CBS6547','CBS6626','DBVPG4002',"CBS3082b","CBS3082a","NRBC1892","NRBC1811","NRBC10955","NRBC10572","NRBC101999","DBVPG3452","CBS6546","CBS10369","68917-2"]
	
	for strain in liStrains:
		print "traite strain %s" % strain
		fileIn1 = "ace%s_63merPrep.scafSeq" % strain
		fileOut1 =  "ace%s_63compl.fasta" % strain
		suff1 = "%s_63compl" % strain

		fileIn2 = "ace%s_75merPrep.scafSeq" % strain
		fileOut2 =  "ace%s_75compl.fasta" % strain
		suff2 = "%s_75compl" % strain

		if os.path.isfile(fileIn1):
			fasta.ajoutSuffHeader(fileIn1,suff1,fileOut1)
		if os.path.isfile(fileIn2):
			fasta.ajoutSuffHeader(fileIn2,suff2,fileOut2)