def dbCpG(inFileName): """ deprecated """ vec_len=[] file_db=open(inFileName) seq=Bioseq() numseq=0 while 1: seq.read(file_db) if seq.sequence==None: break numseq=numseq+1 lr,cg,c,g=seq.countCpG() print 'sequence #',numseq,'=',seq.header[0:40],'...]',lr,cg,c,g file_db.close()
