Mercurial > repos > mingchen0919 > aurora_skewer

annotate shell-script-template.sh @ 8:7c9170b7d9f6 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
author mingchen0919
date Thu, 26 Apr 2018 13:38:15 -0400
parents c5c297020c36
children c0e57c994deb
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
5
2c663ddb8488 planemo upload commit f0ae5255b6c75f69c85319bba8da742be83f27d2-dirty
mingchen0919
parents:
diff changeset
1 # run SHELL_SCRIPT within tool outputs directory
2c663ddb8488 planemo upload commit f0ae5255b6c75f69c85319bba8da742be83f27d2-dirty
mingchen0919
parents:
diff changeset
2 cd ${REPORT_FILES_PATH}
2c663ddb8488 planemo upload commit f0ae5255b6c75f69c85319bba8da742be83f27d2-dirty
mingchen0919
parents:
diff changeset
3
8
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
4 #--------- index genome --------
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
5 # create genome directory for genome indexes
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
6 mkdir -p ${X_d}/genomeDir
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
7
5
2c663ddb8488 planemo upload commit f0ae5255b6c75f69c85319bba8da742be83f27d2-dirty
mingchen0919
parents:
diff changeset
8 cat >temp.sh <<EOF
8
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
9 STAR \\
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
10 --runMode genomeGenerate \\
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
11 --genomeDir ${X_d}/genomeDir \\
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
12 --genomeFastaFiles $( echo ${X_A} | sed 's/,/ /g' ) \\
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
13 --sjdbGTFfile ${X_B} \\
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
14 --sjdbOverhang ${X_C} \\
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
15 > /dev/null 2>&1
5
2c663ddb8488 planemo upload commit f0ae5255b6c75f69c85319bba8da742be83f27d2-dirty
mingchen0919
parents:
diff changeset
16 EOF
2c663ddb8488 planemo upload commit f0ae5255b6c75f69c85319bba8da742be83f27d2-dirty
mingchen0919
parents:
diff changeset
17
8
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
18 grep -v None temp.sh > index-genome.sh
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
19
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
20 # run star
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
21 sh index-genome.sh
5
2c663ddb8488 planemo upload commit f0ae5255b6c75f69c85319bba8da742be83f27d2-dirty
mingchen0919
parents:
diff changeset
22
8
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
23 #---- mapping ---------
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
24 cat >temp.sh <<EOF
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
25 STAR \\
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
26 --genomeDir ${X_d}/genomeDir \\
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
27 --readFilesIn \\
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
28 ${X_F} \\
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
29 ${X_R} \\
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
30 > /dev/null 2>&1
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
31 EOF
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
32
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
33 grep -v None temp.sh > mapping.sh
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
34
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
35 # run mapping
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
36 sh mapping.sh
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
37
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
38 # remove temp.sh
5
2c663ddb8488 planemo upload commit f0ae5255b6c75f69c85319bba8da742be83f27d2-dirty
mingchen0919
parents:
diff changeset
39 rm temp.sh
2c663ddb8488 planemo upload commit f0ae5255b6c75f69c85319bba8da742be83f27d2-dirty
mingchen0919
parents:
diff changeset
40
8
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
41 #----- SAM to sorted BAM ------
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
42 echo "samtools sort -o Aligned.out.sorted.bam Aligned.out.sam" > sam2bam.sh
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
43 sh sam2bam.sh
5
2c663ddb8488 planemo upload commit f0ae5255b6c75f69c85319bba8da742be83f27d2-dirty
mingchen0919
parents:
diff changeset
44
8
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
45 #----- evaluate mapping -------
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
46 echo "samtools flagstat Aligned.out.sorted.bam > flagstat.txt" > flagstat.sh
7c9170b7d9f6 planemo upload commit 5b03cf1a31307c262701d4f267972bd01be24311
mingchen0919
parents: 7
diff changeset
47 sh flagstat.sh