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)
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)