comparison imgtconvert.py @ 6:5b030e48b308 draft

Uploaded
author davidvanzessen
date Mon, 07 Jul 2014 09:54:42 -0400
parents 387fce4a1dd4
children
comparison
equal deleted inserted replaced
5:387fce4a1dd4 6:5b030e48b308
45 for i in range(0, step): 45 for i in range(0, step):
46 triplets.append((files[i], files[i + step], files[i + step + step])) 46 triplets.append((files[i], files[i + step], files[i + step + step]))
47 47
48 outFile = args.output 48 outFile = args.output
49 49
50 fSummary = pd.read_csv(triplets[0][0], sep="\t", low_memory=False) 50 #fSummary = pd.read_csv(triplets[0][0], sep="\t", low_memory=False)
51 fSequence = pd.read_csv(triplets[0][1], sep="\t", low_memory=False) 51 fSummary = pd.read_csv(triplets[0][0], sep="\t")
52 fJunction = pd.read_csv(triplets[0][2], sep="\t", low_memory=False) 52 #fSequence = pd.read_csv(triplets[0][1], sep="\t", low_memory=False)
53 fSequence = pd.read_csv(triplets[0][1], sep="\t")
54 #fJunction = pd.read_csv(triplets[0][2], sep="\t", low_memory=False)
55 fJunction = pd.read_csv(triplets[0][2], sep="\t")
53 tmp = fSummary[["Sequence ID", "JUNCTION frame", "V-GENE and allele", "D-GENE and allele", "J-GENE and allele"]] 56 tmp = fSummary[["Sequence ID", "JUNCTION frame", "V-GENE and allele", "D-GENE and allele", "J-GENE and allele"]]
54 57
55 tmp["CDR1 Seq"] = fSequence["CDR1-IMGT"] 58 tmp["CDR1 Seq"] = fSequence["CDR1-IMGT"]
56 tmp["CDR1 Length"] = fSummary["CDR1-IMGT length"] 59 tmp["CDR1 Length"] = fSummary["CDR1-IMGT length"]
57 60