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)
