annotate mytools/makebigwig.sh @ 0:39217fa39ff2

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