def statsScaff(): #liStrains = ['55-86_1','62-196','62-1041','67-588','77-1003','CBS2861','CBS4104','CBS4568','CBS5828','CBS6545','CBS6547','CBS6626','CBS10367','CBS10368','dd281a','DBVPG3108','DBVPG4002','NCYC543'] liStrains = ['55-86_1','62-196','62-1041','67-588','77-1003','CBS2861','CBS4104','CBS4568','CBS5828','CBS6545','CBS6547','CBS6626','CBS10367','CBS10368','DBVPG3108','DBVPG4002','NCYC543'] print "strain\tnbSacff\ttotal size\tmax scaff\tN50\tmean size" for strain in liStrains: fic = "ace%s_75part.fasta" % strain #print fic fasta.longSeqMultiFasta(fic,strain)
def statsContigs(fic): fasta.longSeqMultiFasta(fic)