Mercurial > repos > devteam > cuffmerge
diff cuffmerge_wrapper.py @ 10:b6e3849293b1 draft
Uploaded
author | devteam |
---|---|
date | Fri, 19 Dec 2014 11:59:06 -0500 |
parents | 5b285b6e4ee3 |
children | 1707a530e598 |
line wrap: on
line diff
--- a/cuffmerge_wrapper.py Mon Jan 20 11:09:02 2014 -0500 +++ b/cuffmerge_wrapper.py Fri Dec 19 11:59:06 2014 -0500 @@ -1,7 +1,5 @@ #!/usr/bin/env python -# Supports Cuffmerge versions 1.3 and newer. - import optparse, os, shutil, subprocess, sys, tempfile def stop_err( msg ): @@ -15,13 +13,13 @@ parser.add_option( '-s', dest='use_seq_data', action="store_true", help='Causes cuffmerge to look into for fasta files with the underlying genomic sequences (one file per contig) against which your reads were aligned for some optional classification functions. For example, Cufflinks transcripts consisting mostly of lower-case bases are classified as repeats. Note that <seq_dir> must contain one fasta file per reference chromosome, and each file must be named after the chromosome, and have a .fa or .fasta extension.') parser.add_option( '-p', '--num-threads', dest='num_threads', help='Use this many threads to align reads. The default is 1.' ) - # Wrapper / Galaxy options. parser.add_option( '', '--index', dest='index', help='The path of the reference genome' ) parser.add_option( '', '--ref_file', dest='ref_file', help='The reference dataset from the history' ) # Outputs. parser.add_option( '', '--merged-transcripts', dest='merged_transcripts' ) + parser.add_option( '--min-isoform-fraction', dest='min_isoform_fraction' ) (options, args) = parser.parse_args() @@ -68,7 +66,8 @@ cmd += " -g %s " % options.ref_annotation if options.use_seq_data: cmd += " -s %s " % seq_path - + if options.min_isoform_fraction: + cmd += " --min-isoform-fraction %s " % (options.min_isoform_fraction) # Add input files to a file. inputs_file_name = tempfile.NamedTemporaryFile( dir="." ).name inputs_file = open( inputs_file_name, 'w' )