annotate capture_kits.py @ 0:1a23ea467feb default tip

intial commit
author Yusuf Ali <ali@yusuf.email>
date Thu, 26 Mar 2015 09:36:17 -0600
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
1 import os
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
2 import re
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
3 import sys
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
4 import operator
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
5 import csv
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
6 from galaxy import config
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
7
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
8 # get tool-data path
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
9 configur = config.Configuration()
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
10 kitDir = configur.resolve_path("tool-data")
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
11
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
12 # determine if config file exists
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
13 if not os.path.exists( kitDir + "/miseq_bam_variants.loc" ):
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
14 kitDir = "/export/achri_galaxy/dbs/CaptureKits/";
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
15 else:
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
16 with open(kitDir + "/miseq_bam_variants.loc", "r") as tsv:
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
17 for line in csv.reader(tsv, delimiter="\t"):
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
18 if line[0] == 'capture_kits_dir':
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
19 kitDir = line[1]
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
20
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
21 def kit_fileOptions():
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
22 list = os.listdir(kitDir);
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
23 list.sort()
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
24 pattern = re.compile('(.*)$')
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
25 fileOptions = [(s) for s in list if os.path.exists(kitDir + s)]
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
26 ds = [pattern.match(s) for s in fileOptions]
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
27 datasets = [(m.group(1), m.group(1), False) for m in ds if m]
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
28 return datasets
1a23ea467feb intial commit
Yusuf Ali <ali@yusuf.email>
parents:
diff changeset
29