annotate sashimi-plot.py @ 0:9304dd9a16a2 draft

"planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
author artbio
date Fri, 23 Aug 2019 11:38:29 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1 #!/usr/bin/env python
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
2
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
3 # Import modules
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
4 import copy
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
5 import os
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
6 import re
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
7 import subprocess as sp
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
8 import sys
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
9 from argparse import ArgumentParser
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
10 from collections import OrderedDict
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
11
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
12
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
13 def define_options():
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
14 # Argument parsing
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
15 parser = ArgumentParser(description="""Create sashimi plot for a given
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
16 genomic region""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
17 parser.add_argument("-b", "--bam", type=str,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
18 help="""
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
19 Individual bam file or file with a list of bam files.
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
20 In the case of a list of files the format is tsv:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
21 1col: id for bam file,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
22 2col: path of bam file,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
23 3+col: additional columns
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
24 """)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
25 parser.add_argument("-c", "--coordinates", type=str,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
26 help="Genomic region. Format: chr:start-end (1-based)")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
27 parser.add_argument("-o", "--out-prefix", type=str, dest="out_prefix",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
28 default="sashimi",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
29 help="Prefix for plot file name [default=%(default)s]")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
30 parser.add_argument("-S", "--out-strand", type=str, dest="out_strand",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
31 default="both", help="""Only for --strand other than
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
32 'NONE'. Choose which signal strand to plot:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
33 <both> <plus> <minus> [default=%(default)s]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
34 parser.add_argument("-M", "--min-coverage", type=int, default=1,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
35 dest="min_coverage", help="""Minimum number of reads
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
36 supporting a junction to be drawn [default=1]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
37 parser.add_argument("-j", "--junctions-bed", type=str, default="",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
38 dest="junctions_bed", help="""Junction BED file name
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
39 [default=no junction file]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
40 parser.add_argument("-g", "--gtf",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
41 help="Gtf file with annotation (only exons is enough)")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
42 parser.add_argument("-s", "--strand", default="NONE", type=str,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
43 help="""Strand specificity: <NONE> <SENSE> <ANTISENSE>
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
44 <MATE1_SENSE> <MATE2_SENSE> [default=%(default)s]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
45 parser.add_argument("--shrink", action="store_true",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
46 help="""Shrink the junctions by a factor for nicer
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
47 display [default=%(default)s]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
48 parser.add_argument("-O", "--overlay", type=int,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
49 help="Index of column with overlay levels (1-based)")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
50 parser.add_argument("-A", "--aggr", type=str, default="",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
51 help="""Aggregate function for overlay:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
52 <mean> <median> <mean_j> <median_j>.
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
53 Use mean_j | median_j to keep density overlay but
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
54 aggregate junction counts [default=no aggregation]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
55 parser.add_argument("-C", "--color-factor", type=int, dest="color_factor",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
56 help="Index of column with color levels (1-based)")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
57 parser.add_argument("--alpha", type=float, default=0.5,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
58 help="""Transparency level for density histogram
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
59 [default=%(default)s]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
60 parser.add_argument("-P", "--palette", type=str,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
61 help="""Color palette file. tsv file with >=1 columns,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
62 where the color is the first column""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
63 parser.add_argument("-L", "--labels", type=int, dest="labels", default=1,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
64 help="""Index of column with labels (1-based)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
65 [default=%(default)s]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
66 parser.add_argument("--height", type=float, default=2,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
67 help="""Height of the individual signal plot in inches
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
68 [default=%(default)s]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
69 parser.add_argument("--ann-height", type=float, default=1.5,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
70 dest="ann_height", help="""Height of annotation plot in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
71 inches [default=%(default)s]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
72 parser.add_argument("--width", type=float, default=10,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
73 help="""Width of the plot in inches
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
74 [default=%(default)s]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
75 parser.add_argument("--base-size", type=float, default=14,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
76 dest="base_size", help="""Base font size of the plot in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
77 pch [default=%(default)s]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
78 parser.add_argument("-F", "--out-format", type=str, default="pdf",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
79 dest="out_format", help="""Output file format:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
80 <pdf> <svg> <png> <jpeg> <tiff>
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
81 [default=%(default)s]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
82 parser.add_argument("-R", "--out-resolution", type=int, default=300,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
83 dest="out_resolution", help="""Output file resolution in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
84 PPI (pixels per inch). Applies only to raster output
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
85 formats [default=%(default)s]""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
86 return parser
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
87
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
88
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
89 def parse_coordinates(c):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
90 c = c.replace(",", "")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
91 chr = c.split(":")[0]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
92 start, end = c.split(":")[1].split("-")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
93 # Convert to 0-based
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
94 start, end = int(start) - 1, int(end)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
95 return chr, start, end
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
96
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
97
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
98 def count_operator(CIGAR_op, CIGAR_len, pos, start, end, a, junctions):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
99
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
100 # Match
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
101 if CIGAR_op == "M":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
102 for i in range(pos, pos + CIGAR_len):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
103 if i < start or i >= end:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
104 continue
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
105 ind = i - start
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
106 a[ind] += 1
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
107
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
108 # Insertion or Soft-clip
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
109 if CIGAR_op == "I" or CIGAR_op == "S":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
110 return pos
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
111
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
112 # Deletion
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
113 if CIGAR_op == "D":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
114 pass
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
115
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
116 # Junction
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
117 if CIGAR_op == "N":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
118 don = pos
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
119 acc = pos + CIGAR_len
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
120 if don > start and acc < end:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
121 junctions[(don, acc)] = junctions.setdefault((don, acc), 0) + 1
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
122
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
123 pos = pos + CIGAR_len
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
124
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
125 return pos
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
126
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
127
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
128 def flip_read(s, samflag):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
129 if s == "NONE" or s == "SENSE":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
130 return 0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
131 if s == "ANTISENSE":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
132 return 1
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
133 if s == "MATE1_SENSE":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
134 if int(samflag) & 64:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
135 return 0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
136 if int(samflag) & 128:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
137 return 1
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
138 if s == "MATE2_SENSE":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
139 if int(samflag) & 64:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
140 return 1
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
141 if int(samflag) & 128:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
142 return 0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
143
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
144
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
145 def read_bam(f, c, s):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
146
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
147 _, start, end = parse_coordinates(c)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
148
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
149 # Initialize coverage array and junction dict
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
150 a = {"+": [0] * (end - start)}
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
151 junctions = {"+": OrderedDict()}
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
152 if s != "NONE":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
153 a["-"] = [0] * (end - start)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
154 junctions["-"] = OrderedDict()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
155
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
156 p = sp.Popen("samtools view %s %s " % (f, c), shell=True, stdout=sp.PIPE)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
157 for line in p.communicate()[0].decode('utf8').strip().split("\n"):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
158
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
159 if line == "":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
160 continue
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
161
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
162 line_sp = line.strip().split("\t")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
163 samflag, read_start, CIGAR = line_sp[1], int(line_sp[3]), line_sp[5]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
164
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
165 # Ignore reads with more exotic CIGAR operators
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
166 if any(map(lambda x: x in CIGAR, ["H", "P", "X", "="])):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
167 continue
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
168
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
169 read_strand = ["+", "-"][flip_read(s, samflag) ^ bool(int(samflag) &
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
170 16)]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
171 if s == "NONE":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
172 read_strand = "+"
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
173
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
174 CIGAR_lens = re.split("[MIDNS]", CIGAR)[:-1]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
175 CIGAR_ops = re.split("[0-9]+", CIGAR)[1:]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
176
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
177 pos = read_start
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
178
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
179 for n, CIGAR_op in enumerate(CIGAR_ops):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
180 CIGAR_len = int(CIGAR_lens[n])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
181 pos = count_operator(CIGAR_op, CIGAR_len, pos, start, end,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
182 a[read_strand], junctions[read_strand])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
183
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
184 p.stdout.close()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
185 return a, junctions
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
186
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
187
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
188 def get_bam_path(index, path):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
189 if os.path.isabs(path):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
190 return path
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
191 base_dir = os.path.dirname(index)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
192 return os.path.join(base_dir, path)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
193
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
194
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
195 def read_bam_input(f, overlay, color, label):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
196 if f.endswith(".bam"):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
197 bn = f.strip().split("/")[-1].strip(".bam")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
198 yield bn, f, None, None, bn
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
199 return
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
200 with open(f) as openf:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
201 for line in openf:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
202 line_sp = line.strip().split("\t")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
203 bam = get_bam_path(f, line_sp[1])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
204 overlay_level = line_sp[overlay-1] if overlay else None
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
205 color_level = line_sp[color-1] if color else None
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
206 label_text = line_sp[label-1] if label else None
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
207 yield line_sp[0], bam, overlay_level, color_level, label_text
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
208
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
209
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
210 def prepare_for_R(a, junctions, c, m):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
211
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
212 _, start, _ = parse_coordinates(args.coordinates)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
213
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
214 # Convert the array index to genomic coordinates
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
215 x = list(i+start for i in range(len(a)))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
216 y = a
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
217
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
218 # Arrays for R
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
219 dons, accs, yd, ya, counts = [], [], [], [], []
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
220
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
221 # Prepare arrays for junctions (which will be the arcs)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
222 for (don, acc), n in junctions.items():
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
223
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
224 # Do not add junctions with less than defined coverage
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
225 if n < m:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
226 continue
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
227
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
228 dons.append(don)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
229 accs.append(acc)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
230 counts.append(n)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
231
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
232 yd.append(a[don - start - 1])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
233 ya.append(a[acc - start + 1])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
234
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
235 return x, y, dons, accs, yd, ya, counts
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
236
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
237
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
238 def intersect_introns(data):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
239 data = sorted(data)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
240 it = iter(data)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
241 a, b = next(it)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
242 for c, d in it:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
243 if b > c:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
244 # Use `if b > c` if you want (1,2), (2,3) not to be
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
245 # treated as intersection.
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
246 b = min(b, d)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
247 a = max(a, c)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
248 else:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
249 yield a, b
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
250 a, b = c, d
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
251 yield a, b
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
252
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
253
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
254 def shrink_density(x, y, introns):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
255 new_x, new_y = [], []
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
256 shift = 0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
257 start = 0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
258 # introns are already sorted by coordinates
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
259 for a, b in introns:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
260 end = x.index(a)+1
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
261 new_x += [int(i-shift) for i in x[start:end]]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
262 new_y += y[start:end]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
263 start = x.index(b)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
264 L = (b-a)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
265 shift += L-L**0.7
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
266 new_x += [int(i-shift) for i in x[start:]]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
267 new_y += y[start:]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
268 return new_x, new_y
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
269
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
270
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
271 def shrink_junctions(dons, accs, introns):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
272 new_dons, new_accs = [0]*len(dons), [0]*len(accs)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
273 shift_acc = 0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
274 shift_don = 0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
275 s = set()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
276 junctions = list(zip(dons, accs))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
277 for a, b in introns:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
278 L = b - a
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
279 shift_acc += L-int(L**0.7)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
280 for i, (don, acc) in enumerate(junctions):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
281 if a >= don and b <= acc:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
282 if (don, acc) not in s:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
283 new_dons[i] = don - shift_don
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
284 new_accs[i] = acc - shift_acc
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
285 else:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
286 new_accs[i] = acc - shift_acc
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
287 s.add((don, acc))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
288 shift_don = shift_acc
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
289 return new_dons, new_accs
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
290
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
291
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
292 def read_palette(f):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
293 palette = "#ff0000", "#00ff00", "#0000ff", "#000000"
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
294 if f:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
295 with open(f) as openf:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
296 palette = list(line.split("\t")[0].strip() for line in openf)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
297 return palette
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
298
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
299
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
300 def read_gtf(f, c):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
301 exons = OrderedDict()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
302 transcripts = OrderedDict()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
303 chr, start, end = parse_coordinates(c)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
304 end = end - 1
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
305 with open(f) as openf:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
306 for line in openf:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
307 if line.startswith("#"):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
308 continue
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
309 (el_chr, _, el, el_start, el_end, _,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
310 strand, _, tags) = line.strip().split("\t")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
311 if el_chr != chr:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
312 continue
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
313 d = dict(kv.strip().split(" ") for kv in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
314 tags.strip(";").split("; "))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
315 transcript_id = d["transcript_id"]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
316 el_start, el_end = int(el_start) - 1, int(el_end)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
317 strand = '"' + strand + '"'
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
318 if el == "transcript":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
319 if (el_end > start and el_start < end):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
320 transcripts[transcript_id] = (max(start, el_start),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
321 min(end, el_end),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
322 strand)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
323 continue
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
324 if el == "exon":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
325 if (start < el_start < end or start < el_end < end):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
326 exons.setdefault(transcript_id,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
327 []).append((max(el_start, start),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
328 min(end, el_end), strand))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
329
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
330 return transcripts, exons
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
331
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
332
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
333 def make_introns(transcripts, exons, intersected_introns=None):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
334 new_transcripts = copy.deepcopy(transcripts)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
335 new_exons = copy.deepcopy(exons)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
336 introns = OrderedDict()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
337 if intersected_introns:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
338 for tx, (tx_start, tx_end, strand) in new_transcripts.items():
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
339 total_shift = 0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
340 for a, b in intersected_introns:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
341 L = b - a
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
342 shift = L - int(L**0.7)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
343 total_shift += shift
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
344 for i, (exon_start, exon_end, strand) in \
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
345 enumerate(exons.get(tx, [])):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
346 new_exon_start, new_exon_end = new_exons[tx][i][:2]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
347 if a < exon_start:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
348 if b > exon_end:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
349 if i == len(exons[tx])-1:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
350 total_shift = total_shift - shift + \
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
351 (exon_start - a)*(1-int(L**-0.3))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
352 shift = (exon_start - a)*(1-int(L**-0.3))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
353 new_exon_end = new_exons[tx][i][1] - shift
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
354 new_exon_start = new_exons[tx][i][0] - shift
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
355 if b <= exon_end:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
356 new_exon_end = new_exons[tx][i][1] - shift
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
357 new_exons[tx][i] = (new_exon_start, new_exon_end, strand)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
358 tx_start = min(tx_start,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
359 sorted(new_exons.get(tx, [[sys.maxsize]]))[0][0])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
360 new_transcripts[tx] = (tx_start, tx_end - total_shift, strand)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
361
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
362 for tx, (tx_start, tx_end, strand) in new_transcripts.items():
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
363 intron_start = tx_start
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
364 ex_end = 0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
365 for ex_start, ex_end, strand in sorted(new_exons.get(tx, [])):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
366 intron_end = ex_start
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
367 if tx_start < ex_start:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
368 introns.setdefault(tx, []).append((intron_start, intron_end,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
369 strand))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
370 intron_start = ex_end
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
371 if tx_end > ex_end:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
372 introns.setdefault(tx, []).append((intron_start, tx_end, strand))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
373 d = {'transcripts': new_transcripts,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
374 'exons': new_exons,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
375 'introns': introns}
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
376 return d
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
377
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
378
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
379 def gtf_for_ggplot(annotation, start, end, arrow_bins):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
380 arrow_space = int((end - start)/arrow_bins)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
381 s = """
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
382
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
383 # data table with exons
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
384 ann_list = list(
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
385 "exons" = data.table(),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
386 "introns" = data.table()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
387 )
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
388 """
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
389
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
390 if annotation["exons"]:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
391
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
392 s += """
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
393 ann_list[['exons']] = data.table(
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
394 tx = rep(c(%(tx_exons)s), c(%(n_exons)s)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
395 start = c(%(exon_start)s),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
396 end = c(%(exon_end)s),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
397 strand = c(%(strand)s)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
398 )
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
399 """ % ({
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
400 "tx_exons": ",".join(annotation["exons"].keys()),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
401 "n_exons": ",".join(map(str, map(len,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
402 annotation["exons"].values()))),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
403 "exon_start": ",".join(map(str, (v[0] for vs in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
404 annotation["exons"].values() for v in vs))),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
405 "exon_end": ",".join(map(str, (v[1] for vs in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
406 annotation["exons"].values() for v in vs))),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
407 "strand": ",".join(map(str, (v[2] for vs in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
408 annotation["exons"].values() for v in vs))),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
409 })
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
410
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
411 if annotation["introns"]:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
412
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
413 s += """
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
414 ann_list[['introns']] = data.table(
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
415 tx = rep(c(%(tx_introns)s), c(%(n_introns)s)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
416 start = c(%(intron_start)s),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
417 end = c(%(intron_end)s),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
418 strand = c(%(strand)s)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
419 )
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
420 # Create data table for strand arrows
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
421 txarrows = data.table()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
422 introns = ann_list[['introns']]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
423 # Add right-pointing arrows for plus strand
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
424 if ("+" %%in%% introns$strand) {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
425 txarrows = rbind(
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
426 txarrows,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
427 introns[strand=="+" & end-start>5, list(
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
428 seq(start+4,end,by=%(arrow_space)s)-1,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
429 seq(start+4,end,by=%(arrow_space)s)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
430 ), by=.(tx,start,end)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
431 ]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
432 )
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
433 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
434 # Add left-pointing arrows for minus strand
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
435 if ("-" %%in%% introns$strand) {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
436 txarrows = rbind(txarrows,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
437 introns[strand=="-" & end-start>5,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
438 list(seq(start,max(start+1, end-4),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
439 by=%(arrow_space)s),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
440 seq(start,max(start+1, end-4),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
441 by=%(arrow_space)s)-1
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
442 ),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
443 by=.(tx,start,end)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
444 ]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
445 )
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
446 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
447 """ % ({
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
448 "tx_introns": ",".join(annotation["introns"].keys()),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
449 "n_introns": ",".join(map(str, map(len,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
450 annotation["introns"].values()))),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
451 "intron_start": ",".join(map(str, (v[0] for vs in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
452 annotation["introns"].values() for v in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
453 vs))),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
454 "intron_end": ",".join(map(str, (v[1] for vs in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
455 annotation["introns"].values() for v in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
456 vs))),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
457 "strand": ",".join(map(str, (v[2] for vs in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
458 annotation["introns"].values() for v in vs))),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
459 "arrow_space": arrow_space,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
460 })
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
461
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
462 s += """
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
463
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
464 gtfp = ggplot()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
465 if (length(ann_list[['introns']]) > 0) {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
466 gtfp = gtfp + geom_segment(data = ann_list[['introns']],
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
467 aes(x = start,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
468 xend = end,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
469 y = tx,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
470 yend = tx),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
471 size = 0.3)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
472 gtfp = gtfp + geom_segment(data = txarrows,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
473 aes(x = V1,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
474 xend = V2,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
475 y = tx,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
476 yend = tx),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
477 arrow = arrow(length = unit(0.02, "npc")))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
478 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
479 if (length(ann_list[['exons']]) > 0) {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
480 gtfp = gtfp + geom_segment(data = ann_list[['exons']],
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
481 aes(x = start,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
482 xend = end,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
483 y = tx,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
484 yend = tx),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
485 size = 5,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
486 alpha = 1)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
487 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
488 gtfp = gtfp + scale_y_discrete(expand = c(0, 0.5))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
489 gtfp = gtfp + scale_x_continuous(expand = c(0, 0.25),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
490 limits = c( %s,% s))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
491 gtfp = gtfp + labs(y = NULL)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
492 gtfp = gtfp + theme(axis.line = element_blank(),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
493 axis.text.x = element_blank(),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
494 axis.ticks = element_blank())
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
495 """ % (start, end)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
496
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
497 return s
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
498
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
499
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
500 def setup_R_script(h, w, b, label_dict):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
501 s = """
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
502 library(ggplot2)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
503 library(grid)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
504 library(gridExtra)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
505 library(data.table)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
506 library(gtable)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
507
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
508 scale_lwd = function(r) {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
509 lmin = 0.1
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
510 lmax = 4
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
511 return( r*(lmax-lmin)+lmin )
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
512 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
513
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
514 base_size = %(b)s
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
515 height = ( %(h)s + base_size*0.352777778/67 ) * 1.02
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
516 width = %(w)s
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
517 theme_set(theme_bw(base_size=base_size))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
518 theme_update(
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
519 plot.margin = unit(c(15,15,15,15), "pt"),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
520 panel.grid = element_blank(),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
521 panel.border = element_blank(),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
522 axis.line = element_line(size=0.5),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
523 axis.title.x = element_blank(),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
524 axis.title.y = element_text(angle=0, vjust=0.5)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
525 )
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
526
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
527 labels = list(%(labels)s)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
528
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
529 density_list = list()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
530 junction_list = list()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
531
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
532 """ % ({
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
533 'h': h,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
534 'w': w,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
535 'b': b,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
536 'labels': ",".join(('"%s"="%s"' % (id, lab) for id, lab in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
537 label_dict.items())),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
538 })
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
539 return s
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
540
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
541
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
542 def median(lst):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
543 quotient, remainder = divmod(len(lst), 2)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
544 if remainder:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
545 return sorted(lst)[quotient]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
546 return sum(sorted(lst)[quotient - 1:quotient + 1]) / 2.
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
547
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
548
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
549 def mean(lst):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
550 return sum(lst)/len(lst)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
551
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
552
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
553 def make_R_lists(id_list, d, overlay_dict, aggr, intersected_introns):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
554 s = ""
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
555 aggr_f = {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
556 "mean": mean,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
557 "median": median,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
558 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
559 id_list = id_list if not overlay_dict else overlay_dict.keys()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
560 # Iterate over ids to get bam signal and junctions
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
561 for k in id_list:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
562 x, y, dons, accs, yd, ya, counts = [], [], [], [], [], [], []
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
563 if not overlay_dict:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
564 x, y, dons, accs, yd, ya, counts = d[k]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
565 if intersected_introns:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
566 x, y = shrink_density(x, y, intersected_introns)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
567 dons, accs = shrink_junctions(dons, accs, intersected_introns)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
568 else:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
569 for id in overlay_dict[k]:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
570 xid, yid, donsid, accsid, ydid, yaid, countsid = d[id]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
571 if intersected_introns:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
572 xid, yid = shrink_density(xid, yid, intersected_introns)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
573 donsid, accsid = shrink_junctions(donsid, accsid,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
574 intersected_introns)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
575 x += xid
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
576 y += yid
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
577 dons += donsid
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
578 accs += accsid
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
579 yd += ydid
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
580 ya += yaid
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
581 counts += countsid
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
582 if aggr and "_j" not in aggr:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
583 x = d[overlay_dict[k][0]][0]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
584 y = list(map(aggr_f[aggr], zip(*(d[id][1] for id in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
585 overlay_dict[k]))))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
586 if intersected_introns:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
587 x, y = shrink_density(x, y, intersected_introns)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
588 # dons, accs, yd, ya, counts = [], [], [], [], []
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
589 s += """
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
590 density_list[["%(id)s"]] = data.frame(x = c(%(x)s), y = c(%(y)s))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
591 junction_list[["%(id)s"]] = data.frame(x = c(%(dons)s),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
592 xend=c(%(accs)s),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
593 y=c(%(yd)s),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
594 yend=c(%(ya)s),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
595 count=c(%(counts)s))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
596 """ % ({
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
597 "id": k,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
598 'x': ",".join(map(str, x)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
599 'y': ",".join(map(str, y)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
600 'dons': ",".join(map(str, dons)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
601 'accs': ",".join(map(str, accs)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
602 'yd': ",".join(map(str, yd)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
603 'ya': ",".join(map(str, ya)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
604 'counts': ",".join(map(str, counts)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
605 })
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
606 return s
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
607
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
608
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
609 def plot(R_script):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
610 p = sp.Popen("R --vanilla --slave", shell=True, stdin=sp.PIPE)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
611 p.communicate(input=R_script.encode())
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
612 p.stdin.close()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
613 p.wait()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
614 return
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
615
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
616
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
617 def colorize(d, p, color_factor):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
618 levels = sorted(set(d.values()))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
619 n = len(levels)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
620 if n > len(p):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
621 p = (p*n)[:n]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
622 if color_factor:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
623 s = "color_list = list(%s)\n" % (",".join('%s="%s"' % (k,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
624 p[levels.index(v)]) for k, v in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
625 d.items()))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
626 else:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
627 s = "color_list = list(%s)\n" % (",".join('%s="%s"' % (k, "grey") for
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
628 k, v in d.items()))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
629 return s
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
630
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
631
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
632 if __name__ == "__main__":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
633
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
634 strand_dict = {"plus": "+", "minus": "-"}
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
635
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
636 parser = define_options()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
637 if len(sys.argv) == 1:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
638 parser.print_help()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
639 sys.exit(1)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
640 args = parser.parse_args()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
641
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
642 if args.aggr and not args.overlay:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
643 print("""ERROR: Cannot apply aggregate function
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
644 if overlay is not selected.""")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
645 exit(1)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
646
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
647 palette = read_palette(args.palette)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
648
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
649 (bam_dict, overlay_dict, color_dict,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
650 id_list, label_dict) = ({"+": OrderedDict()}, OrderedDict(),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
651 OrderedDict(), [], OrderedDict())
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
652 if args.strand != "NONE":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
653 bam_dict["-"] = OrderedDict()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
654 if args.junctions_bed != "":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
655 junctions_list = []
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
656
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
657 for (id, bam, overlay_level,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
658 color_level, label_text) in read_bam_input(args.bam,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
659 args.overlay,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
660 args.color_factor,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
661 args.labels):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
662 if not os.path.isfile(bam):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
663 continue
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
664 id_list.append(id)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
665 label_dict[id] = label_text
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
666 a, junctions = read_bam(bam, args.coordinates, args.strand)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
667 if a.keys() == ["+"] and all(map(lambda x: x == 0,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
668 list(a.values()[0]))):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
669 print("ERROR: No reads in the specified area.")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
670 exit(1)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
671 for strand in a:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
672 # Store junction information
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
673 if args.junctions_bed:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
674 for k, v in zip(junctions[strand].keys(),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
675 junctions[strand].values()):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
676 if v > args.min_coverage:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
677 junctions_list.append('\t'.join([args.coordinates.split
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
678 (':')[0], str(k[0]), str(k[1]),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
679 id, str(v), strand]))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
680 bam_dict[strand][id] = prepare_for_R(a[strand],
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
681 junctions[strand],
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
682 args.coordinates,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
683 args.min_coverage)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
684 if color_level is None:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
685 color_dict.setdefault(id, id)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
686 if overlay_level is not None:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
687 overlay_dict.setdefault(overlay_level, []).append(id)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
688 label_dict[overlay_level] = overlay_level
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
689 color_dict.setdefault(overlay_level, overlay_level)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
690 if overlay_level is None:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
691 color_dict.setdefault(id, color_level)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
692
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
693 # No bam files
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
694 if not bam_dict["+"]:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
695 print("ERROR: No available bam files.")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
696 exit(1)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
697
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
698 # Write junctions to BED
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
699 if args.junctions_bed:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
700 if not args.junctions_bed.endswith('.bed'):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
701 args.junctions_bed = args.junctions_bed + '.bed'
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
702 jbed = open(args.junctions_bed, 'w')
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
703 jbed.write('\n'.join(sorted(junctions_list)))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
704 jbed.close()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
705
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
706 if args.gtf:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
707 transcripts, exons = read_gtf(args.gtf, args.coordinates)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
708
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
709 if args.out_format not in ('pdf', 'png', 'svg', 'tiff', 'jpeg'):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
710 print("""ERROR: Provided output format '%s' is not available.
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
711 Please select among 'pdf', 'png', 'svg',
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
712 'tiff' or 'jpeg'""" % args.out_format)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
713 exit(1)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
714
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
715 # Iterate for plus and minus strand
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
716 for strand in bam_dict:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
717
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
718 # Output file name (allow tiff/tif and jpeg/jpg extensions)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
719 if args.out_prefix.endswith(('.pdf', '.png', '.svg', '.tiff',
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
720 '.tif', '.jpeg', '.jpg')):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
721 out_split = os.path.splitext(args.out_prefix)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
722 if (args.out_format == out_split[1][1:] or
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
723 args.out_format == 'tiff'
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
724 and out_split[1] in ('.tiff', '.tif') or
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
725 args.out_format == 'jpeg'
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
726 and out_split[1] in ('.jpeg', '.jpg')):
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
727 args.out_prefix = out_split[0]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
728 out_suffix = out_split[1][1:]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
729 else:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
730 out_suffix = args.out_format
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
731 else:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
732 out_suffix = args.out_format
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
733 out_prefix = args.out_prefix + "_" + strand
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
734 if args.strand == "NONE":
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
735 out_prefix = args.out_prefix
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
736 else:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
737 if args.out_strand != "both" \
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
738 and strand != strand_dict[args.out_strand]:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
739 continue
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
740
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
741 # Find set of junctions to perform shrink
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
742 intersected_introns = None
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
743 if args.shrink:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
744 introns = (v for vs in bam_dict[strand].values() for v in
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
745 zip(vs[2], vs[3]))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
746 intersected_introns = list(intersect_introns(introns))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
747
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
748 # *** PLOT *** Define plot height
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
749 bam_height = args.height * len(id_list)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
750 if args.overlay:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
751 bam_height = args.height * len(overlay_dict)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
752 if args.gtf:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
753 bam_height += args.ann_height
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
754
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
755 # *** PLOT *** Start R script by loading libraries,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
756 # initializing variables, etc...
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
757 R_script = setup_R_script(bam_height, args.width,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
758 args.base_size, label_dict)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
759
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
760 R_script += colorize(color_dict, palette, args.color_factor)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
761
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
762 # *** PLOT *** Prepare annotation plot only for the first bam file
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
763 arrow_bins = 50
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
764 if args.gtf:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
765 # Make introns from annotation (they are shrunk if required)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
766 annotation = make_introns(transcripts, exons, intersected_introns)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
767 x = list(bam_dict[strand].values())[0][0]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
768 if args.shrink:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
769 x, _ = shrink_density(x, x, intersected_introns)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
770 R_script += gtf_for_ggplot(annotation, x[0], x[-1], arrow_bins)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
771
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
772 R_script += make_R_lists(id_list, bam_dict[strand], overlay_dict,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
773 args.aggr, intersected_introns)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
774
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
775 R_script += """
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
776
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
777 pdf(NULL) # just to remove the blank pdf produced by ggplotGrob
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
778 # fix problems with ggplot2 vs >3.0.0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
779 if(packageVersion('ggplot2') >= '3.0.0'){
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
780 vs = 1
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
781 } else {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
782 vs = 0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
783 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
784
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
785 density_grobs = list();
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
786
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
787 for (bam_index in 1:length(density_list)) {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
788
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
789 id = names(density_list)[bam_index]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
790 d = data.table(density_list[[id]])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
791 junctions = data.table(junction_list[[id]])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
792
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
793 maxheight = max(d[['y']])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
794
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
795 # Density plot
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
796 gp = ggplot(d) + geom_bar(aes(x, y), width=1,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
797 position='identity',
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
798 stat='identity',
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
799 fill=color_list[[id]],
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
800 alpha=%(alpha)s)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
801 gp = gp + labs(y=labels[[id]])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
802 gp = gp + scale_x_continuous(expand=c(0,0.2))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
803
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
804 # fix problems with ggplot2 vs >3.0.0
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
805 if(packageVersion('ggplot2') >= '3.0.0') {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
806 gp = gp +
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
807 scale_y_continuous(breaks =
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
808 ggplot_build(gp
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
809 )$layout$panel_params[[1]]$y.major_source)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
810 } else {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
811 gp = gp +
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
812 scale_y_continuous(breaks =
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
813 ggplot_build(gp
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
814 )$layout$panel_ranges[[1]]$y.major_source)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
815 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
816
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
817 # Aggregate junction counts
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
818 row_i = c()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
819 if (nrow(junctions) >0 ) {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
820
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
821 junctions$jlabel = as.character(junctions$count)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
822 junctions = setNames(junctions[,.(max(y),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
823 max(yend),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
824 round(mean(count)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
825 paste(jlabel,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
826 collapse=",")),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
827 keyby=.(x,xend)],
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
828 names(junctions))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
829 if ("%(args.aggr)s" != "") {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
830 junctions = setNames(
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
831 junctions[,.(max(y),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
832 max(yend),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
833 round(%(args.aggr)s(count)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
834 round(%(args.aggr)s(count))),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
835 keyby=.(x,xend)],
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
836 names(junctions))
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
837 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
838 # The number of rows (unique junctions per bam) has to be
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
839 # calculated after aggregation
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
840 row_i = 1:nrow(junctions)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
841 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
842
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
843
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
844 for (i in row_i) {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
845
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
846 j_tot_counts = sum(junctions[['count']])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
847
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
848 j = as.numeric(junctions[i,1:5])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
849
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
850 if ("%(args.aggr)s" != "") {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
851 j[3] = as.numeric(d[x==j[1]-1,y])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
852 j[4] = as.numeric(d[x==j[2]+1,y])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
853 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
854
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
855 # Find intron midpoint
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
856 xmid = round(mean(j[1:2]), 1)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
857 ymid = max(j[3:4]) * 1.2
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
858
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
859 # Thickness of the arch
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
860 lwd = scale_lwd(j[5]/j_tot_counts)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
861
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
862 curve_par = gpar(lwd=lwd, col=color_list[[id]])
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
863
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
864 # Arc grobs
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
865
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
866 # Choose position of the arch (top or bottom)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
867 nss = i
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
868 if (nss%%%%2 == 0) { #bottom
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
869 ymid = -0.3 * maxheight
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
870 # Draw the arcs
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
871 # Left
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
872 curve = xsplineGrob(x = c(0, 0, 1, 1),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
873 y = c(1, 0, 0, 0),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
874 shape = 1,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
875 gp = curve_par)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
876 gp = gp + annotation_custom(grob = curve,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
877 j[1], xmid, 0, ymid)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
878 # Right
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
879 curve = xsplineGrob(x = c(1, 1, 0, 0),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
880 y = c(1, 0, 0, 0),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
881 shape = 1,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
882 gp = curve_par)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
883 gp = gp + annotation_custom(grob = curve,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
884 xmid,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
885 j[2],
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
886 0, ymid)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
887 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
888
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
889 if (nss%%%%2 != 0) { #top
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
890 # Draw the arcs
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
891 # Left
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
892 curve = xsplineGrob(x = c(0, 0, 1, 1),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
893 y = c(0, 1, 1, 1),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
894 shape = 1,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
895 gp = curve_par)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
896 gp = gp + annotation_custom(grob = curve,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
897 j[1], xmid, j[3], ymid)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
898 # Right
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
899 curve = xsplineGrob(x = c(1, 1, 0, 0),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
900 y = c(0, 1, 1, 1),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
901 shape = 1,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
902 gp = curve_par)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
903 gp = gp + annotation_custom(grob = curve,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
904 xmid, j[2], j[4], ymid)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
905 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
906
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
907 # Add junction labels
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
908 gp = gp + annotate("label", x = xmid, y = ymid,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
909 label = as.character(junctions[i,6]),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
910 vjust = 0.5, hjust = 0.5,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
911 label.padding = unit(0.01, "lines"),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
912 label.size = NA,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
913 size = (base_size*0.352777778)*0.6)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
914
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
915 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
916
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
917 gpGrob = ggplotGrob(gp);
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
918 gpGrob$layout$clip[gpGrob$layout$name=="panel"] <- "off"
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
919 if (bam_index == 1) {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
920 # fix problems ggplot2 vs
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
921 maxWidth = gpGrob$widths[2+vs] + gpGrob$widths[3+vs];
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
922 maxYtextWidth = gpGrob$widths[3+vs];
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
923 # Extract x axis grob (trim=F --> keep empty cells)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
924 xaxisGrob <- gtable_filter(gpGrob, "axis-b", trim=F)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
925 # fix problems ggplot2 vs
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
926 xaxisGrob$heights[8+vs] = gpGrob$heights[1]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
927 x.axis.height = gpGrob$heights[7+vs] + gpGrob$heights[1]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
928 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
929
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
930
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
931 # Remove x axis from all density plots
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
932 kept_names = gpGrob$layout$name[gpGrob$layout$name != "axis-b"]
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
933 gpGrob <- gtable_filter(gpGrob,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
934 paste(kept_names, sep = "",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
935 collapse = "|"),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
936 trim=F)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
937
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
938 # Find max width of y text and y label and max width of y text
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
939 # fix problems ggplot2 vs
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
940 maxWidth = grid::unit.pmax(maxWidth,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
941 gpGrob$widths[2+vs] +
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
942 gpGrob$widths[3+vs]);
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
943 maxYtextWidth = grid::unit.pmax(maxYtextWidth,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
944 gpGrob$widths[3+vs]);
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
945 density_grobs[[id]] = gpGrob;
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
946 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
947
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
948 # Add x axis grob after density grobs BEFORE annotation grob
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
949 density_grobs[["xaxis"]] = xaxisGrob
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
950
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
951 # Annotation grob
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
952 if (%(args.gtf)s == 1) {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
953 gtfGrob = ggplotGrob(gtfp);
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
954 maxWidth = grid::unit.pmax(maxWidth,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
955 gtfGrob$widths[2+vs] +
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
956 gtfGrob$widths[3+vs]);
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
957 density_grobs[['gtf']] = gtfGrob;
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
958 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
959
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
960 # Reassign grob widths to align the plots
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
961 for (id in names(density_grobs)) {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
962 density_grobs[[id]]$widths[1] <-
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
963 density_grobs[[id]]$widths[1] +
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
964 maxWidth - (density_grobs[[id]]$widths[2 + vs] +
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
965 maxYtextWidth)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
966 # fix problems ggplot2 vs
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
967 density_grobs[[id]]$widths[3 + vs] <-
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
968 maxYtextWidth # fix problems ggplot2 vs
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
969 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
970
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
971 # Heights for density, x axis and annotation
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
972 heights = unit.c(
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
973 unit(rep(%(signal_height)s,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
974 length(density_list)), "in"),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
975 x.axis.height,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
976 unit(%(ann_height)s*%(args.gtf)s, "in")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
977 )
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
978
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
979 # Arrange grobs
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
980 argrobs = arrangeGrob(
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
981 grobs=density_grobs,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
982 ncol=1,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
983 heights = heights,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
984 );
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
985
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
986 # Save plot to file in the requested format
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
987 if ("%(out_format)s" == "tiff"){
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
988 # TIFF images will be lzw-compressed
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
989 ggsave("%(out)s",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
990 plot = argrobs,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
991 device = "tiff",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
992 width = width,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
993 height = height,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
994 units = "in",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
995 dpi = %(out_resolution)s,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
996 compression = "lzw")
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
997 } else {
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
998 ggsave("%(out)s",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
999 plot = argrobs,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1000 device = "%(out_format)s",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1001 width = width,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1002 height = height,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1003 units = "in",
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1004 dpi = %(out_resolution)s)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1005 }
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1006
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1007 dev.log = dev.off()
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1008
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1009 """ % ({
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1010 "out": "%s.%s" % (out_prefix, out_suffix),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1011 "out_format": args.out_format,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1012 "out_resolution": args.out_resolution,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1013 "args.gtf": float(bool(args.gtf)),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1014 "args.aggr": args.aggr.rstrip("_j"),
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1015 "signal_height": args.height,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1016 "ann_height": args.ann_height,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1017 "alpha": args.alpha,
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1018 })
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1019 if os.getenv('GGSASHIMI_DEBUG') is not None:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1020 with open("R_script", 'w') as r:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1021 r.write(R_script)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1022 else:
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1023 plot(R_script)
9304dd9a16a2 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sashimi_plot commit 746c03a1187e1d708af8628920a0c615cddcdacc"
artbio
parents:
diff changeset
1024 exit()