小生这厢有礼了(BioFaceBook Personal Blog)

#!/usr/bin/env python import os,re,sys,string,commands,getopt,subprocess,glob,csv import prettytable as pt from os import path #SL335752_kneaddata_paired_1_kneaddata_paired_1.fastq.gz_rep.mpa.txt_bracken.txt #SL311013_1_kneaddata_paired_1_kneaddata_paired_1.fastq.gz_rep_mpa.txt_bracken.txt def main(): dic = {} unique_speciesname = [] speciesname = [] samplenames = [] for d in os.listdir(‘.’): print(d) a = [] a = re.split(“_kneaddata_paired_1_”,d) if len(a) >=2: samplename = a[0] dic[a[0]] = {} #os.rename(d, newname) fh = open(d, ‘r’) fhlines […]