Mercurial > repos > iuc > maf_stats

annotate maf_stats.py @ 1:a4dcd42b407b draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit 95632be9aca24769ac192f2eb6a0454af21eff9e"
author iuc
date Tue, 10 Nov 2020 18:47:45 +0000
parents 6d4d297433e9
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
1 #!/usr/bin/env python
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
2 # Dan Blankenberg
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
3 """
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
4 Reads a list of intervals and a maf. Outputs a new set of intervals with statistics appended.
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
5 """
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
6 from __future__ import print_function
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
7
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
8 import sys
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
9
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
10 import bx.intervals.io
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
11 from bx.bitset import BitSet
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
12 from galaxy.tools.util import maf_utilities
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
13
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
14
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
15 def __main__():
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
16 maf_source_type = sys.argv.pop(1)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
17 input_maf_filename = sys.argv[1].strip()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
18 input_interval_filename = sys.argv[2].strip()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
19 output_filename = sys.argv[3].strip()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
20 dbkey = sys.argv[4].strip()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
21 try:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
22 chr_col = int(sys.argv[5].strip()) - 1
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
23 start_col = int(sys.argv[6].strip()) - 1
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
24 end_col = int(sys.argv[7].strip()) - 1
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
25 except Exception:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
26 print("You appear to be missing metadata. You can specify your metadata by clicking on the pencil icon associated with your interval file.", file=sys.stderr)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
27 sys.exit()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
28 summary = sys.argv[8].strip()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
29 if summary.lower() == "true":
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
30 summary = True
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
31 else:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
32 summary = False
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
33
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
34 mafIndexFile = "%s/maf_index.loc" % sys.argv[9]
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
35 try:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
36 maf_index_filename = sys.argv[10].strip()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
37 except Exception:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
38 maf_index_filename = None
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
39 index = index_filename = None
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
40 if maf_source_type == "user":
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
41 # index maf for use here
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
42 index, index_filename = maf_utilities.open_or_build_maf_index(input_maf_filename, maf_index_filename, species=[dbkey])
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
43 if index is None:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
44 print("Your MAF file appears to be malformed.", file=sys.stderr)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
45 sys.exit()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
46 elif maf_source_type == "cached":
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
47 # access existing indexes
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
48 index = maf_utilities.maf_index_by_uid(input_maf_filename, mafIndexFile)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
49 if index is None:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
50 print("The MAF source specified (%s) appears to be invalid." % (input_maf_filename), file=sys.stderr)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
51 sys.exit()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
52 else:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
53 print('Invalid source type specified: %s' % maf_source_type, file=sys.stdout)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
54 sys.exit()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
55
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
56 out = open(output_filename, 'w')
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
57
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
58 num_region = None
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
59 num_bad_region = 0
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
60 species_summary = {}
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
61 total_length = 0
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
62 # loop through interval file
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
63 for num_region, region in enumerate(bx.intervals.io.NiceReaderWrapper(open(input_interval_filename, 'r'), chrom_col=chr_col, start_col=start_col, end_col=end_col, fix_strand=True, return_header=False, return_comments=False)):
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
64 src = "%s.%s" % (dbkey, region.chrom)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
65 region_length = region.end - region.start
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
66 if region_length < 1:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
67 num_bad_region += 1
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
68 continue
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
69 total_length += region_length
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
70 coverage = {dbkey: BitSet(region_length)}
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
71
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
72 for block in index.get_as_iterator(src, region.start, region.end):
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
73 for spec in maf_utilities.get_species_in_block(block):
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
74 if spec not in coverage:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
75 coverage[spec] = BitSet(region_length)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
76 for block in maf_utilities.iter_blocks_split_by_species(block):
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
77 if maf_utilities.component_overlaps_region(block.get_component_by_src(src), region):
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
78 # need to chop and orient the block
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
79 block = maf_utilities.orient_block_by_region(maf_utilities.chop_block_by_region(block, src, region), src, region, force_strand='+')
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
80 start_offset, alignment = maf_utilities.reduce_block_by_primary_genome(block, dbkey, region.chrom, region.start)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
81 for i in range(len(alignment[dbkey])):
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
82 for spec, text in alignment.items():
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
83 if text[i] != '-':
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
84 coverage[spec].set(start_offset + i)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
85 if summary:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
86 # record summary
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
87 for key in coverage.keys():
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
88 if key not in species_summary:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
89 species_summary[key] = 0
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
90 species_summary[key] = species_summary[key] + coverage[key].count_range()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
91 else:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
92 # print coverage for interval
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
93 coverage_sum = coverage[dbkey].count_range()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
94 out.write("%s\t%s\t%s\t%s\n" % ("\t".join(region.fields), dbkey, coverage_sum, region_length - coverage_sum))
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
95 keys = list(coverage.keys())
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
96 keys.remove(dbkey)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
97 keys.sort()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
98 for key in keys:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
99 coverage_sum = coverage[key].count_range()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
100 out.write("%s\t%s\t%s\t%s\n" % ("\t".join(region.fields), key, coverage_sum, region_length - coverage_sum))
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
101 if summary:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
102 out.write("#species\tnucleotides\tcoverage\n")
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
103 for spec in species_summary:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
104 out.write("%s\t%s\t%.4f\n" % (spec, species_summary[spec], float(species_summary[spec]) / total_length))
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
105 out.close()
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
106 if num_region is not None:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
107 print("%i regions were processed with a total length of %i." % (num_region + 1, total_length))
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
108 if num_bad_region:
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
109 print("%i regions were invalid." % (num_bad_region))
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
110 maf_utilities.remove_temp_index_file(index_filename)
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
111
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
112
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
113 if __name__ == "__main__":
6d4d297433e9 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/interval2maf/ commit f7b04f8c6edca90b1e9bc867a1569996aec76d7a"
iuc
parents:
diff changeset
114 __main__()