Mercurial > repos > yhoogstrate > featurecounts

annotate featurecounts2bed.sh @ 0:ab0323782bb5 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
author yhoogstrate
date Mon, 18 May 2015 04:49:19 -0400
parents
children b5c93611d2c5
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
1 #!/bin/bash
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
2
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
3 # featurecounts2bed - converts featureCounts output to BED format
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
4
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
5 # Copyright 2013-2014, Youri Hoogstrate
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
6
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
7 # This program is free software: you can redistribute it and/or modify
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
8 # it under the terms of the GNU General Public License as published by
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
9 # the Free Software Foundation, either version 3 of the License, or
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
10 # (at your option) any later version.
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
11
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
12 # This program is distributed in the hope that it will be useful,
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
13 # but WITHOUT ANY WARRANTY; without even the implied warranty of
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
14 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
15 # GNU General Public License at <http://www.gnu.org/licenses/> for
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
16 # more details.
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
17
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
18 # This tool has been written by Youri Hoogstrate from the Erasmus
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
19 # Medical Center (Rotterdam, Netherlands) on behalf of the Translational
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
20 # Research IT (TraIT) project:
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
21 # http://www.ctmm.nl/en/programmas/infrastructuren/traitprojecttranslationeleresearch
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
22 #
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
23 # More tools by the Translational Research IT (TraIT) project can be
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
24 # found in the following repository:
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
25 # http://toolshed.dtls.nl/
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
26
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
27
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
28 exon_level="true"
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
29 filename=""
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
30
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
31 # Parse parameters
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
32 while getopts e:f: option
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
33 do
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
34 case "${option}"
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
35 in
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
36 e) exon_level=${OPTARG};;
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
37 f) filename=$OPTARG;;
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
38 esac
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
39 done
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
40
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
41 # Convert the file
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
42 if [ $filename == "" ]; then
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
43 echo "Usage:"
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
44 echo " -e [true, false] true = entry for every exon; false = line for genes first exon"
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
45 echo " -f FILENAME from featureCounts"
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
46 else
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
47 while read line; do
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
48 first=${line:0:1}
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
49 if [ $first != "#" ]; then
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
50 columns=($line)
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
51 uid=${columns[@]:0:1}
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
52 if [ $uid != "Geneid" ]; then
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
53 chr=${columns[@]:1:1}
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
54 start=${columns[@]:2:1}
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
55 stop=${columns[@]:3:1}
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
56 direction=${columns[@]:4:1}
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
57 length=${columns[@]:5:1}
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
58 count=${columns[@]:6:1}
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
59
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
60 chr_splitted=($(echo $chr | tr ";" "\n"))
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
61 start_splitted=($(echo $start | tr ";" "\n"))
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
62 stop_splitted=($(echo $stop | tr ";" "\n"))
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
63 strand_splitted=($(echo $direction | tr ";" "\n"))
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
64
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
65 if [ $exon_level == "true" ]; then
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
66 n=${#chr_splitted[@]}
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
67 else
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
68 n=1
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
69 fi
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
70
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
71 for (( i=0; i<$n; i++ ))
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
72 do
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
73 echo ${chr_splitted[@]:$i:1}" "${start_splitted[@]:$i:1}" "${stop_splitted[@]:$i:1}" "$uid" ("$((${stop_splitted[@]:$i:1}-${start_splitted[@]:$i:1}))"/"$length"nt) "$count" "${strand_splitted[@]:$i:1}
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
74 done
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
75 fi
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
76 fi
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
77 done < $filename
ab0323782bb5 planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/featurecounts commit cc900436bad9c6cca1f73d438c1f158d3bfc4318-dirty
yhoogstrate
parents:
diff changeset
78 fi