annotate tools/mytools/makebigwig.sh @ 0:9071e359b9a3

Uploaded
author xuebing
date Fri, 09 Mar 2012 19:37:19 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
1 # use of output: move to public_html and visualize in ucsc genome browser with the following:
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
2 # track name="xxx" color=0,0,255 type=bigWig bigDataUrl=http://rous.mit.edu/~wuxbl/xxx.bw
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
3
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
4 if [ $# -lt 6 ]
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
5 then
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
6 echo "./bigwig.sh infile outtag bam/bed sorted/none genome strand/none [-split]"
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
7 exit
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
8 fi
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
9
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
10 f=$1
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
11 outf=$2
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
12 extension=$3
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
13 sorted=$4
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
14 genome=$5
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
15 strand=$6
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
16 split=$7
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
17 i=i
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
18 if [ $extension = bam ]
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
19 then
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
20 i=ibam
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
21 if [ $sorted != sorted ]
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
22 then
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
23 echo 'sorting bam file...=>' $f.sorted.bam
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
24 samtools sort $f $f.sorted
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
25 f=$f.sorted.bam
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
26 fi
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
27 else
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
28 if [ $sorted != sorted ]
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
29 then
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
30 echo 'sorting bed file...=>' $f.sorted.bed
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
31 sort -k1,1 $f > $f.sorted.bed
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
32 f=$f.sorted.bed
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
33 fi
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
34 fi
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
35
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
36 echo 'making bedgraph file...=>' $f.bedgraph
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
37 if [ $strand != strand ]
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
38 then
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
39 genomeCoverageBed -bg -$i $f -g $genome $split > $f.bedgraph
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
40 echo 'making bigwig file...=>' $outf.bw
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
41 bedGraphToBigWig $f.bedgraph $genome $outf
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
42 else
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
43 genomeCoverageBed -bg -$i $f -g $genome $split -strand + > $f+.bedgraph
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
44 genomeCoverageBed -bg -$i $f -g $genome $split -strand - > $f-.bedgraph
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
45 echo 'making bigwig file for + strand...=>' $outf+.bw
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
46 bedGraphToBigWig $f+.bedgraph $genome $outf+.bw
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
47 echo 'making bigwig file for - strand...=>' $outf-.bw
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
48 bedGraphToBigWig $f-.bedgraph $genome $outf-.bw
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
49 fi
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
50
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
51 # remove intermediate files
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
52 if [ $sorted != sorted ]
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
53 then
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
54 rm $f
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
55 fi
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
56 rm $f*.bedgraph
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
57