annotate fimo_wrapper.py @ 0:fd522a964017 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
author iuc
date Tue, 22 Dec 2015 17:01:51 -0500
parents
children eca84de658b0
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
1 #!/usr/bin/env python
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
2 import argparse
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
3 import os
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
4 import shutil
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
5 import string
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
6 import subprocess
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
7 import sys
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
8 import tempfile
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
9
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
10 BUFFSIZE = 1048576
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
11 # Translation table for reverse Complement, with ambiguity codes.
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
12 DNA_COMPLEMENT = string.maketrans("ACGTRYKMBDHVacgtrykmbdhv", "TGCAYRMKVHDBtgcayrmkvhdb")
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
13
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
14
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
15 def reverse(sequence):
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
16 # Reverse sequence string.
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
17 return sequence[::-1]
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
18
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
19
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
20 def dna_complement(sequence):
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
21 # Complement DNA sequence string.
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
22 return sequence.translate(DNA_COMPLEMENT)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
23
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
24
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
25 def dna_reverse_complement(sequence):
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
26 # Returns the reverse complement of the sequence.
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
27 sequence = reverse(sequence)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
28 return dna_complement(sequence)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
29
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
30
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
31 def stop_err(msg):
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
32 sys.stderr.write(msg)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
33 sys.exit(1)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
34
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
35 parser = argparse.ArgumentParser()
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
36 parser.add_argument('--input_motifs', dest='input_motifs', help='MEME output formatted files for input to fimo')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
37 parser.add_argument('--input_fasta', dest='input_fasta', help='Fassta sequence file')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
38 parser.add_argument('--options_type', dest='options_type', help='Basic or Advance options')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
39 parser.add_argument('--input_psp', dest='input_psp', default=None, help='File containing position specific priors')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
40 parser.add_argument('--input_prior_dist', dest='input_prior_dist', default=None, help='File containing binned distribution of priors')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
41 parser.add_argument('--alpha', dest='alpha', type=float, default=1.0, help='The alpha parameter for calculating position specific priors')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
42 parser.add_argument('--bgfile', dest='bgfile', default=None, help='Background file type, used only if not "default"')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
43 parser.add_argument('--max_strand', action='store_true', help='If matches on both strands at a given position satisfy the output threshold, only report the match for the strand with the higher score')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
44 parser.add_argument('--max_stored_scores', dest='max_stored_scores', type=int, help='Maximum score count to store')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
45 parser.add_argument('--motif', dest='motifs', action='append', default=[], help='Specify motif by id')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
46 parser.add_argument('--motif_pseudo', dest='motif_pseudo', type=float, default=0.1, help='Pseudocount to add to counts in motif matrix')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
47 parser.add_argument('--no_qvalue', action='store_true', help='Do not compute a q-value for each p-value')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
48 parser.add_argument('--norc', action='store_true', help='Do not score the reverse complement DNA strand')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
49 parser.add_argument('--output_path', dest='output_path', help='Output files directory')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
50 parser.add_argument('--parse_genomic_coord', action='store_true', help='Check each sequence header for UCSC style genomic coordinates')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
51 parser.add_argument('--qv_thresh', action='store_true', help='Use q-values for the output threshold')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
52 parser.add_argument('--thresh', dest='thresh', type=float, help='p-value threshold')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
53 parser.add_argument('--gff_output', dest='gff_output', help='Gff output file')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
54 parser.add_argument('--html_output', dest='html_output', help='HTML output file')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
55 parser.add_argument('--interval_output', dest='interval_output', help='Interval output file')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
56 parser.add_argument('--txt_output', dest='txt_output', help='Text output file')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
57 parser.add_argument('--xml_output', dest='xml_output', help='XML output file')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
58 args = parser.parse_args()
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
59
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
60 fimo_cmd_list = ['fimo']
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
61 if args.options_type == 'advanced':
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
62 fimo_cmd_list.append('--alpha %4f' % args.alpha)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
63 if args.bgfile is not None:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
64 fimo_cmd_list.append('--bgfile "%s"' % args.bgfile)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
65 if args.max_strand:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
66 fimo_cmd_list.append('--max-strand')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
67 fimo_cmd_list.append('--max-stored-scores %d' % args.max_stored_scores)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
68 if len(args.motifs) > 0:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
69 for motif in args.motifs:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
70 fimo_cmd_list.append('--motif "%s"' % motif)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
71 fimo_cmd_list.append('--motif-pseudo %4f' % args.motif_pseudo)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
72 if args.no_qvalue:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
73 fimo_cmd_list.append('--no-qvalue')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
74 if args.norc:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
75 fimo_cmd_list.append('--norc')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
76 if args.parse_genomic_coord:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
77 fimo_cmd_list.append('--parse-genomic-coord')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
78 if args.qv_thresh:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
79 fimo_cmd_list.append('--qv-thresh')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
80 fimo_cmd_list.append('--thresh %4f' % args.thresh)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
81 if args.input_psp is not None:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
82 fimo_cmd_list.append('--psp "%s"' % args.input_psp)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
83 if args.input_prior_dist is not None:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
84 fimo_cmd_list.append('--prior-dist "%s"' % args.input_prior_dist)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
85 fimo_cmd_list.append('--o "%s"' % (args.output_path))
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
86 fimo_cmd_list.append('--verbosity 1')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
87 fimo_cmd_list.append(args.input_motifs)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
88 fimo_cmd_list.append(args.input_fasta)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
89
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
90 fimo_cmd = ' '.join(fimo_cmd_list)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
91
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
92 try:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
93 tmp_stderr = tempfile.NamedTemporaryFile()
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
94 proc = subprocess.Popen(args=fimo_cmd, shell=True, stderr=tmp_stderr)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
95 returncode = proc.wait()
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
96 tmp_stderr.seek(0)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
97 stderr = ''
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
98 try:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
99 while True:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
100 stderr += tmp_stderr.read(BUFFSIZE)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
101 if not stderr or len(stderr) % BUFFSIZE != 0:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
102 break
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
103 except OverflowError:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
104 pass
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
105 if returncode != 0:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
106 stop_err(stderr)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
107 except Exception, e:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
108 stop_err('Error running FIMO:\n%s' % str(e))
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
109
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
110 shutil.move(os.path.join(args.output_path, 'fimo.txt'), args.txt_output)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
111 shutil.move(os.path.join(args.output_path, 'fimo.gff'), args.gff_output)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
112 shutil.move(os.path.join(args.output_path, 'fimo.xml'), args.xml_output)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
113 shutil.move(os.path.join(args.output_path, 'fimo.html'), args.html_output)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
114
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
115 out_file = open(args.interval_output, 'wb')
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
116 out_file.write("#%s\n" % "\t".join(("chr", "start", "end", "pattern name", "score", "strand", "matched sequence", "p-value", "q-value")))
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
117 for line in open(args.txt_output):
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
118 if line.startswith('#'):
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
119 continue
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
120 fields = line.rstrip("\n\r").split("\t")
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
121 start, end = int(fields[2]), int(fields[3])
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
122 sequence = fields[7]
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
123 if start > end:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
124 # Flip start and end and set strand.
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
125 start, end = end, start
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
126 strand = "-"
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
127 # We want sequences relative to strand; FIMO always provides + stranded sequence.
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
128 sequence = dna_reverse_complement(sequence)
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
129 else:
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
130 strand = "+"
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
131 # Make 0-based start position.
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
132 start -= 1
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
133 out_file.write("%s\n" % "\t".join([fields[1], str(start), str(end), fields[0], fields[4], strand, sequence, fields[5], fields[6]]))
fd522a964017 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit e96df94dba60050fa28aaf55b5bb095717a5f260
iuc
parents:
diff changeset
134 out_file.close()