# HG changeset patch # User Jim Johnson # Date 1381874543 18000 # Node ID 3099cec648e7b9f8cdd338ddeab3c5497ca0a781 # Parent 3e3ebdecb0e1fe62ea6fb92ef33f03c7e90a2d72 Update tool_dependencies, add help diff -r 3e3ebdecb0e1 -r 3099cec648e7 create_reference_dataset.xml --- a/create_reference_dataset.xml Fri Sep 06 08:20:41 2013 -0500 +++ b/create_reference_dataset.xml Tue Oct 15 17:02:23 2013 -0500 @@ -89,7 +89,12 @@ - + + + Examples: homo_sapiens, mus_musculus, rattus_norvegicus + ftp://ftp.ensembl.org/pub/release-$ensembl_version/fasta/$ensembl_organism/dna/$ensembl_prefix.$ensembl_genome_version.$ensembl_version.dna.chromosome.$chromosome.fa.gz + + diff -r 3e3ebdecb0e1 -r 3099cec648e7 defuse_results_to_vcf.py --- a/defuse_results_to_vcf.py Fri Sep 06 08:20:41 2013 -0500 +++ b/defuse_results_to_vcf.py Tue Oct 15 17:02:23 2013 -0500 @@ -97,6 +97,8 @@ ##INFO= ##INFO= ##INFO= +##INFO= +##INFO= ##INFO= ##INFO= ##INFO= @@ -255,8 +257,8 @@ info.append('EVERSION') if inversion: info.append('INVERSION') - info1 = [svtype,'MATEID=%s' % mate_id2] + info - info2 = [svtype,'MATEID=%s' % mate_id1] + info + info1 = [svtype,'MATEID=%s;MATELOC=%s:%d' % (mate_id2,gene_chromosome2,genomic_break_pos2)] + info + info2 = [svtype,'MATEID=%s;MATELOC=%s:%d' % (mate_id1,gene_chromosome1,genomic_break_pos1)] + info qual = int(float(fields[columns.index('probability')]) * 255) if columns.index('probability') else '.' vcf1 = '%s\t%d\t%s\t%s\t%s\t%s\t%s\t%s'% (gene_chromosome1,genomic_break_pos1, mate_id1, ref1, alt1, qual, filt, ';'.join(info1) ) vcf2 = '%s\t%d\t%s\t%s\t%s\t%s\t%s\t%s'% (gene_chromosome2,genomic_break_pos2, mate_id2, ref2, alt2, qual, filt, ';'.join(info2) ) @@ -265,7 +267,7 @@ write_vcf() except Exception, e: print >> sys.stderr, "failed: %s" % e - exit(1) + sys.exit(1) if __name__ == "__main__" : __main__() diff -r 3e3ebdecb0e1 -r 3099cec648e7 tool_dependencies.xml --- a/tool_dependencies.xml Fri Sep 06 08:20:41 2013 -0500 +++ b/tool_dependencies.xml Tue Oct 15 17:02:23 2013 -0500 @@ -1,10 +1,10 @@ - + - +