Mercurial > repos > alermine > nebula
comparison [APliBio]Nebula tools suite/Nebula/CreateControlSubSet/createControlSubSet_wrapper.sh @ 0:2ec3ba0e9e70 draft
Uploaded
| author | alermine |
|---|---|
| date | Thu, 25 Oct 2012 08:18:25 -0400 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| -1:000000000000 | 0:2ec3ba0e9e70 |
|---|---|
| 1 #:t:::::::::::::::::g@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ | |
| 2 #:t::::::::::::::;@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ | |
| 3 #:::::::::::::z;@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ | |
| 4 #::::::::::::i@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ | |
| 5 #::::::::::::@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@$@@@@ | |
| 6 #:::::::::::3@@@@@@@@@@@@@@@@@@@@@@@@@B@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ | |
| 7 #::::::::::3@@@@@@@@@@@@@@@@@@@@@BEEESSE5EEEEBBM@@@@@@@@@@@@@@@@@@@@@@@@@@ | |
| 8 #::::::::::3@@@@@@@@@@@@@@@@@@@@BEEEEEE35EE55E2355E5SBMB@@@@@@@@@@@@@@@@@$ | |
| 9 #::::::::::@@@@@@@@@@@@@@@@@@@EEEE55533t3tttt::::::!!!!7755E755SBBMMM@@@MM | |
| 10 #::::::::::3@@@@@@@@@@@@@@@@@@EEEE2t3ttttt:::::::::::::::::::::::!7?5225EE | |
| 11 #::::::::::3@@@@@@@@@@@@@@@@@@EEEEE31t::::::::::::::::::::::::::::::::3E5@ | |
| 12 #::::::::::3@@@@@@@@@@@@@@@@@@EEEEEEtt:::::::::::::::::::::::::::::::::353 | |
| 13 #::::::::::3@@@@@@@@@@@@@@@@@@EEEEEE1ttz::::::::::::::::::::::::::::::::35 | |
| 14 #:::::::::::@@@@@@@@@@@@@@@@@@EEEEEEEtz1::::::::::::::::::::::::::::::::t: | |
| 15 #:::::::::!3@@@@@@@@@@@@@@@@@@@EEEEEttt::::::::::::::::::::::::::::::::;zz | |
| 16 #::::::::::@@@@@@@@@@@@@@@@@@@@EEEEEttt:::::z;z:::::::::::::::::::::::::13 | |
| 17 #::::::::::3B@@@@@@@@@@@@@@@@@@EEEEEEE3tt:czzztti;:::::::::::::::::::::::3 | |
| 18 #::::ttt::::3@@@@@@@@@@@@@@@@EEEEE5EE25Ezt1EEEz5Etzzz;;;;::::::::::::::::: | |
| 19 #:::::::::::I9@@@@@@@@@@@@@@@@@@@@@@@@@@EEEEEE@@@@@@@@@@@@@@Ez;::::::::::: | |
| 20 #:::::::::::::E@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@Ez:::::: | |
| 21 #::::::::::::::E@@@@@@@@@@@@@@@@@@@@@@@@@@@@@BE5EBB@@@@@@@@@@@@@@@EEE::::: | |
| 22 #:::::::::::::::@@@@@@@@@@@@@@@@@@@@@@@@@@@@E1::35@@@@@@@@@@ME3MMME2:::::: | |
| 23 #:::::::::::::::?@@@@@@@@@@@@@@@@@@M@@@@@@@EE:::::3SB@@BBESEEt:::::::::::: | |
| 24 #::::::::::::::::J$@@@@@@@B@@@@@@@@@@@@@@@@EE:::::::!35E33t::::::::::::::: | |
| 25 #:::::::::::::::::3@E@@@EE5EESE5EESE@@@@@@@Et::::::::::::tz::::::::::::::: | |
| 26 #:::::::::::::::::J@E$@EEE5133555SE@@@@@@@@Et::::::::::::::::::::::::::::: | |
| 27 #::::::::::::::::::E@E@EEEEtt3523EEE@@@@@@@E:::::::::::::::::::::::::::::: | |
| 28 #:t::::::::::::::::JEE3@@@EEEEEEEEEE@@@@@@@E:::::::::t;::::::::::::::::::: | |
| 29 #:t:::::::::::::::::!5ES@EEEEEEEEES@@@@@@@@@E;:::;;;:3Ez:::::::::::::::::: | |
| 30 #:t::::::::::::::::::::JE@@EEEEEEE@@@@@@@@@@@@@@@@ME!:::;::::::::::::::::: | |
| 31 #:tz::::::::::::::::::::JE@@@EEEE@@@@@@@@@@@@@@EE!:::::::t:::::::::::::::: | |
| 32 #:t::::::::::::::::::::::3@@@@@@@@@@@@@@@@@@ESBE:::::::::::::::::::::::::: | |
| 33 #:::::::::::::::::::::::::Q@@@@@@@@@@@@@@@@EE3EE;:::::zzzz:::::::::::::::: | |
| 34 #:::::::::::::::::::::::::3@@@@@@@@@@@@@@@@@@@@@@NN@@@@@@Ez::::::::::::::: | |
| 35 #:zt:::::::::::::::::::::::3@@@@EE@@@@@@@@@@EEEEt::;z113E5t::::::::::::::: | |
| 36 #::tt:::::::::::::::::::::::3@@@E@@@@@@@@@@@@@@@@BEt::::::::::::::::t::::: | |
| 37 #:tt:t:::::::::::::::::::::::?S@@@@@@@@@@@BBEEE51!::::::::::::::zzzEt::::: | |
| 38 #::::::::::::::::::::::::::::::3Q@@@@@@@BEEEEEt:::::::::::::;zz@@@EE:::::: | |
| 39 #::::::::::::::::::::::::::::::::75B@@@@@EEEtt;:::::::::;zz@@@@BEEEtz::::: | |
| 40 #::::::::::::::::::::::::::::::::::::?9@@@@@@@@@@@E2Ezg@@@@@B@@@EEEE1t:::: | |
| 41 #:::::::::::::::::::::::::::::::::::::::3@@@@@@@@@@@@@@@@@@@E@EEEEEEEzzz:: | |
| 42 #::::::::::::::::::::::::::::::::::::;@@@@@@@@@@@@@@@@@@@@@@@EEEEEEE5ttttt | |
| 43 #:::::::::::::::::::::::::::::::;g@@@@@@@@@@@@@@@@@@@@@@@@@@EEEEEEEEEEEtzt | |
| 44 #::::::::::::::::::::::::::::;@@@@@@@@@@@@@@@@@@@@@@@@@@E@@EEEEEEEEEEEE@@@ | |
| 45 #::::::::::::::::::::::::::g@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@EEEE3EEEE@@@@@@@ | |
| 46 #:::::::::::::::::::::;;g@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@EEEt33@@@@@@@@@@ | |
| 47 #:::::::::::::::::;g@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@E@@@@@@EEEtg@@@@@@@@@@@@ | |
| 48 #::::::::::::::;@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@EEEE@@@@@@@@@@@@@@@@@@@@@@@@ | |
| 49 #:::::::::::::@@@@@@@@@@@@@@@@@$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ | |
| 50 #::::::::::;@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ | |
| 51 # | |
| 52 # Copyleft ↄ⃝ 2012 Institut Curie | |
| 53 # Author(s): Valentina Boeva, Alban Lermine (Institut Curie) 2012 | |
| 54 # Contact: valentina.boeva@curie.fr, alban.lermine@curie.fr | |
| 55 # This software is distributed under the terms of the GNU General | |
| 56 # Public License, either Version 2, June 1991 or Version 3, June 2007. | |
| 57 | |
| 58 #!/bin/bash | |
| 59 | |
| 60 CONTROLFILE="" | |
| 61 | |
| 62 while getopts "f:c:t:o:s:q:p:" optionName; do | |
| 63 case "$optionName" in | |
| 64 | |
| 65 f) CHIPFILE="$OPTARG";; | |
| 66 c) CONTROLFILE="$OPTARG";; | |
| 67 t) TYPE="$OPTARG";; | |
| 68 s) CHIPOUTPUT="$OPTARG";; | |
| 69 o) CONTROLOUTPUT="$OPTARG";; | |
| 70 q) IFPROCSAMPLE="$OPTARG";; | |
| 71 p) IFOUTSAM="$OPTARG";; | |
| 72 | |
| 73 esac | |
| 74 done | |
| 75 | |
| 76 | |
| 77 LOCAL_PATH=`( cd -P $(dirname $0); pwd)` | |
| 78 | |
| 79 if [[ $IFPROCSAMPLE == "Yes" ]] | |
| 80 then | |
| 81 perl $LOCAL_PATH/createControlSubSet.pl -f $CHIPFILE -c $CONTROLFILE -t $TYPE -s $CHIPOUTPUT -o $CONTROLOUTPUT >/dev/null 2>>/dev/null | |
| 82 else | |
| 83 perl $LOCAL_PATH/createControlSubSet.pl -f $CHIPFILE -c $CONTROLFILE -t $TYPE -o $CONTROLOUTPUT >>/dev/null 2>>/dev/null | |
| 84 fi | |
| 85 | |
| 86 if [[ $TYPE == "bam" && $IFOUTSAM == "No" ]] | |
| 87 then | |
| 88 | |
| 89 samtools view -S -b $CONTROLOUTPUT 2>>/dev/null | samtools sort -m 4000000000 - $CONTROLOUTPUT 2>>/dev/null | |
| 90 rm $CONTROLOUTPUT | |
| 91 mv $CONTROLOUTPUT.bam $CONTROLOUTPUT | |
| 92 | |
| 93 if [[ $IFPROCSAMPLE == "Yes" ]] | |
| 94 then | |
| 95 samtools view -S -b $CHIPOUTPUT 2>>/dev/null | samtools sort -m 4000000000 - $CHIPOUTPUT 2>>/dev/null | |
| 96 rm $CHIPOUTPUT | |
| 97 mv $CHIPOUTPUT.bam $CHIPOUTPUT | |
| 98 fi | |
| 99 fi |
